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ROAD MAP 


REVISIONS AND ADDITIONS TO FEBRUARY SUBMITTAL, 

Volume 2,ERTS SYSTEM STUDIES, initially issued m February, 
presented the results of two study areas 

• Booster vehicle choice 

• Orbital analysis 

As stated m the introduction and summary of that volume, the 
results of two additional areas 

• Improving GDHS performance through observatory 
modifications'’'- 

9 Total system accuracy 
were to be presented m April. 

This volume presents m detail the results of these two additional 
studies In addition, the introduction and summary of the previous 
volume has been expanded to include the new area of study Furthermore, 
typographical errors of the previous volume have been corrected. 


In that volume this study was identified as Potential for improving 
GDHS performance through sensor or spacecraft modifications 
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1. INTRODUCTION AND SUMMARY 

The results of those Phase B/C studies of the ERTS system that 
encompass elements of both the observatory and the ground data handling 
system or that can affect factors in the program beyond these two systems 
are reported in this volume. Four specific studies are involved 

• Booster vehicle choice 

• Orbital analysis 

• Improving GDHS performance through observatory 
modifications 

® Total System Accuracy 

Section 2 contains a detailed examination of the relative advantages 
of the Thor-Delta and Thor-Agena launch vehicles It is concluded that 
they are both equally acceptable. Included m the analysis were weight- 
carrying capability, injection accuracy, vibration loads, launch facilities, 
and shroud Because no compatibility -as so ciated costs arose in this 
analysis, vehicle costs are essentially those chargeable to any other 
program Thus the Thor-Delta as recommended by GSFC m the design 
study specification is acceptable. 

Orbital analysis is the subject of Section 3. Orbital parameters 
for the ERTS mission have been refined to ensure a repetitive ground 
trace every 18 days The resultant sun-synchronous circular orbit has 
the following basic elements 

• Semi-major axis, a = 3936 5 n mi 

• Right ascension of ascending node at vernal equinox, 

= 142. 5 deg 

• Inclination, l = 99. 0848 deg. 

Analyses have been performed on the effects of luni-solar, geopotential, 
and atmospheric-drag perturbations on the orbit and consequently .the 
resultant images Grouping the results of these perturbations into short 
periodic (equal to or less than an orbital period), long periodic, and 
secular, it is found that ..while they are appreciable, all periodic effects 
can be compensated A secular effect due to the gravity gradient of the 
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sun acting on the orbit plane is found to induce a growing inclination 
error and consequently deviation from sun -synchronism Approximately 

17 feet per second of velocity correction is required per year to offset 
this effect* However, by initially biasing the inclination, a year of tol- 
erable operation after initial correction is possible without velocity trim 
An additional fraction of a foot per second per year of velocity correction 
is required to remove the cumulative effect of atmospheric drag 

Ground lighting has been calculated for a nominal launch and for one 
that slips one-half hour* While it was our intention at the start of the 
study to perform calculations for an orbit slightly off from sun-sychronism, 
improvements m the launch vehicle injection accuracy have been so great 
that it is unnecessary to be concerned with drifting from sun -synchronism 
as the result of booster injection errors 

Computer runs have been made for orbital conditions arising from 
launches m the springs of 1972 and 1973 Drag models, lunar and solar 
attractions, and a suitable geopotential model were included m the pro- 
gram* The results lend credence to the analyses and affirm that all 
orbital requirements of the ERTS mission can be met 

Section 4 is concerned with improving GDHS performance through 
observatory modifications. Specific subjects - that have been examined 
during the study period are on-board light source with precisely known 
intensity levels that can be turned on and off by commands m order to 
obtain accurate photometric calibration of RBV images, projection of 
variable intensity periodic reseau patterns to obtain improved geometric 
and photometric accuracy over the opaque cross reseaus, incorporation 
of RBV shutter time in narrowband PCM for increased flexibility m the 
time line operations of the GDHS, and improved attitude determination 
with a two gyro system. 

A synopsis of the first three study areas is given A more detailed 
examination of these areas is found m Volume 17 The remainder of 
Section 4 is concerned with attitude determination earth horizon surface 
model, system dynamical equations, and two gyro performance analysis. 

The roll and pitch attitude of the spacecraft will be obtained from 
earth horizon scanners In order to achieve the 2 nautical mile accuracy 
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using horizon scanners, a software program to periodically update the 
radiance model will be required Included m Section 4 are the derivation 
of the equations which will form the basis for the operational software 
for the updating of the radiance Model 

To obtain highly accurate attitude rate information using either a 
one or two rate gyro configuration, ground data processing of the horizon 
scanner, rate gyro, and reaction wheel data are required A starting 
point in deriving the filtering equations for processing the attitude and 
attitude rate information is to start with the equation of motions govern- 
ing the attitude of the spacecraft Also the relationship between the 
observed data and the elements of the attitude state vector must be de- 
termined Developed m detail are the dynamical and measurement 
equations 

Considered next m this volume is an analysis of the horizon scanner- 
two gyro system For this system the two mile objective is easily sat- 
isfied while the internal consistency error is 600 feet. 

Section 5 is concerned with the total mapping accuracy and photo- 
metric performance for the selected bulk and precision processing modes 
- Bulk I, Bulk II, and Precision I. The mapping performance, both 
absolute and internal consistency, are examined separately for the RBY 
and MSS images because of the significance differences m their sensitivity 
to attitude errors and m the methods of image processing 

For all image processing modes, the LBR is the first image pro- 
cessing black box through which both the RBV and MSS data flow. The 
mapping performance, as is done here, can therefore be conveniently 
examined m terms of pre-LBR and post- LBR 

MSS pre-LBR mapping accuracy is first examined All significant 
error sources are identified, e g , attiude, ephemeris, sensor Given 
the major error sources, the sensitivity matrix of the cross track and 
down track errors to each of the error sources is determined numerically 
The sensitivity matrix is used m a covariance analysis to determine 
absolute location performance and internal consistency performance in 
which various time and spatial correlation of errors are considered 
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Next, the RBV pre-BR mapping accuracy is considered. Here, 
the analysis is accomplished m two steps. Pre-reseau and reseau 
estimation accuracy The pre-reseau analysis is accomplished, as m 
the MSS situation by the utilization of the appropriate sensitivity matrix 
The reseau estimation accuracy analysis determines the mapping errors 
due to least squares fitting a biquadratic to a 3 x 3 set of reseau points 
The total pre-CBR RBV error is the rss value obtained from the two 
analyses To obtain the post-L.BR performance, the error sources and 
an explanation of how these errors arise for two types of LBR's are 
given A complete tabulation of the magnitudes of the significant sources 

of errors is also given. Based on these values, the mapping errors due 
to an LBR are then determined The L.BR errors are independent of the 

errors previously considered m the MSS and RBV analyses The overall 
geometric errors is therefore easily obtained by the rss method 

The section closes with a detailed MTF and photometric analysis 
For the MTF analysis, seven MTF are identified and evaluated This is 
followed by the photometric analysis. The results of the analysis was 
generated by utilizing a Monte Carlo computer program m which a large 
number of photons were generated, their paths traced through the 
atmosphere, reflected or observed at the surface, and the returns traced 
The output is the fractional photon returned 
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2 BOOSTER VEHICLE CHOICE 


Whereas m the spring of 1969 some advantage might have been 
ascribed to the use of a Thor-Agena, this is no longer the case. The Thor- 
Delta and Thor-Agena are equally acceptable technically for ERTS-A and -B 
Since we have found no technical reason why the Thor -Delta launch vehicle 
should not be used, and since the space craft- induced costs are approxi- 
mately the same for Thor-Delta and Thor-Agena, we have made no cost 
comparison between Agena and Delta stages on a Thor 

Three OGO's have been successfully launched by Thor-Agena, 
assuring the acceptability of Thor-Agena as a launch vehicle The prob- 
lem of this study then was to ascertain what differences exist between 
Thor-Agena and Thor-Delta that could affect the ERTS. Data on vehicle 
performance, physical and electrical interfaces, boost environments, and 
launch complex have been obtained from McDonnell-Douglas and Lockheed 
and reviewed in the light of ERTS requirements 

2 t LIFT 

A dramatic improvement has been announced m the performance of 
the Thor-Delta planned for 1972 launches, resulting from a new inertial 
guidance system and a new Delta stage engine (N~>0^) designated DSV- 
3L/1D (^O^) McDonnell-Douglas estimates a lifting capability of 1600 
pounds into an ERTS orbit NASA/GSFC has established a more exact 
definition of spacecraft weight -lifting capabilities as shown below 

Spacecraft Weight Thor-Vehicle 

(pounds) Solid Rockets 

1480 3 

1900 6 

2100 9 

The three -solid version of the Delta is most nearly matched to the ERTS 
need, which is about 1400 pounds for ERTS-A and 1425 for ERTS-B 

2 2 INJECTION ERROR 

The Delta improvement program includes the development of a 
strapped-down inertial system based on inertial measurement units of the 
Apollo abort guidance section, providing a significant improvement in 
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guidance accuracy Delta engineers are confident that the improved per- 
formance will be available m the ERTS 1972 launch. The values listed 
below are supplied by NASA/GSFC from a McDonnell-Douglas orbital 
study of the ERTS mission 

99 Percent Probable 
Minimum Maximum 


Apogee altitude deviation (n mi) 

-3 

11 

Perigee altitude deviation (n mi) 

-11 

3 

Orbit period deviation (mm) 

-0.3 

0.3 

Orbit eccentricity deviation 

0.000 

0.001 

Orbit inclination deviation (deg) 

-0.03 

0. 03 


The relationship of these errors is also established m that study. Our 
estimate of ERTS orbit adjustment requirements based on the above is 
given in Section 3. 4. The conclusion is that orbit error can be removed 
satisfactorily with a relatively simple pneumatics system (see Volume 3, 
Section 10). 



DELTA 5TA 
724 3 
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Figure 2-1 

THE OBSERVATORY FITS wrthm the 



2.3 B OOST ENVIRONMENT 


At the time of writing our June 1969 ERTS proposal we understood 
the Thor -Delta vibration environment was somewhat more severe than 
the Thor-Agena. In fact, as established from data from the launch of 
OGO 6 and other sources, the vibration environment imparted to the 
spacecraft is approximately the same, stemming m both cases principally 
from the burn of the Thor during liftoff and transonic flight Thus we can 
establish no sensible difference between Thor-Agena and Thor-Delta 
spacecraft vibration levels, and we conclude that the ERTS can survive 
the environment provided on both boost vehicles. 

2.4 SHROUD INTERFACE 

The available Thor-Delta shroud ientified as Figure II of the Study 
Specification and also the DACS s hr our are both compatible with our pro- 
posed ERTS. The specification shroud is identified by Me Donne 11 -Douglas 
engineers as the Nimbus shroud, and the design has been used for OGO's 
1 and 3. Figure 2-1 shows ERTS m the Thor-Delta shroud. No changes 

have been made m ERTS that affect shroud compatibility Array corners 
are unchanged from OGO. a urmer interference points of concern, the 
experiment packages on OGO, have been eliminated 




Thor-Delta shroud with adequate margins 


2-3 



2. 5 RING-INTERSTAGE INTERFACE 

The OGO-Agena interface was in the form of a two -inch -thick ring 
60 inches in diameter The Delta ring defined m the Study Specification 
is two inches thick and 60.5 to 59.875 inches m diameter. To adapt the 
OGO interstage to the McDonnell-Douglas ring, all that is needed is to 
adjust the bolt holes slightly (Figure 2-2). 

2.6 ELECTRICAL INTERFACE 

The Delta umbilical disconnect is quite different from the Agena 
design For ERTS purposes it is superior since we have found it desir- 
able to install equipment outside the +Y end of the spacecraft where the 
Agena spacecraft disconnect was located on OGO. 

Delta places a 30 -wire limit on the umbilical cable, which is accep- 
table for ERTS In addition to the umbilical, ERTS requires two coaxial 
cables to be brought into the spacecraft while it is on the launch stand 
These can be disconnected manually before stand removal McDonnell- 
Douglas states that these cables can be readily connected through a shroud 
access hole. The coaxial lines would carry PCM data and unified S-band 
baseband to permit tests during radio silence at Vandenberg Air Force 



Figure 2-1 (concluded) 
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Base. The Agena design provides umbilical disconnect of coaxial lines, 
although OGO has used only one such circuit m the past. Other electrical 
changes m launcher, 450-foot tower, and assembly building are expected, 
but since these will also be needed for Agena launches they are not rele- 
vant to this comparison. 

2.7 LAUNCH STAND 

McDonnell-Douglas reports that launch stand SLG-2W will be used 
for NASA Delta launches in 1972. There is no aspect of the physical 
characteristic of this stand which would impede work m installing or 
servicing ERTS. 



Figure 2-2 

INTERSTAGE RING between Thor-Delta and ERTS is the same 
as that between OGO and Thor-Agena except For the location 
of the bolt holes 
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3. ORBITAL ANALYSIS 


The goals of the orbital analysis have been to 

• Establish basic orbit elements 

• Find the effects of perturbations on the nominal orbit 

• Determine image distortion resulting from the orbit 

• Analyze sub satellite illumination conditions 

• Verify illumination angle at the satellite 

• Determine orbit injection errors and define strategem 
for their removal 

• Analyze tracking coverage. 

In the Phase B/C proposal, TRW indicated the extent to which 
ERTS orbital studies had progressed by the spring of 1969. In that pro- 
posal some advantages of making the orbit ground trace repeat exactly 
on a 17- or 18-day cycle were indicated. Since the 18-day repetition is 
now a requirement for the ERTS program, certain parameters associated 
with the mission are altered and it has been necessary to repeat some 
earlier analyses first to obtain these new parameters and then to use 
them to revise the perturbation analyses, ground lighting conditions, and 
coverage calculations. 

3. 1 BASIC ORBITAL ELEMENTS 

The ERTS orbit is nominally circular and sun -synchronous. Sun- 
synchronism exists when the right ascension of the ascending node, £2, 
advances 360 degrees m one tropical year, or £2 = 0.985, 569,42 deg/day. 
This will be true for the circular orbit when the semi-major axis, a, and 
the inclination, l, relate in the following way 

£2 = 0. 985, 569 = -9. 96468 (a^aj^^ cos i (3-1) 

where a . is the earth's equatorial radius. 

ERTS orbits the earth not quite 14 times a day such that after 
18 days it repeats the ground trace. The rate of advance of the longitude 
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of the ascending node of the mean satellite, is computed as follows 5 *' 


\ N = (M+ »)/S Q - (6 - £2) 


(3-2) 


where M is the mean anomaly of the orbit, oj is the argument of perigee, 
S q = 14 is the nominal number of orbits per day, 0 is the right ascension 
of Greenwich, and £2 is the right ascension of the ascending node. 

In 18 days the observatory makes exactly 18 x 14 -1 orbits in repeat- 
ing the ground trace. Thus the change AX-^j m X.^ m 18 days is 


From this, 


= 360 x 


251/S 



-360/ 14 degrees 


= AK N^ 18 = ” 1 * 42857 deg/day 


was used as an input to TRW’s Rapid Orbit Prediction Program 
(ROPP). ROPP adjusted automatically the nominal semi-major axis to 
yield the required X.^ value The result was 

a = 3936. 5 n mi 

or expressed as altitude of the circular orbit above the equator, 

h — 492. 6 n mi 1 " 1 

Using this value for a. Equation (3-1) gives 

1 = 99 085 degrees 

For a launch that does not slip, the local time of the descending 
node is 9 30 AM This corresponds to a right ascension of the ascending 
node at vernal equinox, of 142 5 degrees 


G S. Gedeon, ’’Tesseral Resonance Effects on Satellite Orbits, 11 
Celestial Mechanics 1, No 2, 167 (1969) 

Computed m terms of SFADATS variables 
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3 2 EFFECTS OF PERTURBATIONS 


The orbital elements established m the previous sections are in 
reality only approximations to actual orbital motions. The nonsphericity 
of the earth, the gravitation of the sun and moon, and atmospheric drag 
all tend to alter the two-body orbital motion. Orbital perturbations 
manifest in three ways as short periodic (equal to or less than an 
orbital period), long periodic, and secular. The effect of these pertur- 
bations is to alter the position and velocity of the spacecraft in orbit and 
thus to cause size and centering shifts m both RBY and MSS images and 
anamorphic stretch m MSS pictures. Analysis of these perturbations 
follows and their effect on image quality is discussed in Section 3 4. 

3 2 1 Short Periodic Perturbations 

Short periodic perturbations are predominantly caused by since 
it is a thousand times larger than the next largest harmonic in the 
potential field Thus m this analysis only is used. 

For circular orbits these perturbation expressions become 

Aa = (3/2) a (a /a )^ J 7 sm^i cos 2 (M + co) 

© « 

Ae = 0 

2 

Ai = - (3/4) (a e /a) sin i cos i cos 2 (M + to) 

An = (a g /a)^ cos l sin 2 (M + w) 

Aco = \ (a /a)^ J_ [l + (3/2) sin^ i] sin 2 (M + co) 

7 2 

AM = - 6 (a /a) J, sin i sin 2 (M + w) 

© L 


G S Gedeon, "Analytic Partials, " TRW memorandum 3412.4 J84 
(September 1968). 
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Substituting the values of Section 3. 1, 


Aa = 4. 75 cos 2 (M + co) n mi 
Ae = 0 

Ai = - 0. 0055 cos 2 (M + w) degrees 

AJ2 = 0. 0057 sm 2 (M + co) degrees 

AtD = 0. 086 sm 2 (M + co) degrees 

AM - -0. 277 sin 2 (M + co) degrees 

These short periodic perturbations (with twice the orbital frequency) 
due to J£ produce a variation m the semi -major axis with an amplitude of 
4.75 n mi, and an along-track oscillation resulting from the combined Aw 
and AM effects with an amplitude of 0. 191 degree. The semi -major axis 
is maximum above the equator, along-track variation is maximum at the 
apex of the orbit (maximum latitudes). Small amplitude across- track 
oscillations arise from the Ai and A S2 terms. 

3.2.2 Long Periodic Perturbations 

The term also produces a long periodic rotation of the line of 
apsides. Apsidal rate is calculated from 

A = -4 98234 (a e /a) 7/2 (i - 5 cos 2 i) 

When the appropriate values of a and l are substituted, w is found to 
equal 2.75 deg/day. Thus it takes 131 days to complete an apsidal cycle. 

This apsidal motion interacts with an eccentricity variation 

introduced by the J, term (pear-shaped earth) to produce long periodic 

$ 

(131 days) variations m the eccentricity Kaula's Equation 3. 77 can be 
put in the following form 

Ae = (i/2)(J 3 /J2)( a e / a ) sin (cj q + 4> t) 


' W. M Kaula, Theory of Satellite Geodesy (Blaisdell Publishing Co. , 
Waltham, Massachusetts, 1966). 
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Substituting the values of a and <o gives 

Ae = 0 001025 sm (a> 0 + 2 75t) 

Since the observatory makes almost 14 revolutions per day, it is 
subjected to resonance with the fourteenth-order tesseral harmonics 
Lagrange's planetary equations are used to evaluate the perturbations 
due to the tes serais The rate of change of the semi-major axis is 


da/dt = (2/na)(SF/SM) 


where F is the perturbing function and n = 5020 deg/day is the mean 
motion. It is difficult to take the partials of the potential function with 
respect to Kepler elements if the potential function is expressed by the 
standard associated Legendre polynomials. Kaula, however, made this 
task very simple. He rotated the reference plane from the equator to an 
arbitrary orbital inclination and then in this plane executed a Fourier 
senes expansion to introduce the mean anomaly (i.e., the time). The 
result of this transformation is a potential function expressed by the 
orbital elements 


V = |i/r + 


CO A 

£ E 

A = 2 m= o 


£ S V. 
p = o q=°o 


where 


t 


V 


fmpq 


-£( i) 

a \ a ' 


F #mpq ^ G fpq ^ ^ 


cos v 
sm 


(1 - m) even 
(1 - m) odd 


and 


{ (i - 2p) u + (A - 2p + q) M + m (Q - 0 - X ) \ 

Am 

In this expression •^ nl pq( 1 } lS a function depending only on the in- 
clination, and G^p^(e) is a function depending on the eccentricity alone. 
These functions are given by Kaula in the forms of both equations and 
tables 
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Though the above expression looks quite complicated, it has two 
important advantages The first is that partials needed m Lagrange 1 s 
planetary equations can be immediately taken. The second is that the 
disturbing function can be specified to be secular, long periodic, 

or short periodic by simply specifying the indices H, m, p, q, e. g. , 


V 

2010 

secular 

V V 

3021, v 303-1 

long periodic 

^ 20 pq 

(p and q / 1 , 0 ) short periodic 


The fourteenth- order tesserals produce the following long periodic 
perturbing functions 


V 15, 14, 7, 0 V 17, 14, 8,0 V 19, 14, 9, 0 * * * 


All these can be verified by substituting the values of J2mpq into 4*. 

For the sake of simplicity we restrict our analytic investigation to the 
V 15 j q Fourier component, for which* 

-22 

Jig = 0. 69 x 10 and \^g = It. 3° 

Now Lagrange's planetary equation can be integrated with 
to the first order if we assume that on the right-hand side the variations 
are due to the dominant V^q ^ q term which produces the well-known secular 
perturbations of the angular elements M, w, The result is 


a = a + — — 2p * q V. 

o na tmpq 


R. J Anderle, "Observations of Resonance Effects on Satellite Orbits 
Arising from the Thirteenth and Fourteenth-Order Tes serai Gravitational 
Coefficients, " Journal of Geophysical Research 7 0, 2453 (May 15, 1965) 
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Z 3 

Substituting 7 q and observing that p. = n a C 

a = a o + (Zna/«i») (*>) 15 U 7 (i) G 15j7> Q (e) J ^ u am * 

The argument + can be expressed as 

* 3 m(K N - W ~ ** 

where as defined earlier is the longitude of the ascending node of the 

"mean" satellite. In the case of circular orbits, q = 0 and G„ =1 The 

fpq 

inclination function for 99 085 degrees computed by a special computer 
program is- 

F 15, 14, 7 085 °) = 1 95 x 10 14 

Since \ N = -i. 42857 deg/day from the required ground trace shift, 
ip = mX. n =-20 deg/day 

and the period of this perturbation is 360 /«|j= 18 days. The variation of 
the semi-major axis becomes A a = 22 5 feet 

Because of the cyclic variation of the semi-major axis there will 
be a periodic variation of the ground trace drift rate This can be 
obtained r by taking the derivative of Equation ( 3 - 2 ) 

= (l/S o ) (M + w) - (6 - S2) 

The amplitude of this variation is 

6 K N = (1/S Q ) 6 M = - ( 3 / 2 ) (n/S o ) (Sa/a) 

For the ERTS orbit, &V N = 0 0005 deg/day 
Gedeon, op cit 
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The period of this variation is 18 days. Thus the amplitude of the 
variations of X.^ is obtained by integration as 

A\ n = (P/2ir) 6 X = 0. 00145° 

This, then, is the amphtude of the periodic deviation from the nominal 
ground trace shift 

Resonance also causes an along-track perturbation which can be 
calculated from 1 ' 

F imp ( l > J im f 2(i +l)-3Ejs i „ + 

Substituting numerical values, A= 0 0186 sm $ This corresponds to a 
1. 28 n mi alOng-track oscillation, or 320-millisecond station error 

3. 2. 3 Secular Perturbations 

Secular perturbations of the line of nodes and line of apsides are 
represented by 


£2 = 


-9.96468 (a e /a} 7/2 cos i 


and 


01 = -4.98234 (a e /a) 7 ^ 2 (1-5 cos i) 

Secular perturbations of the semi-major axis are introduced by drag 
A coefficient of drag for the spacecraft was derived using free molecular 
flow theory and a coefficient of normal and tangential momentum exchange 
of 0. 8 This yields = 2. 4. The frontal area of the spacecraft was 
taken to be A = 20 sq ft and the spacecraft weight 1400 pounds. These lead 
to a ballistic coefficient of 


C d A/W = 0.0343 


^Gedeon, op cit 
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Atmospheric density is determined using the TJ S 1966 Standard 
Atmosphere and setting the solar activity index, F ^ - 100, for a 
March 1972 launch and F^g -, = 88 for a March 1973 launch A geomag- 
netic planetary index of 2 0 was used for both years. Although approxi- 
mate analytic results of atmospheric drag have been calculated, the 
complexity of the density model indicates that computer runs are more 
reliable Additionally, m the computer runs luni- solar perturbations 
can also be included. 

Because the orbit is nominally sun- synchronous, a gravity gradient 
effect from the sun induces a change m orbital inclination and consequently 
a drift from sun- synchronism 

3 2. 4 Computer Runs 

To integrate the orbit m the presence of the variety of perturba- 
tions discussed m the previous section, orbital parameters as generated 
m Section 3. 1 were used Additionally, the eccentricity {which is 
nominally zero) was set at e = 0 001 to more nearly characterize launch 
and velocity correction residuals The initial location of perigee was at 
the ascending node (to Q = 0) Lunar and solar perturbations were 
included as were J 2 , Jy 3^, J 15 J 17 and The 

results of the computer run i r 1972 and for 1973 are presented in 
Figures 3-1, 3-2, and 3-3. 

Figures 3-1 (a) and 3-1 (b) show the semi-major axis and the 
eccentricity variations for a March 22, 1972 and a March 22, 197 3 launch 
The semi-major axis decay is about 136 ft/year for the 1972 trajectory and 
111 ft/year for the 1973 trajectory Resonance causes an oscillation with 
a 38-foot amplitude and 18-day period. The amplitude obtained by ROPP 
is higher than the analytic value, since two more harmonics were included 
m the ROPP run. 

Eccentricity variation has an amplitude of 0 00102 and a period of 
132 days, which agree very well with the analytic values presented m 
Section 3 2 

Figure 3-2 shows the variation of the mean node in 18-day intervals. 
The maximum deviation is about ±3 5 n mi This graph was obtained aftei 
several corrections were made to the initial semi-major axis and to the 
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Figure 3—1 (a) 

SEMI-MAJOR AXIS AND ECCENTRICITY AS FUNCTIONS OF TIME, 
March 22, 1972 launch 





ViA RATION (FEET) 



MARCH 22, 1973 T |M E (DAYS) 


Figure 3-1 (b) 

SEMI-MAJOR AXIS AND ECCENTRICITY AS FUNCTIONS OF TIME, 
March 22 , 1973 launch 
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Figure 3-2 

VARIATION OF THE MEAN NODE 


initial inclination to offset the effect of a strong solar perturbation which 
affects the rate of change of the right ascension of the node, Q 

The gradient in the sun' s gravity field exerts a net torque on the 
spacecraft about the earth’s center With the inclination and right 
ascension of the node that prevail for ERTS, this torque can be shown to 
be nearly parallel to the earth's spin axis This torque causes a change 
in the angular momentum Since this change is northward, the orbital 
inclination must decrease, as shown in Figure 3-3. 

Now the earth's bulge causes a precession in the line of nodes that 
is proportional to the cosine of the inclination. Hence as the sun reduces 
the inclination from the initial value, the cosine diminishes in magnitude 
and consequently the nodal velocity diminishes This also can be seen m 
Figure 3-3 (Note the "beat" which is due to the 18 period resonance and 
the ~15 day period lunar perturbation ) 

The effect of solar perturbations is to displace the line of nodes by 
0. 9 degree or 54 n mi after one year, m a westwardly direction For- 
tunately the effect of drag is to partly counteract this, since orbital decay 
by drag causes the nodal crossing to move eastward However, the 
effect is only half as large as that due to the sun The decay during one 
year was about 136 feet, producing an eastward shift m the equatorial 


3-12 



crossing (in one year) of 27 n mi Hence due to solar perturbations and 
drag the position error at the equatorial crossing after one year is 
27 n mi (westward) This effect can be largely negated by starting the 
orbit at a slightly higher inclination, 0 013 degree higher than required 
for sun- synchronism, i e , at 99 098 degrees Then the resulting 
excursion m the nodal crossing is only 7 n mi Finally, by starting the 
orbit 3. 5 n mi westward from the nominal (at the first equatorial 
crossing), the maximum excursion m the equatorial crossing from a 
nominal one will be only ±3. 5 n mi This is shown m Figure 3-2, 

3 3 IMAGE DISTORTION RESULTING FROM ORBIT PERTURBATIONS 

Along- and across-track orbital perturbations result m centering 
shifts of images. Radial changes in orbital position as well as the 
earth's oblateness lead to image size errors. 



TIME (DAYS) 

Figure 3-3 

VARIATION OF THE INCLINATION AND THE RATE OF CHANGE 
OF THE LINE OF THE NODES 
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INCLINATION (DEGREES) 



For nearly circular orbits we can write 


r = a (1 + e cos v) 

by expanding the conic equation to the first order m e. 

From Section 3. 2, the semi-major axis variation due to oblateness 
is 


a = a^ + Aa cos 2 (M + w) 

where Aa = 4 75 n mi. Substituting this m the first equation, and neglect- 
ing higher order terms, we get 

r ^ a^ + Aa cos 2 (M + o>) - a^e cos v 

A cross section of the oblate earth which is inclined by 99 degrees 
to the equatorial plane can be written as* 

R = [3444 2 cos 2 6 + (3444 - ll) 2 sm 2 0] 1/2 

where 0 is a central angle measured from the equatorial plane which 
also equals (M + co). Expanding R to the first power 

R = 3438.5 + 5.5 cos 20 
Then the geodetic height is 

h = (a^ - 3438 5) + (Aa - 5. 5) cos 20 - a^e cos v 

We have already seen that the harmonic causes a periodic variation 
of the eccentricity with 

Ae = 0. 001 sm w 

Then noting thatv = 0 - «, the geodetic height with the ERTS data becomes 
h = 498 - 0 75 cos 20 - 3 94 cos (0 - o>) 
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Over the ranges of 9 and co, the largest negative variations are 

Equator 45° Latitude Near the Pole 

~4 7 -3 94 -3 20 

Figure 3-4 shows the variation of the altitude above the reference 
ellipsoid and that above a sphere with radius equal to the equatorial 
radius of the earth. The graph was obtained by integrating an orbit with 
all perturbations present, but with an eccentricity of only 0. 0001. 

Given a side-to-side RBY field of view of 11. 5 degrees, the maxi- 
mum altitude variation causes a variation of 0 94 n mi m picture siae 
The short periodic along-track oscillations due to and e cause an image 
slippage which can be obtained from the derivative of 

2 e sin (M + u) - (0. 1 91 /57 3) sin 2 (M + co) 



TIME FROM DESCENDING NODAL CROSSING (MINUTES) 


494 

493 

492 

491 

490 


Figure 3-4 

VARIATION OF ALTITUDE AND RADIUS DURING ONE ORBIT 
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RADIUS MINUS EARTH'S EQUATORIAL RADIUS (N Ml) 



when e = 0.002, which amounts to a maximum of 0.0000108 rad/sec. 

Since RBV images are snapped every 25 seconds, the slippage can amount 
to 0. 0000108 x 25 x 3444 = 0. 935 n mi. For the MSS the stretch or com- 
pression m one picture 100 miles square could be as large as 1.03 n mi. 

Finally, the side overlap variation due to resonance is (0. 0024/57. 3) 
x 492. 6 = 0. 02 n mi. 

The combined secular effect of drag and luni- solar gravity is 
shown m Section 3 2 to lead to a ±3. 5 mile cross-track variation in the 
course of a year. Should it be desirable to reduce this to a smaller 
value or should the mission from all other considerations endure beyond 
a year, an out-of plane and an in-plane correction would be required to 
restore the orbit to its original state. The magnitude of the out-of- 
plane correction would be approximately 17 ft/ sec per year and the 
magnitude of the in-plane correction would be about 0 14 ft/sec per year. 

3 4 ORBIT INJECTION ERRORS 

The ERTS mission requires that an orbit of very tight tolerance be 
established and maintained It is required that the orbit be circular as 
nearly as possible, with correct semi -major axis, and that orbital incli- 
nation be such as to result m sun- synchronism; The semi-major axis is 
determined by requiring a repeating ground trace after exactly 18 (solar) 
days. The sun- synchronism requirement demands a retrograde orbit 
with an inclination of 99. 098 degrees, and the repeating ground trace 
requirement then demands that the mean semi-major axis be 3936 5 n mi. 
The eccentricity is established by requiring a "minimum altitude varia- 
tion orbit. 11 The mean eccentricity is therefore small, but is not quite 
zero because of the asphencity of the earth 

Orbital tolerances are determined by the desire to obtain complete 
photographic coverage without any gaps between adjacent picture frames 
This requires that the eccentricity be less than 0. 003 and that the semi- 
major axis be correct within 20 feet in order that the track after one year 
comes within 10 n mi of the first track of day zero. 

The tight orbital tolerances cannot be established at injection with 
any of the guidance systems considered, and it is therefore necessary to 
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perform a series of orbital correction maneuvers by the spacecraft Until 
quite recently two different guidance systems were m principal contention 
The first with relatively large injection errors, is the Western Electric 
Company guidance system (WECO) This is an open-loop radio command 
system The second contender (and that adopted for this study) is the 
Delta Inertial Guidance System (DIGS) As the name implies, it is a closed- 
loop, on-board inertial system, with errors considerably below the WECO 
system 

Although the DIGS system will be used on ERTS, previous analyses 
of the WECO system are also reviewed here Upon analyzing the latest 
WECO covariance matrix it was found that m approximately 33 percent of 
the cases the total velocity required for removing all the significant injec- 
tion errors ( 1 . e , errors m apogee, perigee, and inclination) exceeded a 
nominal velocity correction capability of 100 ft /sec. Therefore, in about 
33 percent of the cases the desired orbit could not be established In such 
cases it was assumed that the least critical of the ERTS requirements was 
the sun- synchronism Consequently, where all the errors could not be 
removed, a residual error m the inclination was allowed However, m 
such cases the targeted value nr the semi-major axis must be slightly 
changed or a repeating grout trace is not obtained To illustrate, with 
a spacecraft velocity correction capability of 100 ft/sec the expected value 
in the (absolute) residual inclination error using the WECO system is about 
0 03 degree. This transfers into an error m the nodal velocity of 0 003 
degree per day If the orbit were targeted to the nominal semi-major 
axis, the error in the nodal crossing after 100 days would be 0 3 degree, 
or about 18 n mi, an intolerable error. To prevent this, the orbital period 
would have to be changed by 0. 05 second. The required change in the semi- 
major axis is 130 feet, and the required velocity necessary to change by 
this amount is 0. 06 ft /sec, an insignificant amount from the standpoint of 
total velocity correction requirements Notice that with this scheme the 
ground tract would remain periodic (repeating) but the period of coverage 
cycles would change by 1 2 minutes, i e , would be 18 days ±12 minutes 
The orbit would therefore deviate from sun-synchronism by this amount 
every 18 days or 24 3 minutes per year 
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The critical orbital parameters for the ERTS mission are a, e, and 

1 Instead of a and e, however, a more fruitful way of looking at the 

problem is to say that we must target for specified values of i, R , and 

a 

The last two terms denote the radius at apogee and perigee, respec- 
tively For a given set of injection conditions the purpose of the correc- 
tion maneuver is to remove the errors m 1 , R , and R To establish 

a p 

the minimum velocity correction required for the removal of these errors 

is a difficult problem since a nonplanar transfer must be considered In 

general the minimum required correction velocity depends not only on 

Ai, AR , and AR (the errors in 1 , R , and R ), but also on the argument 
a p a p 

of perigee* However, for the ERTS mission other constraints enter in, 
to rule out certain types of transfers. Thus the desire to exclude com-* 
plicated attitude maneuvers has led to the decision to use separate 
thrusters for the m-plane and out-of-plane corrections Effectively, this 
means that the xn-plane and out-of-plane corrections decouple, and the 
total velocity correction requirement is the sum of the two Of course 
a velocity correction penalty must be paid which in some instances is as 
high as 50 percent over the velocity requirements if a coordinated 
maneuver were made* 

The optimum m-plane transfer is given by Hohmann's bi-tangential 
ellipse The required total m-plane transfer velocity is given by 

AV 1 = (V / 4a) (AR a + AR p ) 

where a is the semi-major axis and Y = (p/a) is the constant orbital 
speed (since the orbit is circular). The subscript 1 in AV^ does not 
refer to the first of the two m-plane impulses, but to the sum of the 
two impulses If we denote by AV 2 the velocity necessary to correct 
the inclination alone, then AV 2 = V Ai The total required velocity, 
denoted by AV^, is the sum of the two 

av 3 = AV 4 + av 2 


The small thrust provided by each thruster, about 0 05 pound, 
prevents the velocity correction from being completed m one orbit It 
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PROBABILITY 


Table 3-1. Injection Covariance Matrix 



Velocity 

Path Angle 

Altitude 

Azimuth 

Latitude 

Longitude 

Velocity 

sec) 

0 3568 x *0 3 

-0 $994 x 10“* 

-0 5206 x 10 6 

0 1084 x 10° 

0 2643 x 40 1 

0 5708 x 4 0° 

Path angle 

(dog) 

Altitude 

(ft) 


0 2393 x i(T 3 

0 9349 x IQ 2 

0 6893 x 10~ 5 

0 3750 x 10’ 3 

0 8129 x 10~ 4 



0 7598 x 10 9 

-0 1594 x 10 3 

-0 3835 x 40 4 

-0 8281 x 1G 3 

Azimuth 

(deg) 




0 2348 x iO" 3 

0 9078 x 10" 3 

0 1887 x 10' 3 

Latitude 

(deg) 


Symmetnx matrix 



0 2251 x 10 _i 

0 4861 x 10' 2 

Longitude 

(deg) 






0 1050 x 10~ 2 



Figure 3-5 
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is necessary to perform the correction in bursts over several orbits, 
each burst contributing about 0. 58 ft/sec. In some cases 20 or more 
orbits, l. e. , up to two days, are required to complete the in-plane 
correction alone. If it is assumed that the in-plane correction is 
completed before the out-of-plane correction is made (although there is 
no particular reason for doing this), then the total maneuver may require 
four days or more. 

The covariance matrix for the Delta inertial guidance system used 
in establishing the orbit correction requirements is shown m Table 3-1 

Table 3-2. Injection Statistics of Selected Parameters 


Parameter 

Mean 

Standard Deviation 

Parameter 

Mean 

Standard Deviation 

Apogee radius 
error (n mi) 

l 88 

3 56 

Out-of-plane 
correction velocity 
(A/ 2 ft/sec) 

4 6 

3 5 

Perigee radius 
error (n mi) 

-l 78 

3 41 

Total correction 
velocity 
(£V 3 ft/sec) 

14 3 

6 4 

Inclination 

(deg) 

0 0 

0 013 

Uncorrectable 
inclination error 
(deg) 

0 0 

0 0 

In-plane 

correction velocity 
ft/sec) 

9 7 

5 4 





3(7 = high apogee error ~ 12 56 n mi 
3(7 - low perigee error = '12 01 n mi 


Although listed for completeness, the variance and covariances of the 
longitude error were not used in the analysis since longitude errors do 
not affect any of the critical parameters Table 3-2 shows how this co- 
variance matrix translates into other parameters of interest The cumu- 
lative distribution functions for the m-plane (AV.), out-of-plane (AV 9 ), 
and total (AV^) velocity requirements are shown in Figure 3-5. The 99. 8 
percent confidence -of-not- exceeding the in-plane and out- of- plane velocity 
requirements (3cr on a normal distribution) are 17 ft/sec and 31 ft/sec, 
respectively. The corresponding figure for the combined requirements 
is 39 ft/sec or about 9 ft/sec less than the sum of the two components 
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4. IMPROVING GDHS PERFORMANCE THROUGH 
OBSERVATORY MODIFICATIONS 

4* 1 INTRODUCTION 

The major portion of the ERTS study effort has been directed toward 
the design of a spacecraft system and a ground data handling system which 
meets the mission specifications Studies have also been undertaken to 
ensure that the interface between the two leads to the desired mission 
performance with a mmimum overall burden to the system This section 
summarizes certain studies related to the spacecraft and payload which 
could improve GDHS performance. 

The effort with the biggest impact has been attitude determination. 
The remaining parts of this section deal with that subject. Other tasks 
were also performed and are covered m detail m other volumes of this 
report. Those results are summarized here. 

An important source of information that can be obtained from the 
RBV ! s is the multispectral imagery For this imagery to have its greates 
value, a degree of color fidelity is required This ensures that identical 
phenomena always appear the same from image to image An evaluation 
of the basic fidelity of the system without calibration indicates that a 10 
to 15 percent photometric accuracy can be held over a year As 
calibration could improve this accuracy, the GDHS implementation and 
the performance improvement of havmg RBV calibration signal* have 
been investigated Although p ref light calibration data can be used, its 
accuracy will degrade over a period of time. To overcome this degrada- 
tion it has been suggested that the RBV's carry a light source with 
accurately known des Crete intensity levels. The resulting light source 
generated RBV image would be processed on the ground system to deter- 
mine the photometric response of the RBV camera in terms of picture 
element location and light intensity. The updated calibration curves would 
then be used in recording the imagery A discussion of the impIerrentatio> 
and use of RBV calibration data is contained m Section 5, Volume 17 

Going from a received video signal to a deliverable photograph 
requires several off-line functions essential but not directly associate^ 
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with the image reproduction. These include the generation of annotation 
tapes and the reduction of attitude data. Improved image processing 
efficiency results from carrying out these off-line functions ahead of 
receiving the video tapes. If the required information exists m the narrow 
band telemetry data which can be forwarded over data lines from Alaska 
and Corpus Christi, the prior preparation can take place before the arrival 
of the shipped video tapes. Key to all information processing is the time 
of RBV shuttering. From this time, spacecraft position, spacecraft 
attitude, sun angle, etc. , are all established. Thus provision has been 
made for the RBV shutter time to be inserted m the narrowband telemetry 
mam frame. 

Reseau points are used to determme the RBV geometric image 
distortions The reseau array format proposed for the ERTS sensors is 
shown m Figure 5. 15 . At the GDHS, the 9x9 reseau grid is subdivided 
into 16 squares with three reseau points to a side for a total of nine points 
each. A two-dimensional quadratic polynomial is then fitted to the 
distortions at the nine points of each subregion Each reseau point, 
however, occupies some space on the photograph and since the reseau 
point is completely opaque, all photographic information beneath the 
reseau points is lost and can only be partly restored by interpolation. 

The number of reseau points and their arrangement is also fixed. 

TRW has investigated an alternate reseau approach where the low- 
density, high contrast pattern normally used would be substituted with a 
high density low contrast pattern such as a finely gridded checkerboard 
of very low contrast. The thought behind this is that no picture information 
would be completely destroyed and, after geometric correction, the barely 
visible checkerboard would be completely removed by an inverse process. 

Reseau extraction would be accomplished using a correlation 
technique similar to that mechanized m the Itek PPR where the reference 
image would be an undistorted reseau pattern. Such a system could handle 
quite high order distortions depending on reseau contrast employed. 

It was visualized that a reseau pattern of this type would be projected 
onto the RBV faceplate through an auxiliary optical system. Using a 
variable intensity projection lamp the reseau could be removed at any time, 
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altered to match scene contrast, and used for photometric calibration 
Details of the tradeoffs between required contrast ratio, correlation 
region, and mesh dimensions is covered m Volume 17, Section 2, 

Three study areas comprise the attitude determination tasks filling 
the balance of this section For ERTS, the absolute pitch and roll attitude 
of the spacecraft are obtained from the earth horizon scanners. Volume 4, 
Section 6 discusses the use of these measurements and states the 
necessity to periodically update the radiance model to achieve the required 
two mile absolute location accuracy. The radiance model is rather slowly 
varying and hence calibration is entirely feasible Details of the horizon 
model calibration and updating are described in Section 4 2. 

The use of Kalman filtering to establish precise relative attitude 
data is also discussed m Volume 4, Section 6 6 Here the discussion 
of attitude determination is completed by going into the mechanization of 
the ground process mg required. The general approach is to smooth the 
observation data by equations of motion. Section 4. 3 discusses the 
dynamical system of equations governing the attitude of the spacecraft. 

The relationship between the observations and the state variables of the 
system equations are also given. 

The complexity of the dynamic modeling required to give accurate 
attitude data, coupled with the discovery of a mechanically and electrically 
interchangeable gyrocompass rate gyro of lower drift rate, stimulated 
further analysis of a multi gyro attitude, determination system A two- 
gyro system (one additional gyro) appeared most attractive. The second 
gyro whose input axis lies along the roll axis permits a direct measure- 
ment of the angular -momentum unbalance of tape recorders, thus easing 
a difficult modeling problem It also permits a more accurate measure- 
ment of absolute yaw through the ability to separate cross-coupled roll 
noise. The two- gyro attitude determination system is discussed m 
Section 4 4. The section gives a complete analysis of the two- gyro 
configuration - one gyro along the negative roll axis and the second one 
along the negative yaw axis The performance of the configuration is 
determined by the covariance equations of the Kalman filter. 
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4 2 EARTH HORIZON SURFACE MODEL 


4 2 1 Introduction 

In order to obtain accurate roll and pitch data from the earth horizon 
trackers, an accurate model of the Earth's horizon surface profile must 
be known. This section presents the formulation of a mathematical model 
which may be used to accurately define the horizon surface profile The 
formulation is based on the fact that from a known spacecraft position and 
a well defined horizon model only two orthogonal horizon measurements 
are required to establish body pitch and roll The measurements of the 
remaining two scanners provide redundant information and will com- 
pletely define the parameters of an ellipsoid which is assumed to charac- 
terize the horizon model. Because of horizon- scanner noise, correlation 
over many measurements is performed m establishing parameters for 
use m characterizing the horizon model This correlation approach per- 
mits an inclusion of important variations in the horizon profile due to the 
Earth's oblateness and the seasonal effect of the atmospheric radiance 

The horizon model is assumed to be an ellipsoid with azimuthal 
symmetry A number of ellipsoids are used to curve-fit the measured 
data, each assigned to a specific band of latitudes - Mathematically 
speaking, the size of the major and minor axes for each ellipsoid is 
determined as a function of the latitude band. 

4 2. 2 Analysis 

In this section the semi- major axis, a, and the semi-minor axis, b, 
of the earth horizon profile ellipsoid will be determined. The data avail- 
able for determining these parameters are the pitch plane central angle 
and roll plane central angle that the ellipsoid subtends from the spacecraft 
orbit position. This is illustrated m Figure 4-1 The central angles are 
obtained from the horizon scanner data which are telemetered to the 
ground station The problem of determining the ellipsoid can be conven- 
iently solved by expressing the ellipsoid in a coordinate system centered 
at the spacecraft. 
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Figure 4-1 

GEOMETRY OF THE CENTRAL ANGLES AND THE ELLIPSOID 


4. 2. 2 1 Equation of Ellipsoid 

In the analysis to follow, the definitions of three coordinate systems 
are needed Earth center inertial (ECI), Earth, center fixed (ECF), and 
orbital reference (OR), The definitions of these coordinate systems are 
given m Appendix A 

An ellipsoid to represent the earth horizon profile model can be 
expressed m ECF coordinates as 



2 

a 



or 


2 2 2 
X F + y F " a 



= 0 


(4-1) 
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Using the transformation from ECF to ECI, Equation (4-1) can be 
rewritten m ECI 


X 2 + Y 2 -a 2 + ^ Z 2 = 0 (4-2) 

b 

Where the identity of the forms of (4-1) and (4-2) is easily explained from 
the symmetry of the ellipsoid about the polar axis 

Next write the equation of the ellipsoid in OR coordinates through 
transformation given by Equation (A-4), Appendix A (Note ascending 
node fl is arbitary) Using the components of T 2 (i, j) with indicating 
column and j the row 

[t 2 (1, i)X R + T 2 (2, 1)Y r + T 2 (3, 1) (Z r -R)J 2 
+ ^T 2 (1,2)X r + T 2 (2,2)Y r + T 2 (3,2) (Z r -R)] 2 -a 2 

+ ^[t 2 (1,3)X r + t 2 (2,3) y r + T 2 (3,3> (Z R - R)] 2 = 0 ^ 


To solve the problem at hand we want to find the angle, J3, between an 
arbitrary ray from the spacecraft just tangent to the ellipsoid and the 
nadir line (see Figure 4-2) The angle a defines the orientation of the 
plane P m which the ray and the nadir axis lie When a = 0 or 180°, 
the ray is in the pitch plane and when n = ±90 °, it is m the roll plane 
We begin more generally by allowing the ray to intersect the ellipsoid 
at a distance p from the nadir axis 

This ray can be written m the OR system as 

x R = p sm p cos a 

y R = p sm p since (4-4) 

z R = p cos p 

Substituting these into (4-3) gives 

2 

(\ x p + fijR) 2 + (\ 2 p + ji 2 R) 2 - a 2 + ^(\ 3 p + [i 3 R) 2 = 0 (4-5) 


4-6 




Figure 4-2 

ORIENTATION OF RAY IN THE ORB I TAR REFERENCE SYSTEM 
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where 


\ i = T 2 (l, i) sin p cos a + T 2 (2, 1) sin p since + T 2 (3, cos P 

\ 2 = T 2 (l, 2 ^ sin P cos a + ^(2, 2) sm p sm ci + T 2 (3, 2) cos p 

\ 3 = T 2 (i, 3) sin p cos a + T 2 (2, 3) sm p sm a + T 2 (3, 3) cos p 

f x 1 = -T 2 (3, i) 

1^2 = 

ix 3 = -T 2 (3,3) 


Multiplying out (4-5) and collecting coefficients of powers of p gives' 
P f X l + X Z + ^2 X 3 j + 2 P R ( x i M ‘l + X Z^Z + ^2 


, d 2/ 2 . 2 . a** 21 2 . 

+ R Ifij + p. 2 + ~2 |x 3 I - a =0 


(4-7) 


This can be solved for p to give 


k(\,(i, + \ 2l i 2 + ^ x 3 M 3 j ±R yj{k 1^1 + X. 2 n 2 + b 2 X 3 tl 3) * + X 2 + b 2 + ^2 h £ " |( z) 

(h 2 * k Z + ^2 X 3 2 ) 


(4-8) 


From (4 -8) two real roots are obtained This may be physically 
interpreted as the distanced between the two points of intersection and the 
nadir axis if the line of sight is tangent to the horizon surface, the two 
points of intersection degenerate into a single point and mathematically 
speaking, the quadratic equation can have only one real root Therefore, 
the condition for line of sight to be tangent to the surface is represented 
by a zero discriminant of the quadratic equation 
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Hence, from Equation (4-8), we have 


(h^ ~ ^i) + fe 2 


2 2 
(\ 2 (J .3 - \ 3 H- 2 ) + (^1^3 " ^3^1) 


a 2 /. 2 .2 , a 2 2 

-z[h 


= 0 


(4-9) 


The above equatxonprovid.es the geometric relations necessary for the 
determination of the horizon surface profile The equation is completely 
general and can be used for any orientations of the orbital plane and the 
tracker's line of sight It is noted from Equation (4-6) that the equation 
of tangency is a function of i, v, a, (3 , a, b, and R 

The angles t- and v and the radial coordinate R specify the position 
of the satellites and are inputs from the spacecraft ephemeris program 

4 2 2 2 Method of Solution 

Using Equation (4-9) we can now solve for the angle (3 in terms 
of a, b, a and the orbital parameters 

The form of Equation (4-9) expressed m powers of tan p is 

A tan 2 p-2Btanp + C = 0 (4-10) 

where A, B, and C are functions of l , v, a, a, and b 
The solution of this equation is 


B 2 - AC ±B 

A (4-11) 

where the sign of one of the roots has been reversed to yield both values 
of p positive The two solutions now correspond to the opposite horizon 
scanners m the same plane 

The sum of the two roots corresponds to the sum of the readings 
(0^ and 0^) of the corresponding horizon scanners independent of the 
actual body orientation 


tan P 1 ? 

1 , L. 


/ 
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Thus 


tan (Gj + © 2 ) = tan ((3 j + fi?) = 


± 


2 - AC 

A + C 


(4-12) 


where positive sign is implied for (0 ^ + 0^) < 90°and negative sign is implied 
for (e 1 + e 2 ) > 90° 

Now considering the measurements m both the roll and pitch planes 
two equations like (4-12) can be generated. From these the values of a 
and b are found 

Since some variations of a and b with the spacecraft latitude 
position are expected, the values of a and b are determined as functions 
of the latitude position If the variations are small, a single ellipsoid 
with appropriately averaged a and b may be used If the variations are 
large, a number of judiciously selected ellipsoids with locally averaged a 
and b may be used Since values of a and b m the pitch plane differ 
from those m the roll plane, two sets of ellipsoids, one for pitch and the 
other for roll, should be used 

Once the values of a and b are determined with sufficient accuracy 
from the foregoing correlation equation, these values may be substituted 
back into Equation (4-11) to determine the two angles and (3^ Now the 
angle between the geocentric axis and the bisector can be determined by 


6 


\h ~ P 2 
2 


(4-13) 


Figure 4-3 is the geometric descrip- 
tion of measured and computed angle 
defining the directions of the geo- 
centric axis, the null axis, and the 
bisector In Appendix A the angle 6 
is used to specify the nominal orien- 
tation, NO, coordinate system 



Figure 4-3 

ATTITUDE ERROR 
determination 
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4 2 2 3 Simplifying Assumptions 

It has been shown m the foregoing analysis that the Earth horizon 
model can be determined by correlation of the mathematical model with 
the measured data from the Earth horizon trackers A number of 
simplifying assumptions have been introduced m the formulation of the 
model The first assumption was made when the horizon surface profile 
was assumed to be an ellipsoid This assumption facilitates the use of an 
analytical solution and is also considered reasonable since the mean sea 
level surface may be approximated by an oblate spheroid. The second 
assumption is that of the azimuthal symmetry This assumption is 
regarded as a reasonable one as evidenced by the Tiros 1 flight measure- 
ments of the Earth’s radiance The third assumption was introduced when 
the effects of the pitch and the roll attitude errors were uncoupled This 
was done when the pitch and the roll planes were assumed to be coplaner 
with the and the y^, z^ planes, respectively If the attitude errors 

m pitch and roll are sufficiently small, then the coupling effects are 
negligible In the case of large errors, an iterative scheme is suggested 
First, the attitude errors are determined by assuming an uncoupled effect 
Then, by transformations of coordinates using the computed values of 
pitch and roll attitude errors, the new attitude errors may be determined 
For each iteration, the angle a is determined from the aforementioned 
transformations 

i* 

4 2 3 Accuracy of Mathematical Model 
4 2 3 1 Introduction 

In Section 4 2 2, the formulation of a mathematical model for the 
earth's horizon surface profile has been presented The horizon model 
consists of a number of ellipsoids each of which represents the local 
horizon surface Presented m this section is an investigation of the 
accuracy of the proposed mathematical model For the basis of compari- 
son, a realistic horizon model has been constructed by using the horizon 
data presented m Figure (6 -6)of the final report, Volume 4 The mathe- 
matical model based on ellipsoidal surfaces is then compared with the 
foregoing realistic horizon model 
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4. 2 3 2 Analysis 


The latitude variations of horizon profile for the months of August 
and January were added to the mean sea level surface profile to obtain a 
realistic horizon model The August curve and the January curve were 
used respectively for Northern and Southern hemispheres by assuming 
a summer season in the Northern hemisphere and a winter season in the 
Southern hemisphere The Hayford International Ellipsoid, 
a = 6378 388 (KM), (b/a) = (296/297), has been assumed for the mean 
sea level surface The resulting horizon model is shown m Figure 4-4 
(solid curve) This curve is assumed to represent the three sigma 
Earth horizon surface 

The foregoing horizon model is approximated by a number of 
ellipsoidal surfaces The least squares fit technique is used to determine 
the mean values of the major and the minor axes. 



LATITUDE, B (DEG) 


Figure 4-4 

RADIAL COORDINATE OF EARTH RADIANCE PROFILE 
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First, using the radial coordinate R(<J> t ' of the horizon model, 
the x and z components are determined 

x i = R(<M cos 4> t (4-14) 

= R(<h) sin 4> t (4-15) 

where 4\ is latitude angle Then, the equation of ellipse may be written 
as 


Ax 2 + Bz 2 - 1 = 0 
1 1 


(4-16) 


where 

A = l/<a> 2 
B = l/<b> 2 

<\a> is a mean of the semi -major axis 

is a mean of the semi -minor axis 

According to the least squares method, the following quantity is 
minimized by setting its first derivatives with respect to A and B equal 
to zero 


Q = 



+ Bz 2 
1 



(4-17) 


where M is total number of data points to be used for the curve-fit 
then. 


9Q 

BA 



+ Bz 2 

l 


9Q 

3B 



+ Bz 


2 


i 


H 2 = ° 


') 2 2 = ° 


(4-18) 


(4-19) 
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Solving the above equations for <a> and<b]>, we get 


o> = 


k 4 )K 4 )-K z, ,f 

K 2 )( s ’ 4 ) - ( s * 2 * 2 )K 2 ) 


( 4 - 20 ) 


<b> = 


K 4 )K 4 ) - K 2 * 2 ) 


2 1 


K 4 )K 2 ) - K 2 1 2 )K 2 ) 


1/2 


(4-21) 


Thus, using the above equations, the ellipsoidal surfaces representing 
the horizon model have been determined 


4 2 3 3 Results 


A number of ellipsoids were used to fit the three -sigma earth 
horizon model One, two, three, and nine ellipsoids with the following 
ranges of latitude were selected 


a) One ellipsoid 

b) Two ellipsoids 

c) Three ellipsoids 

d) Nine ellipsoids 


0 = (-90 ~ +90) 

0 = (_90 ~ 0) 

0 = (0 ~ +90/ 

e = (-90 ~ -30) 
6 = (-30 ~ +30) 
0 = (+30 ~ +90) 

e - (-90 ~ -70) 
e = (-70 ~ -50) 
0 = (-50 ~ -30) 
0 = (-30 ~ -10) 
0 = (_10 ~ + 10 ) 
0 = (+10 ~ +30) 
0 = (+30 ~ +50) 
0 — (+50 ~ +70) 
0 = (+70 ~ +90) 


4-14 



No attempt was made to match the value nor the slope of the radial 
coordinate at the boundary points However, the maximum deviations 
are observed to be 0 014 percent m value and 13 percent m slope The 
values of radial coordinate based on these ellipsoids are shown m 
Figure 4-4 The correlation with the three -sigma horizon model is 
observed to improve with the number of ellipsoids Figure 4-5 shows 
the variation of RMS deviation m radial coordmate with the -selected 
number of ellipsoids It is observed that a single ellipsoid approximates 
the horizon model to ±1 18 (KM) 



TIME 

Figure 4-5 

RMS DEVIATION OF ELLIPSOIDAL MODEL 
from earth radiance profile 
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4 3 SYSTEM DYNAMIC EQUATIONS 


4 3* 1 Introduction 

In section 6 of Volume 4 a discussion of an attitude determination 
technique is presented which makes use of horizon scanner, rate gyro, 
and reaction wheel position data to establish time histories of both absolute 
and relative attitudes of the spacecraft for use m assuring accurate image 
reconstruction and annotation As an aid m carrying out ground processing 
of the spacecraft signals, evaluate total performance and suggest areas of 
improvement, a set dynamic system equations have been developed The 
various coordinate systems and transformations between them necessary 
for the development are given m Appendix A 

4 3 2 Dynamical Equations 

4 3 2 1 Euler Angle Equation 

The attitude of the spacecraft is expressed in terms of Euler angles 
defined between the true body (TB) coordinate system and the nominal 
orientation (NO) coordinate system The dynamical equations are now 
derived for both the spacecraft and the reaction wheels The equation of 
motion for the wheel is conveniently expressed m terms of the TB system 

Inertial Body Rate Equation 

The angular momentum H of the spacecraft in body coordinates is 
given as. 


H = Ico^ + I h 
— — B s — s 


•m* 

E n i & 
“111 


1=1 



where, in body coordinates, 


(4-22) 


I = Total spacecraft moment of inertia tensor including the 
instantaneous static inertia of the solar array, all wheels, 
tape recorders, etc 


w = Inertial angular velocity vector of spacecraft 
B 

= (u) x » ujy.» fx> z ) expressed in body coordinates. 

I = Solar panel moment of inertia expressed m the body 
coordinate system 
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il = Angular velocity vector of solar panels relative to 
spacecraft 

= Unit vector defining ith wheel axis with respect to the 
body axes 

I = Scalar moment of inertia of ith wheel about wheel spin 
axis 

o 

Q = Scalar angular velocity if the ith wheel relative to the 
spacecraft 



Angular momentum vector of the jth tape recorder relative 
to the body 


In the remainder of the system equation development let f2 g = 0 
since, during picture taking, the solar array is not driven It will be 
assumed that the tape recorders angular momentum are adequately com-* 
pensated for, and hence will be neglected If later information indicates 
this is not so, the procedure similar to deriving the reaction wheel equa- 
tions can be used for developing the recorder equations. We will therefore 
proceed to develop the system equations based on the angular momentum 
equation 


H - 




o 

£2 

1 l 


(4-23) 


where I is a known function of the solar array orientation with respect to 
the body 

Now the external torque _T applied to the spacecraft equals the rate 
of change of angular momentum, 

d H c 

T =-g F = (H) + w. B (4-?4) 

where, 

d ( ) • 

— = { ) = total inertial time derivative of ( ) 

o 

( ) = time derivative of ( ) m body 
coordinates. 
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Since co x co =0, it follows that 
Jo Jj 


co 

-B 


dt 


= to 


*f- CO 

-B 


x — B 


= CO—. 

-B 


and since ti is constant in body coordinates 



so that equation 4-24 can be written as 


T = 





(4-25) 


(4-26) 


(4-27) 


where all vectors are in body coordinates and C (w^) is the cross product 
operator m matrix form. 



0 

- w z 

CO 

y 

C(w B ) = u B x = 

CO 

z 

0 

~w 

X 


-CO 

(O 

0 


L y 

z 



Rearranging yields the body rate equation. 




1=1 


+ y n i h 

i i 

1 = 1 


(4-28) 


(4-29) 


Using Wg we can derive the Euler acceleration of the spacecraft with 
respect to the NO system. The inertial rate of spacecraft can be thought 
of as the sum of the two components, one due to rate of the body with 
respect to the NO system and the other due to the rotation rate of the NO 
system with respect to ECI 
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Let w o j = total rotation rate of the NO system relative to the ECI 
system expressed in ECI coordinates. It can be assumed that fl (ascend- 
ing node) and t (inclination) change slowly with respect to v (true anomaly) 
(see Section 3 of this Volume). It can also be tentatively assumed that the 
angular rate of the NO coordinate with respect to the OR coordinates is 
negligible 

Now express j m the OR system by equation (A-4), 

= T 2“o,I = (0> ° )T (4-30) 

Transformation T^ takes into NO coordinates, and takes it into 
body coordinates Thus the body rate relative to the NO system expressed 
in TB coordinates is, 

22 B ' S 1 T 32o (4 - 31 > 

In addition, applying transformation from equation (A-17) gives Euler 

rates. 

e = S 2 [ Sb - Sl T 3£o ] (4-32) 

where the Euler angle rate vector is defined as, 

e= (4>, 0>) T (4-33) 

where nominally, 

<i> - roll 

6= pitch (4-34) 

4* = yaw 

Taking the derivative of (4-32) yields, 

£ = S 2 — B + S 2 -B " S 3 T 3^o (4 ~ 35) 

where, 

S 3 = S 2 S 1 (4-36) 
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Note that equation (4-35) can be integrated twice using equation 
(4-29) directly, or equations (4-29) and (4-32) can be simultaneously 
integrated to give Euler angles. The latter method avoids some , 
matrix differentiations given in Appendix A. The choice of integration 
method depends on a study of numerial integration accuracy and speed. 

From (4-32), note that, 


-B ~ S 2 - + S l T 3 a o 


(4-37) 


Equation (4-35) is the desired system equation for the body To 
integrate (4-35) requires knowing the external torques and the reaction 
wheel equations. The external torques are considered first 

4 3. ?. 2 External Torques 

Some tentative estimates of ERTS external torque magnitudes are 
as follows 

Source Max. Value (lb-ft) 


Residual magnetic moment 

2. 6 x 10 ^ 

Gravity gradient 

1.4 x 10" 5 
/ 

Solar pressure 

2. 0 x 10" b 

Aerodynamic pressure 

- 


The last source is considered to be negligible, the first two are 
easily modeled and comprise over 95 percent of the external torque. Solar 
pressure effects appear too small to warrant the complicated model 
required to describe them. Thus the external torque T is modeled as 


— " — GG + - »M 

(4-38) 

1 RM - 22* 

(4-39) 

r- 3(1 A t A 

— GG - T3 - x l - 

(4-40) 


where, m TB coordinates, 

m = residual magnetic moment vector (ft lb/gauss) 
B>k = Earth magnetic field flux density model (gauss) 


4-20 



[x = Earth gravitational constant (1,4082 x 10*^ £t3/sec^) 
R = distance from the Earth's center to the spacecraft (ft) 
r_ = unit vector along the OR z axis (ft) 

2 

I = spacecraft moment of inertia matrix (lb ft sec ) 


Other torques can be considered as 
state noise. 

Let B^j the magnetic field 
flux density of the spacecraft posi- 
tion, be expressed m the radial, 
tangential unit vector directions 
e r and e^ respectively, as shown 
m Figure 4-6. It is given by 





Figure 4-6 

MAGNETIC COORDINATE 
unit vectors 

— b = -B 0 [(2 e m ) ^ - (cos eje t ] (4-41) 

where 


sm 0 


_ 

m / 2~T 2\ 1/2 

\ m + z m 


(4-42) 


cos 0 


■(- 


. - sm 2 0^ | 
m \ ml 


1/2 


(4-43) 


sm <f> = 

m 


m 


( x m 2 + y m 2 ) 


1/2 


(4-44) 


COS kv, = 
m 


m 


( x m 2 + ym') 


1/2 


(4-45) 
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(4-46) 


e 

-r 


cos 6 m cos 4> m 
cos 6 m sin 4> m 
sin 6 

m 


% 


- I sin e m I cos * m 

- | sm 0 m | sm <j> m 

cos 0 

m 


(4-47) 


Based on the transformations defined m Appendix A the magnetic 
m the TB system is' given by the sequence of matrix operations 

-b = S 1 T 3 T 2 T 1 T 4- (4-48) 


4. 3.2.3 Wheel Rate Equation 

For one wheel, define a cartesian coordinate system fixed relative 
to the spacecraft with one axis along the spin axis of the wheel. The 
angular momentum vector of the wheel H w is"then given as 




where 

A - wheel moment of inertia tensor 

CO = inertial angular wheel rate m wheel coordinates 
“W 

The total torque on the wheel T is given as 
^ w 



(4-49) 


(4-50) 
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where, m wheel coordinates 


T 

— M 
T 

— B 


K = 


CO 


B 

Dw 


motor torque 

bearing torque 

dynamic friction coefficient 

body inertial rotation rate 

unit vector in direction of wheel axis 


since, 


d H w c 

- w ~ dt “ -w + — B X -w 


( 4 - 51 ) 


T + T - K 
-M -B 


u 


‘w 


( rj lri = (Au ) + to 
\— B -w/dwj \ — w/ — 


„ x Aco 
B — w 


( 4 - 52 ) 


Take the dot product of this equation with the wheel axis unit vector 


T„ - K w - t) = I co 
M 1 w — B -w w 


1 = I co - n "(co x Aco | ( 4 - 53 ) 

r J w -to \~B ~ w/ 


Note 


to = <o i) 
w w— w 


( 4 - 54 ) 


and that the dot product 


n T- = 0 
Aw — B 


( 4 - 55 ) 


based on the assumption, quite valid, that the bearing torque is perpen- 
dicular to the spin axis 

Also by symmetry of the wheel, rj is a principal axis of A, i e , 

An = In ( 4 - 56 ) 

Aw A w v ' 


where I = moment of inertia of the wheel with respect to the spin axis 


4-23 



By a vector identity and (4-54) and 4-55) we will show that last term of 
(4-53) is zero. By vector identity, we have 

ri T (to x Aw ) = w t, T (Aw x -n ) 

-w \ B — w / — B \ — w -Lw / 

By utilizing (4-54) and (4-55) we obtain, 

(Au XT] \ = Iw w T/ n x ri ) = 0 (4-57) 

— B \ — w Jw/ w B _Lw/ 


which was to be shown. Thus (4-53) becomes 


M 


KL - (w T r, ) 

[_ w \— B — w/ 


= Id) 

— w 


(4-58) 


Using (4-58) we wish now to obtain equation for S2, the relative angular 
acceleration of the body with respect to the body. The reason for this is 
that the wheel data is relative to body. 


Now recall. 


n = w 


Differentiation gives. 


Thus, 


0 

ft = Angular rate of wheel relative to body 

- (iw T -b) 


°o X 

to = S7 - n to 

w - — B 


o /oo T. \ 

T M- Kns l( n -3w T £- B ) 

Now add subscript i to denote the ith wheel and note. 

T T 

Hw -B = ^B 


(4-59) 


(4-60) 


(4-61) 


(4-62) 


where n and are in body coordinates. This gives the ith wheel equation. 


K T 

o° 1 A Ml T 

Q 1 = ”^1 + T“ ~ — i — B 


(4-63) 
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With the wheel equation now derived we can summarize the complete 
dynamical system of equation 

From (4-35) 

e = S 2 co B + S 2 o> B - S 3 T 3 co o (4-64) 

and 


where 


°p 

n = - 





(4-65) 


Terms 

Given by Equation 

S 2 

A-17 

S 2 

A-22 

S 3 

A-23 

u) 

“B 

(4-29) 

U> ^ 

— B 

(4-37) 

OJ 

(4-30) 

— o 


4 3 3 Measurement Equations 
4 3 3 1 Intr oduction 

From various data telemetered down from the spacecraft very 
precise estimates are to be made of Euler angles and rates These 
Euler angles combined with the orbit coordinate axis directions and 
ephemeris data, provide sufficient information to determine the center 
and corner coordinates of the image frame 

Telemetered attitude data which needs to be related to the Euler 
angles and Euler rates are 

• Horizon sensor angles A, B, C, D 

• Rate gyro readings j 

For angles A, B, C, D see Figure A-l 




The reaction wheel position data are also available. The relationship 
between it and the corresponding wheel 1 of (4-60) is obvious. 

4. 3 3.2 Horizon Scanner Measurement 

The angles A and C are related to roll, and B and D are related 
to pitch — the precise relationships to be derived — the horizon sensor 


angle differences are shown in Figur 

are all zero, the body axis 

aligns with the orbit coordinate 

Z axis For nonzero Euler angles, 
o 

the difference m horizon head angle 

measurements A0^, AG^ can be 

related to any shift of from the 

Z direction 

o 

A0 1 = A ~ C (4-66) 

Ae 2 - V 2 <4 - 6?) 


4-7 If the spacecraft Euler angles 



Figure 4-7 
Z BODY AXIS 

and sensor angle differences 


The effect of yaw does not change the A, B, C, D readings by any 
significant degree. Thus, A ^ A ^ can be related directly to roll (<{>) 
and pitch (6) as is clear from Fig- 
ures 4-7 and 4-8. From Figure 4-8 
note. 


Z, - BODY Z AXIS UNIT VECTOR 

<10 = 8 D 
2 Z 


|4>J 

sm (A0 2 ) 

f-b^xyl 

' cos ( Ae i) 


(4-68) 

A8 1 

(4-69) 


Thus the horizon measure- 
ments m terms of Euler angles are, 

A - C 

2 



Z BODY AXIS 
and euler angles 

A0j = 4> (4-70) 
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4. 3 3 3 Rate Gyro Measurements 


In Figure 4-9 two rate gyro measure- 
ments are shown m the x, z plane If 
the rate gyros have a direction unit vectors 
and 2 ]_q z respectively then. 
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(4-72) Figure 4-9 

RATE GYRO 
me a s ur ement s 


where ( w x > w y> w z ) are the inertial body rate components 


= <"*• V V 


(4-73) 


Recall from Equation (4-37), 


— B = S 2* L i + S l T 3 ~o 


(4-74) 


Substituting above. 


where 
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4 4 TWO -GYRO PERFORMANCE ANALYSIS 


4 4 1 Introduction 

This section presents an analysis of an attitude determination 
system with two horizon scanners and two rate gyros, henceforth referred 
to as the two-gyro system, with respect to mapping accuracy. , First, 
a quantitative summary of the effect of the two- gyro system attitude 
errors on the mapping accuracy of the RBV and MSS is given Next, a 
detailed explanation of the source of the errors is presented A simplified 
error model is derived from the exact equations of motion for the space- 
craft The effect of all assumptions made m obtaining the error model 
are stated A covariance analysis of the Kalman filter mechanization 
employing the model is performed Finally, a summary is given tabulat- 
ing the effect of both the assumptions used m deriving the error model and 
the errors from the Kalman filter convariance analysis 

4 4 2 Development of the Two-Gyro Error Model 

The following paragraphs derive an attitude error model for the two- 
gyro system mechanization The model is derived from the basic equa- 
tions of motion and the Euler angle rate equations. Expressions for the 
Kalman filter measurements are also presented The effect on attitude 
accuracy of all assumptions and approximations is given 

4 4 2 1 Equations of Motion > 

The equations of motion derived m Section 4. 3 are repeated below in 
slightly modified form 

= I B I [ T - Vfe - VA - *T - C ( I b“b + Ws + VA - H t)] 

(4-77) 

« =I -1 [t - K £2 1 - Aq 1 u-o (4-78) 

ii[miiJ0B ' ' 
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where 


Wg = inertial angular rate of S/C m S/C reference 
coordinates — 3x1 vector 

^ = angular rate of reaction wheels about their spin 

1 axes relative to the S/C — 3x1 vector 

= angular rate of solar panels about their spin axis 
relative to the S/C — 3 x 1 vector 

T = external S/C disturbance torques m S/C coordi- 
nates — 3x1 vector 

T = reaction wheel motor torques — 3x1 vector 

m 

H,p = tape recorders total momentum m S/C coordi- 
nates — 3x1 vector 

I_ = inertia tensor of S/C m S/C coordinates — 3x3 

B matrix 

I = diagonal matrix giving inertias of reaction wheels 
about their spin axes — 3x3 matrix 

Ig = inertia tensor of solar panels m S/C coordinates — 
3x1 vector 

K = diagonal matrix giving damping coefficient of reac- 
tion wheels about their spin axes — 3x3 matrix 

A fi = non- orthogonal reaction wheel spin axes mi s align - 
J ment matrix — 3x3 matrix 

and 
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The reaction wheel misalignment matrix is defined as 


1 

0 

0 
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where 

G = misalignment of 1 — reaction wheel spin axis into the j 
1 - 1 axis of the S/C coordinate system 

It has been decided from an operational viewpoint to stop the motion 
of the solar panels relative to the spacecraft body during the time pictures 
are taken Therefore, 


o 



From the relation 
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one can obtain the torque components from (4-77) 
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The equation for involves approximately the same relative mag- 
nitudes as that for so it will not be presented here Whatever simplified 
form is derived for will be used for The following paragraphs will 

examine the magnitude of various terms m (4-79) and (4-80) m order to 
arrive at a simplified though accurate approximation of (4-77) 

The known parameters are 

I = 210 slug ft 2 
X B 

I = 250 slug ft 2 

I = 390 slug ft 2 

Z B 

I = 3 05 x 10" 3 slug ft 2 
w 

I =3 05 xl0 3 slug ft 2 
^w 

I = 12 20 x 10“ 3 slug ft 2 
z ° 

w 
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Approximate RMS values for the other parameters are 


v-5 


o> = co =6x10 rad/sec 
x z 


_ 3 

co =1x10 rad/sec 

y 


£2^ = £2^ = £2^ = 65 rad/sec {1/2 saturation speed) 


T = T = T =1x10 ft lb 
x y z 


0 s4 x 10 rad (4 arc mm 3cr) 
ij 


I = 4 slug ft 
ij 6 


The tape recorders are aligned so that their total momentum vector 

_3 

lies normal to the S/C y axis to within 1 x 10 rad This results m 


H = H s 0 007 ft lb sec 

Z (ji 

H at 7 x 10" 6 ft lb sec 
y T 


The term consists of two separate torques as indicated m equation 
(4-78) The first torque is due to the reaction wheel motors The motors 
are turned on and off by the ACS and produce the following torques 

T = T = 21 x 10' 3 ft lb 

x y 

m J m 

T = 52 x 10" 3 ft lb 
z 

m 


4-32 



The second torque is due to motor windage At 1/2 saturation speed the 
torque values are 


K ft = K ft = 0 54 x 10” 3 ft lb 

xx y y 


K ft = 2 82 x 10" 3 ft lb 
z z 


Using the above parameter values, equation (4-79) can be rearranged as 


T, = /-I ft - H \ - u I ft 
1 I x x x„ y z z 
\ w T/ 7 w 

+ |T - co co /l -I \ + co I 

L x 2 n z b y B j Z ' 

+ [other torques] 


ft - to H 


y w y y z t 
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(4-81) 


where the right-hand side of the equation has been grouped into terms with 
expected values of 


1st term 

<1x10 2 ft 

lb 

2nd term 

^0 8 x 10~ 3 

ft 


-5 


3rd term 

=1x10 ft 

lb 

4th term 

<1 x 10" 5 ft 

lb 


Similarly, equation (4-80) can be rearranged as 


T, = /-I ft 

M 7 


H 


CO I ft - CO I ft 1 
Z X X X z z 

w w / 


+ Ty + (other terms) 


(4-82) 
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where the right-hand side of the equation has been grouped into terms 
with expected values of 


1st term 
2nd term 
3rd term 
4th term 


<1 x 10 -2 ft lb 
=5 x 10" 5 ft lb 
m x 10" 5 ft lb 
<1 x 10' 5 ft lb 


As noted above, involves approximately the same magnitudes as so 
it will not be presented here 

o 

It should be pointed out that the short term effect of the torques 

is to produce the spacecraft limit cycle From the ACS analysis it is 
known that the limit cycle worst case magnitude is 0 18° peak Since the 
motion is very random, assume an RMS of 0 10° The torque terms 
involving the misalignment of the reaction wheels will couple the limit 
cycle motion of another axis into the axis being considered Since the mis- 
alignment angles are about four arc minutes (3<r), the attitude error due to 

_ 5 

misaligned wheels is about 4x10 degrees Therefore, this term was 
included m the category of "other torques 11 

The same limit cycle consideration is applicable when the above 
torque expressions are substituted into equation (4-77) The substitution 
gives 

^ = C 1t i + (l _i T, + I _1 tJ 

U) = I _1 T, + fr 1 ! 1 , + rV ) (4-83) 

y y B 2 \ yx 1 yz 3/ 

w = I" 1 T* + (l~ 1 T 1 + I _1 T,) 
z Zg 3 \ zx 1 zy 2/ 

The products of inertia are m the neighborhood of 2 percent of the moments 
of inertia for the spacecraft If their effect is neglected, they will intro- 
duce an attitude error due to gross axis limit cycles of 0 002° This 


4-34 



error is considered negligible for spacecraft attitude determination m 
view of the larger magnitude of other uncorrectable errors m the overall 
system With this approximation, the above equations reduce to: 
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The two- gyro system mechanization uses a rate gyro along both the 
x and z spacecraft axes to measure relative attitude motion No informa- 
tion concerning the reaction wheel torques m these axes is used There- 
fore, the approximation is made that 
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For purposes of a covariance analysis, the assumption is made that Tj 
and Tg are correlated random error torques with a standard deviation of 
1 x 10” ^ ft lb Their correlation time is short enough to preclude their 
being estimated to any effective extent. That is 
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where 


E [ T;x 2 ] = E [ Tz 2 ] = x 10 2 ft lb) 


The attitude determination system uses reaction wheel revolution 
data to measure the relative pitch attitude motion of the spacecraft, 1 . e. , 
it implicitly determines the angular acceleration £2^. and solves for S> 
using the equation 


o J w o 

y I y 


A 2 percent error m the knowledge of the ratio of inertias will result m a 
scale factor error of 2 percent m estimating the limit cycle of the space- 
craft m pitch attitude This produces the same magnitude of error as 
was introduced by the products of inertia, i e , 0 002° 

In addition to using pitch reaction wheel data to determine pitch rela- 
tive attitude, the mechanization also uses roll and yaw wheel data to esti- 
mate the roll and yaw wheel momentum Since roll and yaw body rates 
are measured directly by the rate gyros, it is assumed that the second 
term of the right-hand side of equation (4-82) can be neglected There- 
fore, the approximation is made that 
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“ -I £2 
2 y y 



+ T 
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For purposes of a covariance analysis, the assumption is made that Ty is 
a correlated random error torque with a standard deviation of 1 x 10“^ ft lb 
The correlation time is short enough to prevent the torque from being esti- 
mated That is 


o o 


-I £2 - H + T 



(4-85) 
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where 


E T 2 = (1 x 10 5 ft lb) 


In summary, equation (4-77) has been reduced to the following 
equations 


where 




(4-86) 




= (1x10 
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ft lb) 
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= (1x10 
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rhe assumptions made were 


1) The solar panels were not moved during attitude estimation. 

2) H , H < 10~ 2 ft lb. 

X T Z T 

3) Gyroscopic torques m the y axis due to x and z axis wheel 
momentum can be estimated better than 1 x 10"^ ft lb. 


The approximations made were 
1) 


Misalignment of the reaction wheels was not corrected. 
This resulted m an attitude error of 4 x 10"^ degrees in 
all three axes. 


2) Products of inertia of the spacecraft were neglected. This 
resulted m an attitude error of 0. 002° m all three axes. 


3) The ratio of wheel to body inertia m the pitch axis was 
known to within 2 percent. This resulted in an attitude 
error m pitch of 0. 002° 


1.4. 2. 2 Euler Angle Rate Equations 

The Euler angle rate equation derived m Section 4. 3 is repeated 
>elow. 


OJ 


B 



+ S. GJ 
1 o 


(4-87) 


vhere 

e = Euler angles (yaw, pitch, roll, m order) describing the 
attitude of S/C reference axes relative to rotating orbital 
axes. 

w = inertial angular rate of orbital axes m orbital axes 
coordinates. 

Sj = direction cosine matrix describing attitude of S/C refer- 
ence axes relative to rotating orbital axes. 

.nd 


coso 

-smo 

cosO 

cos0sin<j> 

-smO 

cos0cos<t> 
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The orbital coordinate system is defined as follows 1) the origin is at 
the center of the S/C reference axes, 2) the z axis points to the center of 
the Earth, 3) the y axis is normal to the orbit plane and in the approxi- 
mate direction of the orbit rate vector, and 4) x xjy = z. 

The inertial rate of the orbital coordinates consists of two compo- 
nents 1) the orbit rate of approximately 210°/hr due to rotation about 
the Earth, and the orbit plane precession rate of approximately 0. 04° /hr 
due to the sun synchronous orbit. Although an attitude determination 
system can account for its presence, the assumption is made here that 
it can be omitted from the performance analysis without affecting the 
predicted results. Therefore, in orbit coordinates 



The largest Euler angle expected is 0.4° and all Euler angle rates 
are expected to be in the same order of magnitude. Therefore, the fol- 
lowing small angle approximations are made for all three Euler angles. 

cos e si 

i 

sm e = d 1=1, 2, 3 e = 4 >, e =0, e_ = 4 ’ 

ii '’12 3 

e » e e 

i i i 

Using this approximation, Equation (4-87) becomes 

co ~ <S> + ^ co 
x o 

co =0+co (4-88) 

y o v 9 

• 

O) s ^ - 4> OJ 
z o 

In summary, Equation (4-87) was reduced to Equation (4-88) using 
the following assumptions 
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1) Small angle approximations are valid. 

2) Omission of orbit precession rate does not affect the per- 
formance analysis accuracy. 

4. 4 2. 3 Measurement Equations 

Now that the simplified dynamical equations have been derived, 
measurement equations, are needed before the covariance analysis can 
be performed. Referring to Fig- 
ure 4-10, it is observed that a per- 
fect Earth sensor m a perfect 
attitude control loop will not point 
the spacecraft z axis at the center 
of the Earth. Instead, its attitude 
will lie somewhere between geo- 
centric and geodetic if the equipo- 
tential radiance surface tracked by 
Earth sensors is an ellipsoid fairly 
close to the Earth For reference, 
the maximum difference between 
geocentric and geodetic latitude on 
the Earth (at 45° latitude) is about 
0 2 °. 

In view of the fact that the maximum 6 shown m Figure 5. 1 will be 
less than 0 2° and that small angle approximations have been assumed 
valid for the Euler angles, then the Earth sensor pitch and roll signals 
are assumed to be 



Figure 4-10 

GEOMETRIC DEFINITION of earth 
sensor measurement m polar plane 


where 


V *n 
6 6’ 6 <l> 


e ESA = 6 - 6 e +e B + e n 


*ESA = * • 6 + + + 


Earth sensor bias errors, excluding misalignment. 

Earth sensor high frequency noise (f Q > 1 cps). 

error m prediction of the pitch and roll components 
of 6. 
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Misalignment errors are not included for either Earth sensor since it is 
assumed an initial calibration of attitude sensor plus camera misalign- 
ments will be made using pictures with ground truth points along with 
attitude estimate data. 

The error m predicting 6 is 0. 022° (l<r). The ESA bias is 0. 028° 
(Itr) and the noise is 0. 04° (lcr). Assuming that 6 is subtracted from the 
measured signal, the ESA signals are redefined as 


6 = 0 
ESA 


+ 0+0 

Jd 


n 


where 


^ESA “ ^ + 




+ <f> 

T n 


r„ 2 i 

r, 21 

e_ 

- E 

l B J 

L y b j 

r 2 ] 

r 21 


= e k 

L n . 

L 11 J 


= (0. 036 °) Z 
= (0 04°) 2 


The roll and yaw gyros measure roll and yaw inertial rates plus 
errors For the same reason as stated above, misalignment angles are 
not considered. The gyros have a bias error of 0 06° /hr (lo-) and a high 
frequency noise of 0. l°/hr (lcr). Their outputs are 


to 


4 > 


= CO + € 
X 


♦ + n $ 


where 


e[« 4 , 2 ] = e[^ 2 ] = (°. 06»/hr) 2 
E jn ^ j = = (0. l”/hr) 2 


At every 45° of the pitch reaction wheel revolution, a pulse is pro- 
vided and recorded and hence a quantization error of 45° is introduced. 
It is easily shown that error has a uniform probability density function 
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with zero mean and a standard deviation of 13° . Therefore, the reaction 
wheel signal, r^^., is 



w <j> = “x + e 4> + n c|> 

w v|, = + + 

\w = 11 e + n RW 

where 

e [ 6 b 2 ] = *[* B Z ] = (o * 0360)2 

E [ e n 2 ] = E [* n 2 ]=(°- 04 ") 2 

E [' 4 > 2 ] * E [ V ] = (0 06 ” /hr > 2 

E [ n 4> 2 ] = E [ n /] = (0 - I ° /hl ' )2 

e [ w ] = < i3 “> 2 

The assumption used in deriving the above equations is that the small 
angle approximations apply to both the Euler angles and 6. 
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4. 4. 2. 4 Error Model 

Equation (4-86), (4-88), and (4-89) are expressed m block diagram 
form m Figure 4-11. This error model will be used m the covariance 
analysis 

4. 4. 3 Covariance Analysis of Error Model 

The previous paragraphs developed a simplified error model 
and a set of measurement equations which made noisy observations of 
various states m the model. This section will determine the variance of 
the error in trying to estimate the three attitude angles m the model by 
processing the measurement data. In addition to estimating absolute 



Figure 4-11 

ERROR BLOCK DIAGRAM! FOR COVARIANCE ANALYSIS 
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pitch, roll, and yaw, an estimate of the relative change m pitch and roll 
during a 28 second period will also be made. 

In determining the estimation error, a Kalman filter mechanization 
was employed The Kalman filter was used to indicate the best results 
achievable without the necessity of a Monte Carlo analysis. 

4 4.3.1 Relative Pitch Attitude 

An error model for relative pitch attitude estimation is shown m Fig- 
ure 4-12. This model is the same as Figure 4-11 except for omission of 
the pitch bias error The bias will not affect the relative attitude estimation 
as long as the assumption that it is a constant during the estimation period 
is valid. Figure 4-12 also shows the tape recorder momentum as a step 
rate input instead of a torque. This simulates the effect of turning on or 
off the tape recorders. 

A Kalman filter mechanization using the error model m Figure 4-12 
was programmed on a digital computer and the covariance of the estima- 
tion error was recorded as a function of time for various error magnitudes 
and initial conditions. 

The first effect noted was that the estimation of relative pitch atti- 
tude was mostly unaffected by reaction wheel motor torques. The reason 
for this is that the motor torques produce a wheel motion identical to the 
body motion except for a scale factor equal to the ratio of wheel and body 



Figure 4-12 

PITCH ATTITUDE ERROR MODEL 
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inertias. The uncertainty m wheel position because of the quantized data 
is insignificant compared to the error due to other spacecraft disturbance 
torques. In starting the filter, initial spacecraft and wheel rates of 
0. 003 '/sec and 58 rad/sec were assumed. With a one second sampling 
rate, the filter converged on a steady state wheel rate estimate within 
20 seconds. Converging on a steady state body rate estimate produced 
the predominate time constant m the estimation process. The steady 
state values for body rate and attitude are a function of the magnitude of 
disturbance torque, T^ , assumed. It was found that variations m the 
disturbance torque time constant, t, produced little effect on attitude 
estimation for values below 100 seconds. In view of the many factors 
producing the disturbance torque, T^ , the assumption of a longer time 
constant was questionable. Figure 4-13 shows the effect of the magnitude 
of T on the standard deviation of the error m estimating 0. 

The maximum tape recorder momentum is 0. 007 ft lb sec. The 
tape recorders are aligned so that their total momentum vector is normal 
to the pitch axis within 1 milliradian. Therefore, turning the tape 
recorders on or off will result m a momentum change in the pitch axis 



Figure 4-13 

PITCH RELATIVE ACCURACY VERSUS DISTURBING TORQUES 
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of 7 X 10 6 ft lb sec. In turn, this will produce a body rate of 1.6 x 10 - ^ 
deg/sec. It was determined that this step change m rate did not signifi- 
cantly affect the estimation of pitch attitude. 

From the derivation of the error model it was found that a reason- 

- 5 

able value for the magnitude of T was 1x10 ft lb. Referring to 

Figure 4-13, this will produce a relative pitch attitude accuracy of 
0. 006° m 150 seconds from the time that the Kalman filter starts process- 
ing data. 

4. 4. 3. 2 Absolute Pitch Attitude 

The error model for absolute pitch attitude estimation is the com- 
plete pitch channel shown m Figure 4-11. As pointed out previously, the 
only difference between the absolute and relative error models is the 

bias 6_. The bias arises from radiance model and Earth sensor errors 
B 

and has a magnitude of 0. 036° (l<r). 

A digital computer analysis was performed using the absolute pitch 
error model m a Kalman filter mechanization. Within seconds the pitch 
attitude error leveled out at approximately 0. 036° This result was pre- 
dictable since the bias is uncorrelated with the relative attitude error and 
their net effect on absolute pitch attitude is the RSS of the two errors. 

4. 4. 3. 3 Relative Roll Attitude 

The error model for relative roll attitude estimation is exactly as 
shown m Figure 4-11 except that the roll Earth sensor bias is omitted. 

The bias was omitted for the same reasons explained m the paragraphs 
on relative pitch estimation. 

A Kalman filter mechanization using this error model was pro- 
grammed on a digital computer and the covariance of the estimation error 
was recorded as a function of time for various error magnitudes and 
initial conditions. 

Figure 4-14 shows the transient response of the error m the esti- 
mation of roll attitude. As would be expected, the rate gyro on the roll 
axis immediately estimates the initial roll rate of the spacecraft down to 
the neighborhood of the gyro noise. However, the filter must still esti- 
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0 SEC 100 SEC 200 SEC 300 SEC 

TIME 

Figure 4-14 

ROLL RELATIVE ACCURACY 

mate the rate introduced m roll by an initial yaw misalignment. This is 
the reason the error remains relatively constant from 50 seconds to 
100 seconds. After 100 seconds the error gradually decreases mostly 
due to the fact that the yaw estimate is also improving. 

Referring to Figure 4-14, after a time period of 300 seconds, the 
error has a value of 0. 005° (Icr). Three hundred seconds is chosen 
because the Kalman filter requires this long to reduce the error m the 
yaw estimate to an acceptable level. 

4. 4. 3. 4 Absolute Roll Attitude 

This error model for absolute roll attitude estimation is the com- 
plete roll and yaw channels shown m Figure 4-11. As pointed out pre- 
viously, the only difference between the absolute and relative error 
models is the bias <j>g. The bias arises from radiance model and Earth 
sensor errors and has a magnitude of 0. 036° (lcr). 

A digital computer analysis was performed using the absolute roll 
error model m a Kalman filter mechanization The result was the same 
as found for the absolute pitch case, 1 . e. , the resultant error was the 
RSS of the bias and the relative roll error, approximate! %' 0 036° (lcr). 
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4.4 3. 5 Absolute Yaw Attitude 


The error model for absolute yaw attitude determination consists 
of the complete yaw and roll channels shown m Figure 4-11. 

An examination of the error model indicates that yaw is estimated 
not only by processing the gyro signals, and but also by observing 
the roll Earth sensor signal. Use of the roll signal is not readily apparent 
m the ACS mechanization since the signal is not used m the yaw attitude 
control electronics. The premise for the yaw control loop, however, is 
that roll will be maintained at approximately zero. For example, if the 
body coordinates and the orbit coordinates are both rotating about their 
y-axis at orbit rate, but are initially misaligned by a yaw angle, then 
as time progresses the roll angle will increase. To counter this increase, 
the roll control channel will command a counter rate about the body roll 
axis. This counter rate will be equal to the yaw misalignment times orbit 
rate. The ACS gyro, which has a component of its input axis along the 
body roll axis, will measure this rate and correct the yaw misalignment. 

The above example also illustrates the basic time constant in the 
transient response of either the ACS or ADS yaw channel. Correction or 
estimation of the yaw misalignment depended on both detecting the roll- 
misalignment, which is proportional to m the presence of roll ESA 
noise and detecting the commanded roll rate, also proportional to u> Q , m 
the presence of gyro noise. Therefore, with an optimum estimation 
scheme, the transient response becomes a function of both orbit rate and 
sensor noise levels. 

> Figure 4-15 illustrates the result of a covariance analysis of the 
yaw and roll channels. The errors are broken up into two components m 
order to illustrate the transient response of the estimation error without 
the roll gyro bias included. At t = 300 seconds the total yaw error is 
approximately the RSS of the transient response and the error due to the 
gyro drift rate, namely 0. 038° (lor). 

4. 5 SUMMARY OF TOTAL ADS ERRORS 

Section 4. 4 2 derived an error model from the basic equations 
describing spacecraft attitude by making assumptions and approximations 
concerning 
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Figure 4-15 

YAW ABSOLUTE ACCURACY _ 

1) Solar panel motion 

2) Tape recorder torques 

3) Correction of gyroscopic torques from, reaction wheel 
momentum 

4) Misalignment of reaction wheels 

5) Spacecraft products of inertia 

6) Accuracy of spacecraft moments of inertia 

7) Small angle approximations 

8) Omission of orbit precession m analysis 

Section 4. 4. 3 used the error model m a Kalman filter mechanization and 
presented the covariance of the error in estimating the three spacecraft 
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attitude angles, both absolute and relative. Additional assumptions were 
required concerning* 

1) The time invariance of the gyro and Earth sensor bias. 

2) The autocorrelation function of the disturbance torques 
modeled. 

3) The statistical independence of one second samples of the 
gyro and Earth sensor high frequency noise terms. 

Table 4-1 presents a summary of the attitude errors resulting from 
both the approximations of Section 4. 2 and the covariance errors predicted 
m Section 4. 3. 

The roll, pitch, and yaw bias errors were increased to 0. 05° (l<r) 
m subsequent analyses of the volume, for the two gyro systems, for the 
purpose of conservative error analyses. 


Table 4-1 Summary of ADS Attitude Errors 


Criteria/Error Sources 

Error Source 
Magnitude 

Attitude Error 

Absolute Attitude Determination 





Pitch 





Radiance model 

l 4 km (iff) 

0 

022 

(ler) 

ESA bias 

0 028 (Iff) 

0 

028 

(iff) 

Limit cycle (aee below) 

— 

0_ 

00? 

_Q£i 



RMS = 0 

036 

(IO 

Roll 





Radiance model 

1 4 Vtn (ler) 

0 

022 

(iff) 

ESA bids 

0 028 (Iff) 

0 

028 

(Iff) 

Limit cycle {see below) 

— 

0_ 

oos 

JM 



RMS - 0 

036 

(iff) 

3 aw 





Gyro drift rate bias 

0 06 /hr (lo-) 

0 

029 

(iff) 

Gyro noise 

0 10 /hr (Iff) 

0 

ooo 

(iff) 

ESA noise (transient error 

0 04 (Iff) 

0 _ 

025 

Alsl 

at 300 sec) 


RMS = 0 

038 

(iff) 

Relative Attitude Determination 





Pitch limit cycle 





Reaction wheel quantization 

8 pul /rev 

0 

000 

(iff) 

Ratio of inertias 

2% 

0 

002 

(i<r) 

Reaction wheel misalignment 

4 arc min max 

0 

000 

(iff) 

S/C disturbance torques 

1 x 1 0" 5 ft lb 

0 

006 

(iff) 

Tape recorder momentum 
change 

0 007 ft 1b sec 

0 

000 

(iff) 

S/C products of inertia 

4 slug ft^ 

0_ 

002 

(iff) 



RMS ^ 0 

007 

(iff) 

Roll limit cycle 





ESA noise (transient error 
at 300 sec) 

0 04 (Iff) 

0 

oos 

(iff) 

Gyro noise 

0 10 /hr (Iff) 

0_ 

ooo 

JLL£l 



0 

005 

(iff) 
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5. TOTAL SYSTEM ACCURACY 


5 1 INTRODUCTION AND SUMMARY 

This section documents the results of the overall ERTS system 
geometric and photometric error analysis investigations Only the error 
analyses of the selected data processing configuration are presented here 
Tradeoff studies supporting the final configuration selection are docu- 
mented m Volume 17. A complete description of the selected data proces- 
sing configuration is given m Volume 14 The data processing is divided 
into the following four modes 

Bulk I 

No geometric correction for MSS or RBV, hard copy imagery pro- 
duced on the LBR with photometric adjustments for gam, gamma and 
shading. 

Bulk II 

Digital RBV reseau measurements used to geometrically correct 
RBV image m PPR (precision photo restitutor), attitude determination 
data used to geometrically correct MSS image m PPR 

Precision I 

Similar to bulk II but all geometric corrections m MSS and RBV 
performed digitally, also attitude data may be improved by ground truth 
measurements 

Precision II 

A "special request" mode, can include image enhancement, cos- 
metic corrections, etc 

Error analyses were performed to obtain overall system mapping 
and photometric accuracies for both the MSS and RBV payloads. The chief 
results of the geometric analysis are given m Table 5-1 The photometric 
analysis results are summarized m the MTF curves and the radiometric 
analysis of Section 5-6 Tradeoff studies which supported the final data 
processing configuration selection are presented separately m Volume 17 
Only the error analyses of the selected configuration are presented here 
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Table 5-1 Geometric Error Analysis Summary 



Picturing Centering Error 
Due to Attitude Determination 

Internal Consistency 
Error 

Total Error 

Processing Mode 

MSS 

RBV 

MSS 

RBV 

MSS 

BBV 

Bulk I 

(No geometric correction) 

11 000 ft (3<r> 

8 000 ft ( 3cr> 

9 000 ft |3c) 

12 000 ft (3<r) 

11 000 ft (3c) 

14 000 ft (3c) 

Bulk II 

(In PPR correct RBV 
using digital reseau 
measurements and MSS 
using attitude data from the 
spacecraft) 

9 000 ft (3c) 

8 000 ft (3c) 

1 900 ft (3c r) 

{2 700 ft <3c) 
one gyro case} 

880 ft (3cr) 

{880 ft (3 c) 

one gyro case} 

9 000 ft (3c) 

8 000 ft (3c) 

* 

Precision 1 

(Digitally correct RBV 
using digital reseau 
measurements digitally 
correct MSS using attitude 
data (a) ground truth (b) 
PPR correlation of RBV 
with MSS to improve esti 
mate of yaw bias or (c) 
both 

Precision 2 
Special request mode 

I" 

11 


(a) 1 500 ft to 

1 900 ft (3cr) 
depends on 
ground truth 
location 
accuracy 

(b) 1 500 ft (3c) 

(c) 1 110 ft to 

1 400 ft (3c) 

320 ft (3<r) 

(a) 1 500 ft 
to 1 900 ft 
(3<r) depends 
on ground 
truth location 
accuracy 
lb) 8 000 ft (3c) 
(c) I 500 ft to 
1,900 ft (3c) 

8 000 ft (3<r) 
930 ft (3c) 
using ground 
truth 


The accuracy numbers are errors m the digital tape output and do not include LBR or film processing errors 


5 2 MSS ERROR ANALYSIS 
5 2 1 Summary 

In digital processing of the MSS picture, satellite attitude and 
ephemeris data is used to correct for geometric distortions of the pic- 
ture. The MSS scans out, on an angle versus time basis, the full 
100 nmz X 100 nmi picture To understand the relationship between the 
scanner image and scanner attitude see Figure 5-1 The vehicle moves 
along the orbital path as shown during the first scan At During this 
time period a corresponding series of terrain points are swept out by 
the scanner, these points correspond to the intersection of the instan- 
taneous scan vector j3 with the earth's surface. In digital data proces- 
sing this is done by a knowledge of vehicle attitude and satellite motion 
during the scan The latitude -longitude is computed by the intersection 
of p, originating from the satellite location, with an "adjusted" geoidal 
surface representing the mean altitude of the terrain within one MSS 
picture. Because of attitude and ephemeris errors the wrong terrain 
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SATELLITE ORBIT 



Figure 5-1 

MSS PICTURE Geometry During Scanning 


points will be computed Obviously earth rotation must be considered m 
order to accurately reconstruct longitude from the data 

Volume 17 contains a discussion of the picture distortions produced 
by various error sources The principal effects, however, are pitch and 
roll bias (producing a MSS picture centering error), yaw bias (producing 
a skewing of the picture), and vehicle attitude time variations that produce 
internal consistency errors within the picture Due to computing time 
involved, the precision ground point is not computed for each of the 11 
million picture elements, instead, only some 25 reference points are 
computed within the picture, the rest are computed by a "pseudo-reseau" 
technique similar to that used m the RBV data processing 

Equations have been derived for the absolute position erroi sensi- 
tivities These equations give the sensitivity of picture erroi s to geo- 
metric and sensor error sources A computer program for MSS error 
analysis has been developed using these sensitivities m a covariance 
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analysis of both, absolute mapping and internal consistency errors. 
Internal consistency errors are of more importance than absolute posi- 
tion errors since, if the internal consistency is good, one ground truth 
point can establish the absolute position of all other points m the 
100 nmi X 100 nmi picture. Let Ax, Ay represent the cross-scan and 
along-scan absolute position errors, respectively, and let ej, . , e n 

represent error sources 1 through n. Then 


> 

Ax 


r 9Ax 
9e x * ' 

9Ax 

3e 

n 


e l 

• 

Ay 

J 


9Ay 
aej • 

3Ay 

3e 

n 

J 


♦ 

e 

n 

^ j 


or, m matrix form 


£ = M e 

— i i — i 


(5-2) 


^ is the vector absolute position error at time t^, _e is the error source 
vector at time t^ and M i is the sensitivity matrix 

Internal consistency is defined as the accuracy with which one 
ground point can be located relative to another. Thus, the internal con- 
sistency, or relative position error, can be found by differencing the 
absolute errors at two different points t^ and t m the MSS picture (taking 
into account the error-correlation between the two points). One of the 
two points, the reference point, is normally taken to be at the center of 
the picture Thus, if A£ represents the internal consistency error vector 



M e - M e 
i-i j -j 


(5-3) 


a covariance analysis of (5-2) gives, for the absolute error 


C 


£ 


M C 
l e 

l 



( 5 - 4 ) 
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where for the relative error, a covariance analysis gives 



C c +C c -M <e e^> 

h 1 - 1 


T T T 

M - M < e e > M 
1 3 - 3-1 1 


(5-5) 


From (5-5) it should be noted that with the exception of uncorrelated 

noise, the internal consistency errors go to zero as t^ 1 (that is, the 

two points approach each other within one scan line) It should also be 

noted from (5-5) that the buildup of internal consistency errors with 

time separation t - t depends critically on the time as well as the spatial 
J 1 ip 

correlation of the error sources ce^^e >. It is expected that the attitude 
determination errors will be exponentially correlated m time as well as 
being cross- correlated m roll and yaw. Terrain variations will be expo- 
nentially correlated m terms of their spatial separation (distance) m the 
picture. MSS alignment errors will be constant m time with zero 
ensemble average. Sensor quantization will be assumed as white noise 
(uncorrelated from one point to the next) 


Thus it can be see's. that the internal consistency error analysis is 
critically dependent upon the assumed error- correlation model The 
fmcil MSS error analysis results are shown m Table 5-2, where the rela- 
tive errors are given as measured from picture point No 13 (center of 
the picture) or picture point No 25 (corner of the picture) As can be 
seen from the totals, the 3- sigma absolute error is on the order of 
1.5 mile while the 3-sigma internal consistency errors are 2700 feet as 
measured from one corner to the diagonally opposite corner Table 5-2 
also shows that ground truth data reduces the absolute errors to the order 
of the relative errors (about 1500 feet 3- sigma for ground truth data on 
the order of 200 feet 1 - sigma) 


The results given m this table apply to the precision digitally 
processed "pseudo reseau" points which can then be used m con 3 unction 
with a quadratic distortion interpolation function to remove the picture 
distortions caused by spacecraft motion and attitude limit cycling. 
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Table 5-2 MSS Error Analysis Summary 


Summary (Errors 

at #1) 

1 - Sigma 
North 

(Cross-scan) 

(feet) 

1 -Sigma 
East 

(Along- scan) 
(feet) 

Attitude Determination 5 ^ 

Absolute 

2680 

2738 

(No ground truth) 

Relative (#13) 

493 

284 

Relative (#25) 

880 

385 

Attitude Determination 

Absolute 

381 

356 

(With ground truth) 
of 200 ft lcr 

Relative (#13) 
Relative (#25) 

316 

282 

Terrain Variations 

Absolute 

33 

230 

<r, = 200 ft 

h 

Relative (#13) 

33 

230 

1/(3^ < 25 miles 

Relative (#25 

51 

324 

Quantization Only 

Absolute 

66 

66 

(230 ft maximum) 

Relative (#13) 

94 

94 


Relative (#25) 

94 

94 

Sensor & Alignment 

Absolute 

200 

365 

+ Quantization 

Relative (#13) 

126 

385 


Relative (#25) 

136 

449 

Totals^ 

Absolute 

2682 

2756 

(no ground truth) 

Relative (#13) 

505 

485 


Relative (#25) 

890 

592 


^Numbers are for the 1-gyro configuration. For the 2-gyro configura- 
tion the 1- sigma total errors are reduced to 2690 ft. absolute and 
653 ft relative (as measured from picture point #25) 


5 2. 2 Computation of Sensitivity Matrix for the MSS Imagery 
5 2. 2. 1 Summary 


The equations for the absolute position error sensitivities are 
presented. These equations give the sensitivities as a function of satel- 
lite position and scan location at a particular time or interest t These 
equations m addition to those of Section 5 2 2 3 and the orbit prediction 
equations were programmed on the time shared computer for use m con- 
ducting MSS systems error analysis The computer program will 
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evaluate both absolute and relative (internal consistency) position errors 
Appendix B contains a derivation of these equations 

5 2. 2 2 Satellite Reference Coordinate Axes 

Compute the altitude vector (measured up from subsatellite point 
to ERTS position along normal to geoid) 


** 1 


Kx , Ky , , 
+ K- + r+K i + 


a 2 /b 2 Kz 
1 + a 2 /b 2 K 


(5-6) 


Where i, _j, k are inertial coordinates with k along the polar axis and l m 
the equatorial plane along the vernal equinox direction, x, y, z are the 
ERTS position vector components m Earth fixed axes, a is the equatorial 
radius and b the polar radius of the geoid, and K is determined by itera- 
tion of the following equations 


K , . = K -f AK , K n = 
n+1 n n 0 


. V 


2 2 2 
x + y ~ + z - a 


AK 

n 


S 

n_ 

d S 
n 

d K 


n 


where 


S = 
n 


2 2 2 
x + y b 

2 + ~~1 

(1 + K r a 


n 


,“ 2 . 2 .2 

(b /a + K n ) 


d S 
n 

d K 
n 


,, 2 , 2 , ,, 2 2 

2(x + y ) 2b z 

(1 + K n ) 3 a 2 (b 2 /a 2 + K n ) 3 


(5-7) 
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_ 7 

The iteration proceeds until AK <10 , this corresponds to an error m 

h less than 10 feet. The satellite reference coordinate axes (which 

represent the gyrocompass attitude of ERTS with zero limit cycle angles) 

are given by equation sets (5-8), (5-9)* and (5-10) for _e , e , e 
■ " z x — y 

respectively. 

for the e yaw axis vector 
— z 


e 

-z 




(5-8) 


Z 


x 




q 

A 


where 


q = 


(1 +K) 

~~Z — 2 — — z 

(b /a + K) 


A 


2 





i 


For the e roll axis vector 
—x 


where 


e = (X , X , X ) 

— X ' x y’ z 


X 

x 



Z 

= 

B 


X 



A 


1_ 

A 


(x x + y y + q q) 


q 


(b 2 /a 2 - l)z K 1 + K 

(b 2 /a 2 +K) 2 (b 2 /a 2 +K) 



2 2 

Z + Z 
x y 


+ Z 


2 

z 


(5-9) 
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(5-10) 


x A x 


*7 — * zz. y ~ y y y — 

*-* ~ ° " A ’ ? A * ^ rr ~ 


y A y 


„ q a q 


x .2 A 
A 


* A C A ’ * A 


For the _e pitch axis vector 


Y = Z X - Z X 
x y z z y 


Y = Z X - Z X 
y z x x z 


(5-11) 


Y = Z X - Z X 
z x y y x 


5. 2. 2. 3 Scan Vector Computation 

The MSS scan vector p^ m satellite roll, pitch, yaw coordinates is 


= Col (0, sm (3^ cos pj 


(5-12) 


The scan single p j is given by 


P = -6 +|| (t - t ) 
K i K At ' i s' 


(5-13) 


Where (3 is the maximum excursion of the scan line, t g is the time of the 
start of the current scan line, At is the scan time (for one line) The 
start of the current scan is 


t = n At + t 
si o 


n = 

i 


Int 


t - t 

_i 

At 


o 
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Where 


n = number of current scan line 
1 

t - start of scan 
o 

At = scan time 

The scan vector is then transformed to Earth fixed coordinates to give 
the LOS vector u^ 

\ = [ T ] O] [ L ] ^ (5-14) 

where the coordinate transformation matrices [rj , [e] , [L] are 


[t] 


cos wt 

1 

- sin wt 

i 


sm wt 

l 

cos wt 


l 


0 


0 


0 

0 

1 


(Earth rotation) 


( 5 - 15 ) 


where time t is measured from zero when the Greenwich meridian pass 
through the vernal equinox line, and w is the Earth’s rotation rate 

The [E] matrix transforms from satellite (normal to surface) 
reference axes to inertial coordinates 


X Y Z 

XXX 


fEl = X Y Z 

L J y y y 


X Y Z 

z z z 


( 5 - 16 ) 
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The roll, pitch, yaw limit cycle matrix is (Small angle approximation) 


M = 


- a 


a 


a 


a 


a 


- a 


(5-17) 


5.2 2 4 Ground Point Computation 

The vector point on the surface of the Earth, where the LOS vector 
u from the satellite intersects the geoid, is the vector sum of the satel- 
lite position plus pu^ or 


= PP X + £ 


(5-18) 


The equation of the geoid is, m terms of ground point position 

P = col (P , P , P ) 

— i ' x y z ; 


? „2 a ^2 2 

P + P + =^P - a =0 

x y b 2 z 


(5-19) 


or, substituting (5-18) into (5-19) 


(p u x + x) 2 + (p u + y) 2 l^lp^ + z) 2 
y b 


= 0 


(5-20) 


Equation (5-20) constitutes a quadratic m the satellite -ground point range 
p that can be solved to give 


P = 


-C 


2 



(5-21) 
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Where 


_ 2 2 a 2 

C, = u „ + u + — y u 
1 x y , 2 z 


C 2 = +y U + -J Z u 2 ) 

b 


_ 2 2 a 2 2 

C, = x + y 1 — y z - a 


Equation (5-21) and (5 -18) give the ground point position P corresponding 
to time point t j on the scan. 

5. 2. 2 5 Computation of Sensitivity Matrix 

The sensitivity matrix M is computed numerically by introducing 
unit error sources into u^, the satellite ephemeris _r , or geoid param- 
eters a and b (corresponding to terrain variations) Thus, each error 
source e^, . . . , e^ has associated with it an erroneous ground point 

P. , , P 

— li — ni 

The position error vector is (for unit error sources Ae^, . , Ae.^) 


^?lx = fix - ?i A *1 


(5-22) 


AP = P - P for Ae 
— m — ni —i n 
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The cross scan Ax and along scan Ay components of the absolute position 

I. x 

error are 


Ax, = AP, . e 
li — li —x 


= • % 


} for Ae 


1 


Ax = AP e 
m — ni —x 


> for Ae 


Ay = AP • e 
m — y 


n 


The numerical sensitivities are then 


Ax 


li 


Ae., 


Ax 


m 


* 'Ae 


n 


M = 


Ay^ 

Ae, 


Ay. 

* A? 


m 


n 




(5-23) 


(5-24) 


5. 2. 3 Mathematical Description of Internal Consistency 
Errors m the MSS Imagery 

5 2 3 1 Summary 

Mathematical expressions giving the internal consistency errors m 
the MSS imagery are derived for various models of the error source cor- 
relation Both time and spatial correlations are considered and will be 
included m the ERTS error analysis program The assumed time corre- 
lation critically effects the growth of the internal consistency position 
errors in that the various time correlation models result m internal 
consistency errors ranging from zero to greater than the absolute position 
error. Section 5.2.4 shows an assignment of error models and 
error sources for use m the internal consistency analysis 
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5.2 3. 2 Effect of Correlated Noise on ERTS Internal Consistency Errors 

Let Ax, Ay represent the cross- scan and along- scan absolute 
position errors, respectively. Let e^, . . , represent error sources 

1 through n Then 



(5-25) 


or 

| = M e (5-26) 

2 i i — i 

where 

|j = col (Ax(t i ), Ay(t i )) 

e = col (e,, . . . , e ) @ t = t 
-l '1 n 7 l 

9 li 

(matrix of sensitivities) 

1 -i 

Internal consistency is defined as the accuracy with which one 
ground point can be located relative to another ground point Thus, the 
internal consistency, or relative position error, can be found by differ- 
encing the absolute errors at two different time points ^ and t^ m the MSS 
picture (talcing into account the error- correlation between the two points) 
One of the two points, the reference point, is normally taken to be at the 
center of the picture Thus, if A| represents the internal consistency 
error vector. 

A£ = £ - £ = M e -Me (5-27) 

— -i -3 i-i J—J 
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The mean and covariance of the absolute error 


ii are 


<£ > = M < e > 
-i 1 —i 


C 


6i 


M C 
i e. 



(5-28) 


where 

C > =2x2 covariance of absolute error 

5 i 

C = nx n error source covariance 

e i 

Equation (5-28) does not depend on the time-behavior of ^ in any sense. 
The mean and covariance of the internal consistency errors are (from 
equation (5-27)). 


<A£> = M <e > - M < e > 
— 1 -! 3 -J 

<A|> = <£ x > - <£ 3 > 


(5-29) 


The covariance is (assuming the have zero mean values) 


C = M <e e T > M T + M <e e T > M T - M <e e T > M T - 
Ag l-i-i i J -3 -J j i -i -j j 

(5-30) 

T T 

M <e e > M x 
J ~J -i i 


or 


C Ae = C f. + C„ - M <e e?> M T - M <e e T > M T (5-31) 

A§ i -i-j j J ~J -i i 

From equations (5-29) and (5-31), note that the internal consistency 
errors go to zero as t — t (that is, the two points approach each other 

1 j 

within one scan line) 
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<A£> 


(5-32) 


= 0 



1 


for t = t 
i J 


Equation (5-32) applies for all error sources with the exception of the 
sensor resolution. Thus, it is obvious that internal consistency errors 
are critically dependent on the time- correlation of the error sources 

5. 2. 3 3 Effects of Three Types of Correlated Noise 

First, consider errors that are constant over the whole picture but 
have zero ensemble average 


C. £ = <{M e - Me) (Me - Me ) T > 
A£ i-i j—y i—i J —y 


(5-33) 


= (M x - ivy <e x e^> (M i - M^)' 


for all t 


t 

J 


Second, consider errors that are constant within one scan line and 
uncorrelated from one scan line to the next. 


and 



(5-34) 



/ 
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Third, consider white noise 


C At = M <e e T > M T + M <e e T > M T for all t , t (5-35) 

1-1-1 1 J -J -J j i j ’ 

hi the first case, equation (5-33), internal consistency errors are 
proportional to the change m the sensitivities from t i to ty and go to zero 
as t^ approaches t^. in the third case, equation (5-35), internal consis- 
tency errors are the statistical sum of the absolute error covariance 

matrices for t and t . 

1 3 

5.2.3 4 Exponentially Correlated Noise Model 

Consider a single error source e^, then equation (5-27) becomes 

A£ = M e - M e (5-36) 

— 1 1 J 1 

Where is now a 2 x 1 sensitivity vector Then, if e is exponentially 
correlated noise'’ 


<e e > 
i J 



(5-37) 


From (5-37) note that as the time separation of the ground points increases, 
their error source cross correlation decreases exponentially The covar- 
iance of the internal consistency error is 


C. £ - M <e^> + M <e^> - M <e e > - 

A§ ill JJl i 1 1 1 

(5-38) 

M <e e > 

111 i 


* Exponentially correlated noise can be thought of as white noise passed 
through a first order filter. 
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or, substituting equation (5-37) into (5-38) 


C A > = <r 

e 


M^" - e‘ (3 | t j " *i I (M x M^)j (5-39) 


Note that for t^ = t^, equation (5-39) gives a zero covariance of internal 
consistency errors. If the time difference is long compared to the corre- 
lation time (white noise case) 


'A£ 


= 0 


for |t^ - tj » -i (5-40) 


Equation (5-40) is simply the sum of the absolute error covariance 

matrices C t and C t 

^1 

Another limiting case of interest is for t - small compared to 
the correlation time 




2 r 

= a 

e L 


(M - M ) (M - M )• 
1 j ' 1 


+ 6 t - t 
F J i 


(M M 
i J 


+ M M 
J 


?>] 


(5-41) 


for p It 


J 


t « 1 

l 


Equation (5-41) gives a square root of time increase of internal con- 
sistency error if the sensitivities are locally constant M = M 

1 J 


c a 4 = 2 % { Pl‘j ' ‘>1 M i M i T } 


(5-42) 


5 2. 4 Error Analysis Results 

The MSS error analysis program was run for the noise model m 
Table 5-3 as a function of the 25 picture locations of Table 5-4. The 
simulation program computes the absolute and relative error covariance 
matrices as a function of picture location The 1- sigma North and East 
position error variances are given m Figures 5-2 through 5-10 for 
various error sources and reference point locations 
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Table 5-3 


Noise Models for Error Analysis 


* 


Error Sources 

Numerical Value 
1- Sigma 

Description 

Attitude Determination 
Bias Error (Radiance 
Mod & Gyro Drift) 

Roll = Pitch 
= 0.05° 

Yaw = 0. 05° + 
0. 707 roll 

Constant error source with 
zero ensemble average 
Roll and pitch uncorrelated 
but roll and yaw correlated 

Attitude Determination 
Time Varying 

Roll = Pitch 
= Yaw = 0 006° 

Time -Exponentially corre- 
lated with time constant 
= 30 sec. 

Alignment 

0 01° max 

Constant- Uncorrelated 

Resolution 

0 00146° 

Ap Random within one scan- 
line 

Scan Rate Variation 

2 psec (3 <r) 

AT Constant within one scan- 
random from one scan'to the 
next 

Dwell Time Variation 

3% 

Ap White noise 

Sync Variation 

40 ft max 

Ap Constant over one scan, 
random from one scan to 
next 

Scan Nonlmearities 

0 2% (3a-) 

Ap Random m one scan 

Terrain Variation 
from Geoid 

er = 2000 ft. 
n 

i/p n = 12, 000 ft. 

Spatially- Exponentially 
correlated (in terms of 
distance) 

Ephemeris Errors 

Negligible immed- 
iately after a state 
vector update 

Time and spatially 
correlated 


Figure 5-2 gives the absolute and relative 1- sigma errors for 
attitude determination error sources only with no ground truth data for 
25 picture points For example at point No 1 the 1- sigma North and 
East absolute errors are 2680 feet and 2738 feet respectively, and the 
1-sigma North and East relative errors (with respect to the center of the 
picture) are 493 feet and 284 feet respectively At picture point No 13, 
as can be seen from Figure 5-2, the relative errors are zero (see equa- 
tion 5-32). The relative error at picture point No 3 is 280 feet x 280 feet, 


1-gyro error model 
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Table 5-4. Location of Picture Points for MSS Error Analysis 


Picture Point Latitude Longitude 


1 

45° 

44' 

56. 0" 

-99° 

46' 

35 0 

2 

45° 

41' 

9 1" 

-100° 

22' 

35 7 

3 

45° 

37' 

12. 3" 

-100° 

58' 

19. 2 

4 

45° 

33' 

4.5" 

-101° 

33' 

58. 3 

5 

45° 

28' 

44. 0 

-102° 

9' 

45. 9 

6 

45° 

19’ 

57 5" 

-99° 

43' 

22 2 

7 

45° 

16' 

12. 2" 

-100° 

19' 

7 0 

8 

45° 

12' 

17. 1" 

-100° 

45' 

34 8 

9 

45° 

8' 

11 1" 

-101° 

29' 

58 5 

10 

45° 

3' 

5 2. 7" 

-102° 

5’ 

30 7 

11 

44° 

54' 

58. 8" 

-99° 

40' 

13.4 

12 

44° 

51' 

14. 9" 

- 100° 

15' 

42. 7 

13 

44° 

47' 

21 5" 

-100° 

50' 

55 2 

14 

44° 

43' 

17 4" 

-101° 

26' 

3 7 

15 

44° 

39' 

0 9" 

J 

V— ' 

O 

o 

r 

20 7 

16 

44° 

29' 

59. 7" 

-99° 

37' 

8.6 

17 

44° 

26' 

17. 3" 

-100° 

12' 

22 7 

18 

44° 

22' 

25. 6" 

-100° 

45' 

20. 1 

19 

44° 

18' 

23. 2" 

-101° 

22' 

13. 7 

20 

44° 

14' 

8.6" 

-101° 

57' 

15 9 

21 

44° 

5' 

0 3" 

-99° 

34' 

7.5 

22 

44° 

1' 

19 4" 

I 

o 

o 

o 

9' 

6. 7 

23 

43° 

57' 

29. 2" 

-100° 

43' 

49 5 

24 

43° 

53' 

28. 6" 

-101° 

18' 

28 4 

25 

43° 

49’ 

15 9" 

-101° 

53' 

16. 1 
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NO 5 

NO 4 

NO 3 

NO 2 

NO 1 

2742,2679 

2702,2671 

2678,2678 

2671,2700 

2680,2738 

498,285 

348,280 

280,280 

344,280 

493,284 

NO 10 

NO 9 

NO 8 

NO 7 

NO 6 

2741,2679 

2701,2670 

2677,2677 

2670,2699 

2680,2737 

409,45 

203,20 

0,0 

203,20 

403,44 

NO 15 

NO 14 

NO 13 

NO 12 

NO 11 

2740,2678 

2700,2670 

2677,2676 

2670,2698 

2679,2736 

409,45 

203,20 

0,0 

203,20 

403,44 

NO 20 

NO 19 

NO 18 

NO 17 

NO 16 

2738,2677 

2699,2669 

2676,2676 

2669,2697 

2678,2735 

451,204 

290,210 

209,210 

292,209 

459,214 

NO 25 

NO 24 

NO 23 

NO 22 

NO 21 

2737,2677 

2698,2668 

2675,2675 

266B, 2697 

2678, 2734 

491,283 

343,280 

280,280 

347,280 

496,284 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(1 


0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 


Figure 5-2 

MSS Absolute and Relative Errors for Attitude Determination Error Sources 
Only No Ground Truth Errors with Respect to Picture Point No 13 
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Figure 5-3 

MSS Absolute and Relative Errors for Attitude Determination Errors Only 
with Two Ground Truth Points on Opposite Sides of Map Relative Errors with 
Respect to Picture Point No 13 
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Figure 5-4 

MSS Absolute and Relative Errors for Attitude Determination Errors Only. 
Relative Errors with Respect to Picture Point No 2 5 
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Figure 5-5 

MSS North, East Absolute (Relative) Errors for Terrain Error Only, 
<7^ = 2000 ft , Errors With Respect to Picture Point No 13 
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Figure 5-6 

MSS Absolute and Relative Errors for Terrain Error Sources Only. Relative 
Errors with Respect to Picture Point 'No 25. 
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Figure 5-7 

MSS Absolute and Relative Errors for Sensor and Alignment Error 
Sources Only Relative Errors with Respect to Picture Point No- 13 
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Figure 5-8 

MSS Absolute and Relative Errors for Sensor and Alignment Error 
Sources Only Relative Errors with Respect to Picture Point No 25 
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Figure 5-9 

MSS Absolute and Relative Errors for All Error Sources. Relative 
Errors with Respect to Picture Point No 13 
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Figure 5-10 

MSS Absolute and Relative Errors for All Error Sources. Relative 
Errors with Respect to Picture Point No 2 5 
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since this point is directly below the spacecraft the error is primarily 
due to time varying pitch and roll errors Note that the errors at points 
No 3, 8, 18 and 23 vary as the square root of the time difference from 
point No. 13 (this was demonstrated by equation (5-42) The north com- 
ponent of internal consistency error at the corners is primarily due to 
the yaw bias error (because the sensitivity to yaw is maximum at the 
corners and zero at the center of the picture) A pictorial representation 
of the relative position error ellipse is also shown m Figure 5-2 

Figure 5-3 shows the absolute and relative errors for attitude 
determination error sources only, assuming that the yaw, pitch, and 
roll error biases are calibrated by means of two ground truth points 
(1- sigma = 200 feet) at opposite ends of a scan line The first thing to 
be noted from Figure 5-3 is that ground truth data makes the absolute 
and relative errors of comparable size (instead of the order of magnitude 
difference m Figure 5-2). It should be noticed, however, that internal 
consistency error is not much reduced over Figure 5-2 (without ground 
truth). 

Figure 5-4 gives the effects of attitude determination error without 
ground truth, with the internal consistency errors referred to picture 
point No 25 The maximum relative error is now increased to 1- sigma 
880 feet north, 385 feet east at the corner of the picture (point No 1) 

Figure 5-5 gives the absolute error and relative error with respect 
to picture point No 13 for spatial- exponentially correlated terrain varia- 
tions of 2000 feet 1- sigma Since the center of the picture was taken as 
the reference point, the absolute and relative errors are equal (the sensi- 
tivity to terrain variations being zero at this point) The error ellipse, m 
this case, is degenerate with the error being along the scan line only. 

Also, it should be noted that with point No. 13 as reference the spatial 
correlation has no effect on the relative errors 

Figure 5-6 gives the effects of terrain variation with point No 25 
(the corner of the picture) as reference. In this case the correlation 
distance does effect the internal consistency However, as long as the 
correlation distance is less than the distance between picture points, the 
error behaves essentially as white noise This accounts for the fact that 
relative error is */2 times the absolute error at the edges of the picture 
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Figures 5-7 and 5-8 give the effects of sensor and alignment error 
sources only. A picture quantization of 230 feet (equal to the resolution) 
is included m the sensor error sources, this is why the relative error 
does not go to zero at the reference point 

Figures 5-9 and 5-10 give the total of all error sources m Table 5-3 
(without ground truth data) The 3- sigma absolute error is on the order 
of 1 5 nmx which easily meets the 2 nmi requirement. 

All the data given here are for the 1-gyro configuration If the 
2-gyro configuration is used, the 1- sigma internal consistency error is 
reduced to 653 ft (from 890 ft m Figure 5-9) 

5 2 5 Pseudo-Reseau Correc tion 

In MSS Bulk I, no geometric corrections are made, and the picture 
is represented m the scan angle ( P x ) versus time (t^) plane MSS preci- 
sion processing involves finding the mapping from the ((3^, t ) plane into 
the ( 0^ cj> ) latitude -longitude plane by means of a knowledge of a small 
number of reference points To do this, 0^ and d> are represented by the 
polynomials 


6 

i 






5 

V J 


(5-43) 
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^ + <V t i) ,(P r P 3 ,2,( Vtj )2 


The distortion vectors a and b can then be computed from (5-43) by a 
knowledge of the precision ground points of six picture elements (obtained 
from attitude and ephemeris data). ' The latitude and longitude compo- 
nents of picture distortion are the differences between the interpolated 
values (5-43) and the precision computed values. If the six points are 
chosen as m Figures 5-11 through 5-14, namely points (1), (3), (5), (11), 
(21) and (25), the modeling errors go to zero at these points and reach a 
maximum between pseudo-reseau points. Figures 5-11 and 5-12 give the 
latitude -longitude distortion errors m radians versus picture location and 
were obtained via computer simulation The errors reach a maximum of 
7 74x10 ^ rad latitude at point No 22 and a maximum longitude erroi of 
-3 936 x 10 5 rad at point No 17. 

Translated into distance errors this is 162 ft North and 5 82 ft East. 
Figures 5-11 and 5-12 also apply for the case of constant attitude errors 
and rates (that is to say that the quadratic distortion polynomial accurately 
models the effects of these errors) Figures 5-13 and 5-14 give the 
modeling errors for a maximum limit cycle of 0. 01° / sec at 0. 1 cps The 
maximum distortion is m Figure 5-14 at picture point No. 19, and is 
1690 ft** 


r-, 's 
b l 


5 


(5-43) 


Volume 17 contains a more detailed discussion of the MSS distortion 
correction algorithms and associated throughput 

•J 

More complicated algorithms have been developed that give even smaller 
distortions The important thing is that the distortions be small compared 
to the interval consistency of the reference points. 
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Figure 5-11 

MSS LATITUDE Distortion vs Picture Location 
Limit Cycle = 0 
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Figure 5-12 

MSS LONGITUDE Distortion vs Picture Location 
Limit Cycle = 0 
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Figure 5-13 

MSS LATITUDE Distortion vs Picture Location 
Limit Cycle is 0 01°/sec, 10 sec period 
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Figure 5-14 

MSS LONGITUDE Distortion vs Picture Location 
Limit Cycle is 0. 0l°/sec, 10 sec period. 
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5 3 RBV GEOMETRIC ERROR ANALYSIS 


The ERTS RBV camera system consists of three RCA return beam 
vidicon cameras which simultaneously photograph the same image m 
three corresponding spectral bands A pattern of fiducial marks, or 
reseau, is interposed between each camera and the scene to provide a 
reference for estimation of the RBV camera distortion The images 
taken by the three cameras are read out sequentially and transmitted to 
the ground data handling system, where each of the res eau points 
positions are determined Based on these positions, the camera dis- 
tortion is determined and the image subsequently corrected, either 
digitally or via analog (Precision Photo -Res titutor). This section will 
determine the image geometric distortion correction accuracy as a 
function of the various system error sources This will include deter- 
mination of the error sources, modeling of the error sources and 
processing algorithms, and evaluation of the algorithm performance 

The RBV system geometric errors can be divided into three dis- 
tinct groups pre-reseau errors, or errors which are caused by space- 
craft attitude and ephemeris, camera alignment, and terrain 
irregularities, and are thus undetected by the res eau measurement, 
post-reseau/pre-processmg errors, which include the camera distortion 
and data lmk/tape recorder errors, and can be estimated through the 
res eau image distortion by the digital processing, and post- digital proc- 
essing errors, including the geometric distortion correction errors and 
the image recording and processing errors The image distortions 
caused by each of the three groups of errors are independent and thus 
may be evaluated separately The post-digital processing errors are 
discussed m Section 5 4 and the post-reseau/preprocessmg errors are 
discussed m Section 5,3 2 

5 3 1 Pre-reseau Errors 

Image distortions due to spacecraft attitude and ephemeris errors 
for the RBV differ from those for the MSS m that the attitude error rates 
do not affect the geometric distortions of the RBV, other than through 
smear, due to the fact that the entire RBV image is exposed at one time 
Keystoning and distortions related to the earths curvature and spacecraft 
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attitude, however, must be considered. To correct for these distortions, a 
linear transformation from spacecraft attitude and ephemeras to geo- 
metric image distortion is determined as a function of position in the 
image Assuming small-angle approximations for roll 0 R , pitch e , 
and yaw 0^, and independence of their effects on image distortions, it 

can easily be shown that apparent x and y distortions 6,6, respectively, 

x y 

are given by 
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where 6 H is spacecraft altitude error, and 6 x q , 6y Q are ephemeris 
errors The errors of interest here are internal consistency errors, 
those which will change distances between points m the image These 
can be expressed as 
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O' = 


H 


4 2. 2 2 22 2 

y or R' } * x y <r p + y 


H 


For o = 05° = and c r = 200 ft , the specifications on the attitude 

p R H 

determination system and ephemens data, the maximum error is 

(3a- ) = (3 o- ) =: 141 ft 

v y x max 

with a correlation coefficient of one, giving (3cr) = 200 ft 

This is a maximum deviation, occurring only at the picture corners 
and diminishing greatly at the interior of the image For example, at the 
picture edge, halfway between two adjacent corners, the 3 or deviation is 
only 141 ft Implementation of this transformation is easily accomplished 
in conjunction with the camera distortion correction algorithm, as dis- 
cussed m the following section 

5 3 2 Reseau Distortion Estimation Analysis 

The images generated by the RBV cameras will contain geometric 
distortions which are primarily a result of imperfect scanning by the 
electron beam of the RBV Data link and tape recorder geometric errors 
will be negligible with respect to the camera errors In the camera, an 
optical system focuses an image on a faceplate, which is scanned by an 

t 

electron beam and converted to a voltage tram to be transmitted to and 
demodulated on the ground Image distortion may arise through non- 
linearities in the scanning beam path or scan rates To obtain an esti- 
mate of this distortion, a known reseau grid (Figure 5-15) consisting of 
81 evenly spaced fiducials m a 9x9 array is placed m front of the camera 
focal plane By measuring the apparent reseau distortion, the camera 
distortion parameters may be estimated 

A comprehensive trade-off study was performed to select the best 
algorithms for estimating the distortion parameters Selection was 
based on a computer study of functional and statistical estimation errors 


5-39 



12 70 »005 


4 4 4 4 
4 4 4 4 
4 4 4 4 


— H [—2 80* 

| TYPIC4 


TYPICAL 
4 4 4 4 

4 + + + 


+ A7YHCAL 
(SEE DETAIL A) 


-i 4 + — 4 + — 4—4 — 4 — 4 — 4 

I 

4 4 4 + 4 

12 70 i 0 05 | + + + + | 

4 4 4 4 4 

I 

»4 4 4 4 4 


j. 

12 7040 05 


4 4 4 4 

4 4 4 4 

4 4 

4 4 4 4 


- 12 70 40 05-H 
* AIL DIMENSIONS Alt IN MILLIMETERS EXCEfT AS NOTED 


25 MICROMETERS 




■ 25 MICROMETERS 


700 MICROMETERS 


| 200 MICROMETERS 

DETAIL A 


25 MICROMETERS 


, B TYPICAL 
4 (SEE DETAIL I) 


25 MICROMETERS 

_L 


250 MICROMETERS 


r+ 


250 MICROMETERS 


TT_ 


DETAIL B 


200 MICAOMETEK 


Figure 5-15 

RESEA U ARRAY Format for RBV 

and on computer timing requirements for each algorithm The results 
of this study are discussed in Section 2. 5 3, Volume 17 Only the error 
analysis results for the recommended system are presented here 

Typical sources of distortion may be modeled as follows First, a 
coordinate system is defined, centered at the image plane center and with 
coordinate axes x and y. Position displacements at a point (x, y) due to 
distortion are given by 6x(x, y) and 6y(x, y) 

An image centering shift can result from incorrect deflection 
voltage bias or misalignment and causes distortion terms of the form 

6x = a Q , 6y = b Q 

The specifications call for a maximum centering error of 1 percent full 
scale. 


5-40 



Image size error can result from incorrect deflection gam and 
cause distortion terms of the form 

fix = a^x , 6y = b.y 

The specifications require these terms to be less than 1 percent full 
scale. 

Skew distortions can arise from sources such as misaligned 
deflection coils and may be represented by terms of the form 

fix = a^y , 6y = b^x 


According to specifications, these terms should be less than 5 percent 
full scale 

Pincushion distortion can arise from magnetic field imperfections 
and gives terms of the form 


fix = x^a 3 (x 2 + y 2 ) + a 4 (x 2 + y 2 ) 2 + J 

Sy = y£bj(x 2 + y 2 ) + b 4 (x 2 + y 2 ) 2 + j 


Symmetric tangential distortion can be caused by the combined 
effect of unwanted radial electric fields and the axial focusing magnetic 
field and gives terms of the form. 

fix = y^a 5 (x 2 + y 2 ) + a^fx 2 + y 2 ) 2 + . j 

5y = x^bg(x 2 + y 2 ) + b 6 (x 2 + y 2 ) 2 + j 


5-41 



Specifications call for the maximum "image distortion" to be less 
than 1 percent full scale Since specifications for skew, size, and 
centering are quoted separately, an "image distortion" is assumed 
referring to all other distortions including pincushion and tangential 
distortion 

The total effect of all of these sources of distortion becomes 

2 2 2 2 2 

6x = a Q + ajX + a 2 y + a^x (x + y ) + a 4 x (x + y ) 

, 2 , 2 , . , 2 , 2,2 

+ a 5 y (x + y } + a 6 y (x + y ) 

6y = b Q + b x y + b 2 x + b 3 y (x 2 + y Z ) + b 4 y (x 2 + y 2 ) 2 

+ b 5 x (x + y ) + b^x (x + y ) 

where only the first two terms m pincushion and tangential have been 
included. No specifications could be found which give the stability of the 
distortions from picture to picture This list of distortion terms is not 
exhaustive, possibly other sources of distortions will be found m the RBV 
Because these terms are not uncommon, they can serve as a specific 
test of proposed distortion correction schemes It can be assumed that 
the distortions are of the form* 

x + Sx = f(x, y) + v 

where f is a smoothly varying function and v is random spatial noise 
A similar equation exists for y distortion The function f, which will be 
called the functional distortion, contains terms arising from all smooth 
distortion sources, such as those listed above The spatial noise terms 
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can arise from sources such as sync-jitter, beam jitter, and 
high frequency distortions It is assumed that the noise is spatially 
uncorrelated so that 


E(v v ) = 6 cr 
1 j lj N 


at points i and j The values of Sx^^ for the reseau points are given as 


6x - f(x , y ) + v + m , i - 1, 2, 
i i ] v 1 1 


, N. 


R 


til 

where represents the measurement error for the i reseau point It 
is assumed that 


2 

E(m m ) = 6 <r and E(v m ) = 0 
1 J ij m li 

The reseau displacements are fit by an interpolation algorithm to 
determine the distortion elsewhere 

The 9x9 reseau grid is subdivided into 16 squares with three reseau 
points to a side and a total of nine points each A two-dimensional 
quadratic polynomial is then fitted to the distortions at the nine points of 
each subregion The form of the biquadratic model for x (or y) distortion 
is 


6x 


■i i 

1=1 j=i 


a x 1 y-* 
ij 


where the a are parameters to be determined In matrix format, this 
becomes 


6x = A(x, y) a 
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where 


1 2 1 2 ' 2 1 2 ' 2 2 ' 
xy i x ly i xy ixyixy 
ill ii 




r i ' i 

A (x, y) = lixlyi 

L ' I l 

r 

-M = [ a oo ,a lo> a ol’ a ll’ a 2o’ a o2 ,a 12’ a 21’ a 22] 


The actual distortion 6x also contains other terms such as first-order 
? 7 

pincushion x(x + y ) a^ o and spatially uncorrelated noise v as discussed 
above Thus, 


5x{x, y) 



N(x. 




-M 


-N 


+ v(x, y) 


where 6 x(x, y) is x-distortion at (x, y), a^ are the unmodeled distortion 
parameters, and N(x, y) are the coefficients of a^ m the distortion 
The distortion can be measured at each of the nine reseau points m each 
subregion and arranged as follows 


6x 


[A N] 


2M 

% 


+ V 


where 


5x = 

'Sxfxj, y 1 )‘ 

, A - 

’A(x 1 , yj)' 

, N = 

“N(x 1 , y 2 ) * 

, V = 

"vfxj, y 2 ) + m x _ 


6x(3c 2 , y 2 ) 


a ( x 2' 


N(x 2* y 2> 
♦ 


v(x 2 , y 2 ) + m 2 


, 5x( vy 


A(X N- V 




/ ( V y N ,+m N. 


where for biquadratic model, N = 9. 
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The distortion estimation equation, based on the biquadratic model, is 
given by 


6x(x,y) = A(x, y ) a = A(x,y)A _1 6x 

which leaves a distortion estimation error e(x, y) 

e(x,y) = e E (x, y) + e M (x, y) 

where e E (x, y) is the statistical error caused by the noise v(x, y) 

€ E ( x >y) = A(x,y) A _1 v -v(x,y) 

and e^(x, y) 1S the functional error caused by neglecting the terms a.^ m 
the distortion estimation 

€ M (X ’ y) = K(x> y * -N 

where 


K(x, y) = A(x, y) A 1 N -N(x, y) 


The term y m j= E (x, y) includes reseau point measurement errors which 
are small with respect to the term v(x, y) and can be root- sum- squared 

A computer study was made of the performance of this algorithm 
with respect to the variance of the spatially-uncorrelated distortions and 
reseau measurement errors 


E|v ( x, yj) 


= <r 


all reseau points (x^, y ) 
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The results of this study are shown in Figures 5-16 and 5-17 
Figure 5-16 presents the ratio of the standard deviation of the statistical 
error due to reseau point distortion to c The grid represents one 
quadrant of the image plane, the center of the image appearing m the 
upper left corner Entries appear twice per reseau point spacing m the 
horizontal direction and four times per reseau point spacing m the verti- 
cal The entries corresponding to reseau locations are encircled It 
can be seen that the statistical error m the far corner of the image has 
uncertainty of 5 219 times the uncertainty of the reseau location due to 
measurement and random camera distortion The sensitivity of the 
residual distortion estimation error to unmodeled distortions is shown 
m Figure 5-17. The format of the grids m this figure is the same as that 
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Figure 5-16 

BIQUADRATIC INTERPOLATION 
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Modeling Biquadratic Interpolation 
(First-Order Pincushion) 
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Figure 5-17 

BIQUADRATIC INTERPOLATION 



m Figure 5-16 The entries in the first grid represent the sensitivity 
to the low-order effects of pincushion -type (second-order) distortions and 
those m the second grid the sensitivity to higher-order effects Thus, if 
1 percent full scale first-order pincushion is present m the image, 

1 04 percent of this distortion or 01 percent full scale will remain m the 
image corner Similarly, 5 52 percent of second-order pincushion dis- 
tortion will remain m the same corner 

If a (small) number of reseau points cannot be detected, their positions 
are estimated in a similar manner in this event, the matrix A(x, y) is 
given by 


f 1 I 1 { 2 1 2 1 2 1 2 • 3 > 3~| 

A(x, y) = It xiy i xy ix ly l xy ixylx iy 

L i I i I I l t I I J 

and the parameter is determined for the entire image using all 81 
reseau points m a least squares algorithm The appropriate estimation 
equation is thus 


6x (x , y ) - A(x, y) (A^A^A^Sx 

J J 

where Sx and A contain all the available reseau distortion measurements 
and coefficients. Here the matrix A is a 81 x 9 The corresponding error 
is e(x, y) as before, with 

e E (x, y) = A(x, y) (A T A) _1 A T v -v(x, y) 

€ M = (A T A) ~ 1 A T N -N(x, y)) a N 

Again, these errors were evaluated over a quadrant of the image plane 
and the results are shown in Figure 5-18 The format of these grids is 
the same as m the previous figures The sensitivity of to first-order 
pincushion is identically zero since first-order pincushion is included m 
the model The errors of interest m this figure are the reseau point 
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Figure 5-18 Global 3rd-Order Polynomial 


errors The standard deviation of the statistical reseau location estima- 
tion error is a maximum of 633 cr^ at the outermost reseau point, where 
the modeling error reaches a maximum of 5 91 percent of the maximum 
second-order pincushion m the image The maximum errors are sum- 
marized in Table 5-5 together with a typical conversion to picture element 
(pixel) magnitudes It can be seen that errors are much greater m the 
image region outside the reseau pattern, indicating the advisability of 
extending the reseau pattern as close to the image edge as possible 

It should be stressed that these are worst case numbers which apply 
when only one image frame is available If the functional distortion is 
fairly constant over subsequent frames, then an average of the coefficients 
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Modeling Error {Second-Order Pincushion) 
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Figure 5-19 

GLOBAL THIRD-ORDER Polynomial 

obtained for past frames will lower the statistical uncertainty by \f~n, 
where n is the number of frames averaged If severe beam-pulling is 
an important factor m reseau position determination, coefficient averaging 
will greatly dimmish this source of error After enough frames are 
averaged, the functional error will dominate and thus be the determining 
factor in accuracy calculations 

The pre-reseau distortion correction algorithm discussed m 
Section 5 3. 1 is incorporated into the camera distortion correction algorithm 
with negligible increase m processing time as follows The parameters 
m the biquadratic interpolation algorithm are augmented by the 
following vector 


= HG 
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Table 5 5 Maximum Error Magnitudes for RBV Processing 



Biquadratic 
Interpolation 
(Entire Image) 

Biquadratic 
Interpolation 
(Interior to Res eau) 

Reseau Location 
Estimation 

Statistical 

or 

5 219 

1 000 

633 

cr 

R 

Sensitivity 
to 1st Order 
Pincushion 

1% 

2% 

0 

Sensitivity 
to 2nd Order 
Pincushion 

5 5% 

85% 

5 9% 

RSS J 

1 ± 8 pixels 

2 ± 1 5 pixels 

5 ± 1 pixel 


Assuming 1 % max 1st order pincushion 


2% max 2nd order pincushion 
1/2 pixel ( 0125%) random distortion (lcr) 

1/2 pixel { 0125%) res eau location error (l<r) 


for x distortion and 




He R ; e y |H &H ;s e pi 0 ;^ 9 r« 


°] 


for y distortion The modified estimation algorithm is then 


6x t (x, y) = A{x, y) = A(z, y) (A~ l 6x - b^.) 


Here, as before, a = (A ^ 5x - b ) must be calculated only once per 

* — M — — —x 7 c 

interpolation subregion, thus increasing the computational load negligibly 
with respect to that for calculating 5x (x, y) at each of the 17 million image 
points {x, y) The pre-reseau and camera distortion errors are 
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independent and may thus be root- sum- squared, giving a total internal 
consistency error within the reseau grid of 360 ft (3<x) under the assump- 
tions of Table 5-5 

5 3 3 Ground Truth Incorporation Analysis 

An important performance parameter of the RBV (and MSS) system 
is the residual absolute and relative image errors when image UTM 
coordinates (ground truth) are known for a number of image points and 
are used to reduce the effects of attitude and ephemeris uncertainties of 
the spacecraft This section will detail a computer analysis of the 
accuracy of the RBV ground truth incorporation algorithm 

It is pointed out m Section 5 3 1 of this volume that the image error 
at a point (x,y) m the image due to spacecraft attitude /ephemeris devi- 
ations can be expressed by* 


6x 

o 

6y o 


■fix' 


H 0 xy (H 2 +x 2 ) 

1 

K 

w 

!* 

r 


1 


5y 

“ H 

0 H (H 2 +y 2 ) ^ 

xH y 


6H 


or 


6 = M(x, y) a 

where x and y are along-and across -track distances from image center, 

H is nominal altitude, 6x , 6y , and SH are ephemeris deviations, and 

o o 

0-^> 0p and 0^ are roll, pitch, and yaw deviations The parameters a 
are measured by the attitude determination and ephemeris systems as 
with specified random errors w 


a J* = a -f w 
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where w is zero mean with covariance W ' Thus, based on these systems 
alone, the attitude/ ephemens pre-reseau errors are given by 

= M(x, y) w 


as shown m Section 5 3 1 

However, if N ground truth measurements 6^‘ at (x , y ) l = 1, 
, N, are available, the estimate of a can be improved through the 
following least-squares approach Define 


-» r 


1 = 

V 

, M = 

1 

£ 

1 — • 

jwal 

,£ = 

‘*1 


.V 


M(x N* y N } 


6' 

. N . 


Then 


6_ = M a 

The measurements 6^ are corrupted by residual internal con- 
sistency image processing errors, geodetic ground control point errors, 
and image reading errors, which are essentially uncorrelated between 

widely-spaced ground-truth points If these errors are given by v , with 

2 

mean zero and variance c r^j or 

v = [v r ..,v N ] , E^yy j = o-y I 

Then 

6 = M a + v 
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The least squares estimate of a^ has an error covariance 


J 

a 


T 1 -1-1 

(M L M ~ + W ) 

°V 


giving a distortion error covariance at (x, y) of 

J 6 = M(x, y) J a M T (x, y) 


The least squares estimate is given by 


6(x, y) = 


M(x, y) J M 
a 


T 1 


(6 - M a ) 


The equations for the pre-reseau distortion estimation error covar- 
iance Jg were solved via digital computer for several system configurations; 
one- and two-gyro attitude determination with digital and PPR-distorti on 
correction The appropriate specifications are seen m Table 5-6 


Table 5-6 Pre-Reseau Error Source Magnitudes 


Ephemens 

Attitude 


Roll -yaw correlation 


Ground Truth Sq T 


^xo’ V’ = 200 ' M 

^ROLL’ '’’PITCH = 05 ° ( ltr ) 

o-y aw = 05° dual gyro (l«r) 

073° single gyro (ler) 

0 dual gyro 
47 single gyro 
312' (digital correction) 

350' (PPR correction) 
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The ground truth error standard deviations are obtained by- 
root- sum- squaring the geodetic map position error, the residual relative 
distortions in the corrected photographs, and the image point location 
measurement errors 

Tables 5-7 and 5-8 present the residual pre-reseau errors both 
without ground truth and with one to three ground truth points The case 
of one ground truth point assumes the point at image center, that of two 
points at opposite image corners, and that of three points at image 
corners The residual errors are shown at the image center and the far 
image corner, thus giving the minimum and maximum of the random 
errors These errors are then summarized with the remaining system 
errors m Table 5-14 It is noted for the single gyro case that roll and yaw 
are correlated Errors m y caused by these angles will add m one-half 
the picture and subtract in the other The tables contain an average of 
the y error m the four corners 


Table 5-7 Digital Processing Ground Truth Incorporation Errors 


No of Ground 
Truth Points 

Absolute Errors* 

Relative Errors 

Picture Center 
x y correl 

Picture Corner 

x y correl 



Picture Corner 
x y correl 


0 

2674 

2674 

0 

2730 

2730 

- 07 

43 

43 

-1 


1 

310 

310 

0 

463 

463 

52 

21 

21 

-1 

Single Gyro 












2 

218 

218 

- 02 

262 

262 

25 

20 

20 

-1 


3 

182 

182 

03 

246 

246 

42 

20 

20 

-1 


0 

2674 

2674 

0 

2714 

2714 

- 01 

43 

43 

-1 


1 

310 

310 

0 

412 

412 

41 

21 

21 

-1 

Dual Gyro 












2 

218 

218 

- 02 

259 

259 

24 

20 

20 

-1 


3 

182 

182 

04 

241 

241 

41 

20 

20 

-1 


^All numbers are lcr values in feet 
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Table 5-8 PPR Processing Ground Truth Incorporation Errors 




Absolute Errors 

T* 


Relative Errors 


Picture Center 

Picture Corner 

Picture 

Corner 

No of Ground 
Truth Points 

X 

y 

correl 

X 

y 

correl 

X 

y 

correl 

0 

2674 

2674 

0 

2730 

2730 

- 06 

43 

43 

-1 

1 

Single Gyro 

347 

347 

0 

489 

489 

45 

21 

21 

-1 

2 

244 

244 

- 02 

292 

292 

24 

20 

20 

-1 

3 

200 

200 

- 01 

256 

256 

34 

20 

20 

-1 

0 

2674 

2674 

0 

2713 

2713 

- 01 

43 

43 

-1 

1 

347 

347 

0 

441 

441 

35 

21 

21 

-1 

Dual Gyro 

2 

244 

244 

- 02 

287 

287 

23 

20 

20 

-1 

3 

204 

204 

03 

267 

267 

40 

20 

20 

-1 


All numbers are lo - values m feet. 
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5 4 POST DIGITAL PROCESSING ERROR ANALYSIS 


After the distortion parameters for the RBY have been estimated 
by the digital computer m the bulk II processing, the parameters are 
used off-line with a precision photo -re stitutor (PPR) to generate an 
image corrected for this distortion In precision mode I, the RBV and 
MSS corrections are applied digitally within the computer and then 
reproduced via a laser beam recorder (LBR) and processed photograph- 
ically The digital distortion correction procedure will contribute less 
than one pixel error, thus indicating the LBR and photo processing as 
the limiting functions m image accuracy This section will provide an 
error analysis of the PPR and the LBR Geometric errors due to 
photographic processing and film shrinkage is estimated at 182 ft 
(1 -sigma) and is independent of the remaining errors 

5 4 1 Precision Photo- Re sti tutor Geometric Error Analysis 

The PPR is used m two distinct modes In the bulk II processing, 
digital computer distortion correction data for each image is input to the 
PPR, which generates an image corrected for geometric distortion In 
the precision I mode, the red RBV channel is used as a control photo to 
which the red channel of the MSS is correlated and registered The 
amount of correction necessary to register the images is used to further 
estimate the yaw bias of the spacecraft attitude and thus decrease the 
MSS error substantially The operation of the PPR is covered m detail 
elsewhere m this report, it is the intent of this section to investigate 
geometric image distortions left by the device 

The PPR system must be controlled online from a computer All 
transformation can be pre-programmed and is available for each frame 
to be procured, but the position of the input and control images are not 
known accurately until they are measured at the PPR Minor adjustment 
to the control equations must be made after frame fiducial marks are 
measured To reduce the online data rates, it may be convenient to 
drive the PPR through a computer interface buffer which is loaded with 
data for one printing strip The PPR will be operated automatically 
after all film rolls are loaded 
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Two primary computer requirements exist for the PPR m its 
operation on RBV data The first set of calculations is needed prior to 
the insertion of the photograph into the PPR, so its results must be 
available from online storage The second set of calculations is the rota- 
tion, translation, and scale change of the first set of data to the coordinate 
system on the film actually m the PPR at that time This set of data is 
determined from three LBR points measured m the PPR m combination 
with the output of the first calculations. The new data is returned to drive 
the PPR The PPR will read the three LBR points, pause for computation, 
and begin to expose the nine strips of the photograph Each strip begins at 
a known center position (x, y), known scale (M), and known rotation (X) 
Increments of the strip are exposed by moving the PPR slit with differ- 
ential incremental values {dx, dy, dM, da) stored in a 4 by 147 buffer 
The computer can proceed strip by strip, calculating values at the begm- 
nin g> inserting these values m the buffer, and calculating again at the 
beginning of the new strip. If all data is known, the computer can initially 
calculate for the entire nine- strip photograph 

The residual errors are based on the relative and absolute errors 
that exist after error is removed by reseau calibration, lens calibration, 
and satellite attitude from ground truth using 1 25, 000 map geodetic 
positions and are also based on the internal characteristics of the PPR 
Since the bulk II images and the master RBV for MSS processing are 
always accomplished on the PPR, the figures for errors will be those 
given m Table 5-9 If the maximum excursion of the satellite is 0 8 
degree and the Earth rotation effect is 20, 000 feet m one-half picture, 
a size of maximum excursion on Ax, Ay, scale, and rotation was 
obtained The maximum distortion per column is 3 degrees, 54 minutes 
rotation, 12 6 percent increase m scale, and a maximum x or y displace- 
ment of 39, 361 feet In a single element of the column (9 x 147) the 
maximum rotation is 23 minutes, scale change is 1 4 percent, and Ax 
incremental change is 229 minutes, Ay incremental changes are not 
significant for sizing the PPR True scale is increased to 105 7 percent 
as seen m the first increment in a column and this decreased to 93 1 
percent at the bottom Appropriate limits have been included m the 
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Table 5-9 Geometric Residual Errors 


Output of Master Copy- 

Geographic Residual Relative Residual 
Error, Feet Error, Feet 

PPR performance 

(0 8- or 0 4-mch slit) 206 150 206 150 


design of the PPR The residual errors are a function of the extreme 
excursions of the PPR error removal capability 

Dynamic errors are those that occur during exposure time with 
resultant loss of restitutor resolution The dynamic error allowance is 
budgeted equally m the mam servos to result m an error allowance of 
20 micrometers per servo This error is added (rms) to the error that 
results from averaging from top to bottom of the 0 5- by 0 8-mch slit 
to obtain the table for the worst condition of distortion (Table 5-10) 

These resolution losses are cut down for the 0 4-mch-wide slit 

The residual errors are the rms sum of relative errors, absolute 
(geographic) errors, LBR errors, film shrinkage errors, and PPR 
induced errors If the RBV master is produced m the PPR and if all 
error removing capabilities are used, the figures given in Table 5-10 
result The errors indicated as PPR errors m Table 5-11 are the 
registration errors m the multispectral image registration 

5 4 2 Laser-Beam Recorder Geometric Error Analysis 

An accurate error analysis requires a sufficiently detailed physical 
description and configuration of the device being investigated Because 
LBR technology is still relatively new and undergoing changes, the 
information available from manufacturers (especially those that relate to 
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Table 5-10 Overall Dynamic Errors 



Servo 

Dynamic 

Total 

0 05 by 0 8 

Total 
0 05 by 4 


K rotator, ft 

112 

70 

132 

118 


Z motion {scale) ft 

100 

78 

127 

108 


x and y input, ft 

100 

70 

122 

106 


x and y motion, ft 

— 

35 

35 

18 



Table 5-11 

Comparison of Errors in 
and Control Images 

Master RBV 


PPR Slit 



0 8 by 0 05 inch 

0 04 by 0 05 inch 

Correlation error 

±60 microns 

±30 microns 

Setting error 



Ax, Ay 

±30 

±30 

6 ± 0 1 degree 

±34 

±17 

M ± 0 3 percent 

±30 

±15 

Linear interpolation 

and blur ±20 

±10 

RSS 

±83 microns 

±49 microns 


±276 ft 

±163 ft 


error sources and their valuation) is meager and quite often sketchy 
For this reason the following assumptions and restrictions are placed 
upon the analysis 

• The error estimates presented are based on published informa- 
tion currently available from the manufacturer and as a result 
of visits and telephone calls 

• Where error data was not forthcoming, estimates were made 
on the basis of experience and analysis of the relevant error 
source 

• The list of error sources presented in the body of the text does 
not attempt to present a comprehensive compilation Instead 
the list shows the most significant error sources on the basis 
of our current understanding of the physical system 
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For the purpose of this analysis, the output image from the LBR 
is defined to he the exposed film prior to its processing This 
suggests that the error analysis does not include any photo- 
metric and geometric changes m the output image (such as film 
shrinkage, film gamma variations) due to photoprocessing 

• The LBR error analysis assumes the same device will be used 
for all three types of input - RBV and MSS type image data as 
well as the precision processed data coming from the NDPF 
computers 

The recorder analysis task selected LBR as the technology most 
suitable to the ERTS requirements at the present time While specifica- 
tions will be prepared for an LBR film recorder m this phase of the ERTS 
program the selection of a unit by a particular manufacturer is left for 
Phase D In the absence of a specific unit, the approach m the error 
analysis has been to categorize the LBR units or models currently 
available m the field and from these categories determine the errors to 
be expected from their component functions and then arrive at composite 
error estimates 

• Swept Beam Scan Type 

This type of LBR generates its scanning motion by sweeping its 
laser beam through reflection off a rotating mirror or prism 
The RCA and Ampex LBR ! s and the CBS Lab LBR use this kind 
of scanning motion Two variations are shown, where m one 
case the spot forming function is performed before and m the 
other case after the scanning mechanism The movement of the 
beam m the longitudinal direction is provided by a film 
transport 

• Swept Film Scan Type 

This type of LBR generates its scanning motion by the use of 
film wrapped around a rotating drum over which a fixed orienta- 
tion laser beam is projected Translation of the beam longi- 
tudinally is provided by a transport mechanism which moves the 
laser modulator and spot forming system parallel to the scanning 
drum. 

The error analysis presented investigates each of the two types and 
uses the results to determine a composite LBR error estimation In the 
absence of a selection of a particular LBR, the results really indicate what 
is currently realizable m LBR performance m the field without its being 
identified with a particular unit 

Each functional block of the LBR has been analyzed m Table 5-12 
listing the error sources with the estimates of their magnitudes given m 
Table 5-13 The error sources are listed which cause spatial distortions 
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m the output image and therefore affect the accuracy of locating areas of 
interest on the ground and their relative distances and orientations They 
are divided into two groups 

1) Scan Coordinate Error (S) - These relate to errors along 
the direction of the scan as it is generated on the output 
image 

2) Longitudinal Errors (L) - These errors relate to the error 
sources m a direction perpendicular to the scan motion 

The error estimates prepared m Table 5-12 are indicated m 
Table 5-13 m summary fashion for every one of the LBR functional 
blocks This manner of presentation permits comparison of the error 
estimates for both basic types of LBR - the swept beams scan type and 
the swept film scan type It is evident from the tables that the differences 
are minor so that it is reasonable to characterize both types of LBR by 
a single set of error estimates This is precisely the intent of the use 
of a so-called "composite" LBR model for error analysis It should be 
noted that for the sake of being conservative, whenever there was a 
choice among estimates, the higher values were used 

Errors for all the LBR functional blocks are combined m rms 
fashion when they are unrelated errors or otherwise by direct addition 
to arrive at the overall combined errors for each LBR type. These two 
sets of estimates are m turn combined to obtain an overall LBR set of 
error estimates They are as follows 

• Longitudinal error 0 05 percent 

• Scan coordinate error 0 02 percent 

The geometric error is the resultant of the errors m the two coordinates 
5 5 OVERALL GEOMETRIC ERRORS 

The overall geometric distortion errors for the RBV and MSS 
imagery are shown for bulk I, bulk II and precision I processing m 
Tables 5-14 and 5-15 Precision processing absolute errors are shown 
with and without ground truth data Applicable uncertainty values are 
shown m the Table All entries are 3 sigma values Precision pro< ess- 
mg also includes digital tape output accuracy, l e , accuracy befoi e LBR 
and photo processing distortion 
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Table 5-12 List of Geometnc Error Sources 


Item No 
1 0 
1 1 

1 2 

2 0 
2 1 

2 2 


3 0 


3 I 


4 0 


4 1 


1 Laser 

1 Variation m beam position due to vibration 
in laser mounting - geometric effect on 
image 

2 Temp changes may cause beam variations 
due to the large heat input probably 
associated with a laser 

2 Modulator 


1 Variation m beam position due to mechani- 
cal vibration of the modulator 

2 Temperature changes may cause variations 
m beam position since a modulator is 
inherently inefficient m converting input 
energy 

3 Spot Forming System 

This is the optical system which operates to form 
the spot It may consist of several lenses and 
reflectors, but for this block it excludes any 
scanning optics if used 

1 Variation is spot position on the film due 
to vibration of the spot forming system 
and/or temperature changes 

4 Swept Beam-Scan Mechanism 


This functional block includes the scanning motor 
1 Longitudinal errors 

a Spacing - optical error due to varia- 
tions m prism or mirror surface 
causing a variation in line spacing 

b Pairing error, if multiple beams are 
used 

c Sequential skew error due to axial 
misalignment of mirror axis with 
motor rotational axis 
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Table 5-12 List of Geometric Error Sources (Cont) 


Item No 

4 2 2 Scan coordinate error 

a Jitter along scan - caused by motor hunting 

b Shift of the scan line along the scan This 

error is caused by the optical error between 
mirror faces - 1 e , azimuth alignment 
errors 

c Fixed skew error due to multiple input beams 
departing from a common plane 

d Scan non-linearity 

(1) Scan velocity error, in contrast to jitter, 
is a longer term error It is also caused 
by variations m motor speed (see sketch 
below) 

(2) Non-uniform film distance to the rotational 
axis will cause deviations m spot velocity 
at the film (see sketch below) This may 
T^e due to rotational axis motion 


5 0 
5 1 



rotational axis 


film plane 


Scan velocity error is due to oj changing 
with 1 = l 1 ^ 

Film distance error is due to 1 changing 
with to - co 1 


5 Film Transport 

1 Table speed deviations from desired values 
result m longitudinal error This is a 
rate error, and depending on whether it is 
open or closed loop would determine its 
characteristics 
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Table 5-12 List of Geometric Error Sources (Cont) 


Item No 


5 2 

5 3 

5 4 

6 0 

6 1 

6 2 

7 0 
7 1 

7 2 


2 Transport banding error due to the same 
cause as 1 The error is manifested as a 
bunching or separation of lines and is 
expressed as a percentage of the line to 
line dimension 

t 

3 Longitudinal jitter - due to drive motor 
random hunting or drive mechanism high 
speed variations The error may be a 
function of the operating speed 

4 The errors shown m the above sketch may 
also be due to transport misalignment 
having similar effects 

6 Film 

Only exposure of film considered, not its photo 
processing 

1 Non-uniform film thickness or tension 
causing slip may result in non-uniform 
film speed resulting m 1 ) longitudinal 
error and 2) scan coordinate error 

2 Wave-like distortions (buckling, etc ) of 
film may result m non- linearities m both 
coordinates 

Swept Film Scan Mechanism (Scanning Drum) 

1 Longitudinal error 

a Scan jitter caused by axial end-play of 
the drum, improper alignment of 
the axis Temperature changes can 
also produce this effect 

2 Scan coordinate error 

a Jitter along scan - caused by motor 
hunting 

b Shift of the scan line - due to the 

timing error between the start of the 
scan and drum position 
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Table 5-12 Last of Geometric Error Sources (Cont) 


Item No 


8 0 
8 1 

8 2 

8 3 


c Scan non- linearity due to surface 

speed variations of the drum caused by 
rotational speed changes and out-of- 
roundnes s of the drum 

Beam Translation (Laser beam, spot forming transport) 

1 Transport errors due to speed variations 
This can cause striations m the output 
image 

2 Non-uniform drum surface distance to 
beam caused by variations m the movement 
of the beam transport 

3 Misalignment of the beam transport axis 
with the drum axis causing a shift of the 
scan line along the scan giving a biased 
skew to the output 


5-66 



Table 5-13 LBR Error Estimates 


Item No 

Error Source 

Error Estimates 
(N - negligible) 

1 0 

Laser 


I 1 

Vibration 

N . 

1 2 

Temperature Changes 

N 

2 0 

Modulator 


2 1 

j 

Vibration 

■N 

2 2 

Temperature Changes 

.. N 

3 0 

Spot Forming System 


3 1 

Geometric Variations 

N 

4 0 

i > • 

Swept Beam Scan 

Amp ex - 0 1% 

4 1 

Longitudinal Errors 

RCA - 0 05% 



CBS , - 0 01% 

4 2 

Scan Coordinate Errors 

Ampex - 0 02% 



RCA - 0 02% 

CBS - 0 01% 

5 0 

Film Transport 


5 1 

Table Speed Variations 

Included with 4 0 

5 2 

Transport Banding Error 

< 0 01% 
r 

5 3 

Longitudinal Jitter 

Included with 4 0 

5 4 

Misalignment _ 

Included with 4 0 

^ ' t K 

6 0 

Film 

i r 

6 1 

Thickness or Slip 

N 

6 2 

Wave -like Distortions 

N 

7 o 

Swept Film Scan 


7 1 

Longitudinal Errors 

1 o oor% 

i * 

7 2 

Scan Coordinate Errors 

0 001% 

8 0 

Beam Translation 

> .. J 
I » 

8 1 

Geometric Error 

<0 1% overall 
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Table 5-14 R.BV Total Errors 


Bulk I 


Internal 

Consistency 

(3 <r) 


Absolute 
(3 cr) 


• Attitude (6000 ft 960 ft S000 ft 

altitude, 0 35° pitch, 
roll, 0 7° yaw 
errors) 


• Camera distortion 12000 ft 12000 ft 

(2% max ) 


• LBR distortion 

328 ft 


328 ft 

• Photo processing 

distortion 

540 ft 


540 ft 

Total 

12000 ft 

14000 ft 

Bulk II 




» Attitude (600 ft 

altitude, 0 15° pitch, 
roll, 0 21 0 yaw 
3 sigma uncertainty) 

200 ft 


8000 ft 

• Camera distortion 

estimation (internal 
to reseau) 

230 ft 


230 ft 

• PPR correction 

accuracy (0 8 x 0 5" 
slit) 

618 ft 


618 ft 

• Photo processing 

distortion 

540 ft 


540 ft 

Total (RSS) 

880 ft 


8000 ft 

Precision I 


With 

Ground 

Truth 

Without 

Ground 

Truth 

• Attitude (See 
Bulk II) 

200 ft 

630 ft 

8000 ft 

• Camera distortion 

estimation (See 
Bulk II) 

230 ft 

230 ft 

230 ft 
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Table 5-14 RBY Total Errors (Cont) 



Internal 

Consistency 

(3<r) 

With. 

Ground 0 -^ 

Truth 

Without 

Ground 

Truth 

• Digital distortion 

correction accuracy 

100 ft 

100 ft 

100 ft 

• L.BR distortion 

328 ft 

328 ft 

328 ft 

• Photo processing 

distortion- 

540 ft 

540 ft 

540 ft 

Total (RSS) 710 ft 

928 ft 

8000 ft 

' Tape only {not includ- 
ing terms asterisked) 

320 ft 

680 ft 

8000 ft 

J "Ground truth location 
600 ft (3 cr), 2 data points 






Table 5-15 MSS Total Errors 


Bulk I 


• Attitude /sens or 
errors (6000 ft 
alt, 0 35° pitch, 
roll, 0 7° yaw) 

• LBR distortion 


Internal 
Consistency 
(3 a-) 

9000 ft 


328 ft 


Photo processing 540 ft 

distortion 


Total (RSS) 9000 ft 


Bulk II 

• Attitude /sens or 
errors (600 ft 
alt, 0 15° pitch, 
roll, 0 21° yaw) 

• PPR correction 
accuracy (0 8 x 
0 05" slit) 

• Photo processing 
distortion 


1900 ft 
(2700 ft)'"' 1 * 

618 ft 


540 ft 


Total (RSS) 2100 ft 
(2700 ft) 


Absolute 

(3or) 

11000 ft 

328 ft 
540 ft 

11000 ft 
9000 ft 

618 ft 

540 ft 
9000 ft 


Precision I 


With Without 

Ground ’ ’ Ground 

Truth Truth 


• Attitude/ sensor 1500 ft 1500 ft 8000 ft 

errors (yaw bias 

determined by PPR 
correlation of 
MSS and RBV) 

• Digital distortion 600 ft 600 ft 600 ft 

correction accuracy 

• LBR distortion^ 328 ft 328 ft 328 ft 
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Table 5-15 MSS Total Errors (Cont) 


Internal 
Consistency 
(3 <r) 

With 

Ground 2 *" * 
Truth 

Without 

Ground 

Truth 

• Photo processing 

distortion «• 

540 ft 

540 ft 

540 ft 

Total (RSS) 

1700 ft 

1700 ft 

8000 ft 

-i'Tape only (not includ- 
ing terms asterisked) 

1600 ft 

1600 ft 

8000 ft 


■^Ground truth location 
600 ft (3 sigma), 

2 data points 

J J ^Numbers m parentheses are for the one gyro configuration 
If there are no numbers m parentheses, they are the same 
for the one gyro configuration as for the two gyro 
configuration 


5-71 




5 6 PHOTOMETRIC ANALYSIS 
5 6 1 MTF Analysis 
5 6 11 Summary 

The payload of the first earth resources technology satellite will be 
three return beam vidicon (RBV) cameras and a multispectral scanner 
(MSS) These sensors provide pictures of the earth in various optical 
bands for identifying and cataloging features of interest It is important 
to determine the expected radiometric performance of the system, m 
order to properly identify and size its operations 

This section computes the S/N ratios and the modulation transfer 
functions of the return beam vidicon cameras and the multispectral 
scanner These sensors are then characterized by their S/N ratios 
viewing a 200 foot bar pattern and by their limiting resolutions Table 
Table 5-16 lists the analysis steps The specifications of the two sensors, 
listed m Tables 5-17 and 5-18, are taken from the '‘Design Study Specifi- 
cations for the Earth Resources Technology Satellite ERTS A and B, " as 
outlined m Attachments I and II of that document Table 5-19 summarizes 
important system parameters The assumed target is the one listed m 
Section 3 3.1 5 of Attachment H Its properties are extrapolated for use 
m the RBV analysis as shown m Table 5-20 Granted the above assumptions, 
the "DC" S/N ratio may be computed and used to normalize the zero fre- 
quency of the modulation transfer function curves From these curves. 
Figures 5-20 through 5-27, signal to noise ratios for various targets may 
be found. Table 5-21 lists the system SPP/NRMS for extended targets 
Table 5-22 summarizes the sensor performance relating to the hypothetical 
target 

A somewhat separate question concerns the information degradation 
of the data link It need not be configured for the high contrast target 
considered above, as very few such targets exist in nature Instead, a 
data link S/N ratio of 33 dB is probably satisfactory. Table 5-23 charac- 
terizes such a data link Finally, Table 5-24 summarizes over-all 
system performance, for viewing targets of moderately high contrast 
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Table 5-16. Design Approach 



RBV 



MSS 

Highlight/towlight radiance 


Highhght/lowiight radiance 

Exposure time 



Detector area 


Vidicon sensitivity 



Detector responsivity 

Highhghfc/lowhght current 


Highlight /lowhght current 

Assume noise-m-; 

signal limit 


Assume noise-in 

-signal limit 

s pp /n rms 1ow 

spatial frequency 


Spp/N RMS low spatial frequency 

Cross Track MTF 

Along Track MTF 

Cross Track MTF 

Along Track MTF 

(Along Scan) 

(Cross Scan) 

(Along Scan) 

(Cross Scan) 

2 Atmospheric 

Atmospheric 


Atmospheric 

Atmospheric 

Turbulence 

Turbulence* 


Turbulence* 

Turbulence 

2 Optics Resolution 

Optics Resolution 

Optics Resolution 

Optics Resolution 

3 Vidicon Beam 

Vidicon Beam 


NA 

NA 

4 Video Amplifier 

NA 


Video Amplifier 

NA 

5 S/C Limit Cycle 

S/C Limit Cycle 

S/C Limit Cycle* 

S/C Limit Cycle 

Oscillation 

Oscillation 


Oscillation 

Oscillation 

6 NA 

Uncompensated Image 

NA 

Uncompensated Image < 


Motion 



Motion 

7 NA 

Raster Line Spacing 

Sampling Aperture 

Sampling Aperture 

Limiting Resolution 





S/N @ 200 feet 

Same 


Same 

Same 


Negligible 



Table 5-17 Return Beam. Vidicon Camera Specifications 

Resolution (at maximum scene 
highlight contrast) 

Edge resolution (percent of center) 

Signal-to-noise ratio (at 10 TVL) 

Dynamic range 

Grey scale (\/2 transmission 
steps) 

Shading (maximum vertical and 
horizontal) 

Residual image (maximum) 

Horizontal scan rate (lines /sec) 

Number of scan lines 

Read out time (seconds) 

Video bandwidth (MHz) 

Time between picture sets (Sec) 

Exposure time (milliseconds) 

Exposure control (milliseconds) 

Image distortion (maximum) 

Deflection skew (maximum) 

Size and centering shift 
(maximum) 

Spectral bandwidth (nanometers) 
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Table 5-18 Multispectral Scanner Specifications 


Sensitivity and Signal/ Noise Ratio - The sensitivity of the MSS shall be 
adequate to provide high signai/noise (S/N) ratios (peak to peak signal/ 
rms noise) for typical scene reflectances The following S/N voltage 
ratio shall be provided for radiance values corresponding to a 230 ft 
equivalent target 




Bands 



1 

2 

3 

4 

High radiance (high 
contrast) watts cm“2 
ster"! 

(25)(10)- 4 

(23)(10)“ 4 

(13)(10)“ 4 

(15)(10)- 4 

Signai/noise 

156 

120 

60 

71 

Low radiance (high 
contrast) watts cm~2 
ster“ 1 

2 5 x 10 4 

2 3 x 10- 4 

1 3 x 10“ 4 

1 5 x 10“ 4 

Signai/noise 

47 

36 

18 

21 


Band 5 shall provide 1°C noise equivalent temperature difference 
(NEAT) for a scene temperature of 300 °K under normal operating 
conditions 


Modulation Transfer Function - The MTF for the MSS shall be adequate 
to provide a 35% response at a spatial frequency corresponding to 
200 feet in Bands 1, 2, 3, and 4 A 35% response at a spatial frequency 
corresponding to 600 feet shall be provided m Band 5 
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Figure 5-21 
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Figure 5-26 

SYSTEM S/N Performance, MSS Band II, High Contrast Bar Targets 
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Figure 5-27 

SYSTEM S/N Performance, MSS Band III, High Contrast Bar Targets 
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5 6 12 Analysis DC S/N 


The standard equation relating irradiance m the focal plane to 
radiance of the target is 

- TN >C 

H x = 2 

4 (F/#) Z 

= spectral irradiance in the focal plane 

N ^ = apparent spectral radiance of the target 
T = optics transmission 

jr/$ = optics focal ratio 

The resultant signal current for the RBV is 

I s = rSH x 

and for the MSS is 


where 


h = A d RH x 


S = vidicon photo cathode sensitivity (A/j/cm^) 
t s vidicon exposure time (sec) 

R = photomultiplier photo cathode res pons lvity (A/W) 
= area of resolution cell of MSS (cm 6 ) 


So, we may define a variable K such that 

X S = KH x 

K = xS for RBV 

K = A,R for MSS 
d 

If the dominant noise is shot noise (which is a good approximation for 
photomultiplier tubes and for vidicons operated m the ruturn beam mode), 
the noise current is 

X N = 

= Noise current 

e •= Electron charge (coulombs) 

Af = Noise bandwidth (Hz) 


5-81 



As our target is a square wave pattern of radiance, the signal to noise 
ratio is 

S/N = (I 2 - I 1 )/N 

S/N = Peak-to-peak signal to root mean square noise ratio 
I = Highlight current 


X 1 

hT = 


Lowlight current 
Average noise current 




2eAf 2 (I 1 + I 2 ) 


Hence, 


S/N = 


h ■ h. 


yj e ^f y 1 2 + 


S/N = a 


^V\ 


where 




KttT 


4eAf(F /#)' 


These system parameters are listed in Table 5-19 Target parameters 
of apparent radiance and approximate reflectance are listed in Table 5-20 

The results of these calculations are the signal/noise ratios of 
both sensing systems in their different spectral bands when viewing a large 
high contrast target 

5 6 13 Analysis MTF 

The overall MTF is the product of the MTF functions of each 
degradation process These are summarized m Table 5-16 and will be 
considered in that order 
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Table 5-19 System Parameters 


Parameter 

Band 1 

Band 2 

Band 3 

x 

RBV 

0 475-0 575|j, 

0 58—0 68(jt 

0 69- 0 83|x 

A. 

MSS 

0 5-*- 0 6p 

0 6— 0 7(jl 

0 7~ 0 8(i 

T 

RBV 

0 525)1 

0 63|jl 

0.76p 

A. 

MSS 

0 55[j. 

0 65)i 

0 7 5 [jl 

S 

RBV 

20 mA/j/cm^ 

2 

36 mA/j/cm 

7 5 mA/j/cm^ 

R 

MSS 

43 5 mA/w 

22 mA/w 

6 mA/w 

T 

RBV 


12 msec 


A d 

MSS 


{5 84 x 10~ 3 cm) 2 


Af 

RBV 


3 5 MHz 


MSS 


43 1 KHz 


F/# 

RBV 

MSS 


2 66 
3 52 


r r 

RBV 


69% (Spectral) 


i 

MSS 


170% (Spectral) 
(90% (Obscuration) 


a 

RBV 

5 77 x 10 3 

7 74 x 10 3 

3 53 x 10 3 

MSS 

2 94 x 10 3 

2 10 x 10 3 

1 22 x 10 3 



cm w“ ^ / 2 

cm w“ l/ 2 

-1/2 

cm w 


1 Atmospheric Turbulence 

The MTF defining the effect of turbulence within the 
atmosphere in the spatial frequency response of the 
cameras is defined by the following equation * 

Tj(k) = exp (-2ir 2 f 2 cr 2 k 2 ) 

T^k) = MTF function 

k = corner frequency (cycles /mm) 

f = focal length 

cr - seeing strength (radians) 

^Hufnagel, R E , “Random Wavefront Effects, 11 Perkm-Elmer Symposium 
on Modulation Transfer Functions, March 6, 1963 
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Table 5-20 Target Characteristics 


In Band 
Radiance J 

0 5— >-0 6(i 

0 6—0 7[i 

0 7—0. 8p 

0 8—1 l(i 


2 5 

2 3 

1 3 

1 5 


0 25 

0 23 

0 13 

0 15 


'mwatts cm~2 s ter"^ 


MSS target radiance (from NASA-GFSC Specification) 


Spectral 

Reflectivity 

0 55 

0 65 

0 75 

0 95 

_ Clear 
Hl Haze 

91% 

94% 

66% 


92% 

97% 

70% 


Clear 

L °Haze 

5% 

5% 

6% 


0% 

1% 

4% 



Target reflectance (for 60° Solar Zenith Angle) 


In Band 
Radiance 

0 475-0 575(i 

0 58-—0 6 8fjt 

0 69—0 83(i 

Hilite 

2 65 

2 3 

1 75 

Lolite 

0 26 

0 23 

.. I 

0 17 

_2 

m watts cm 

. -1 
ster 




RBV target reflectance (for similar reflectivity) 

Table 5-21 Peak-to-Peak Signal/Root Mean Square Noise 


Sensor 

Band 1 

Band 2 

Band 3 

RBV 

257 

318 

117 


48 dB 

50 dB 

41 dB 

MSS 

126 

87 

38 


42 dB 

39 dB 

32 dB 
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= 0. 30 /altitude m feet 
= 10 radians 

The effect of atmospheric turbulence is negligible as the 
MTF is 1 00 for spatial frequencies up to 1000 cycles/mm 


2 Optics Resolution 


The optics of the two systems are quite different An MTF 
for the MSS optics is supplied so it will be used ^ The 
RBV cameras utilize refractive optics which may be 
characterized as follows 3 


T 


2 





T 2 (k) = MTF function 

k - spatial frequency (cycles /mm) 
k-^ - diffraction limit spatial frequency 


- (XF/#)" 1 


This function of course assumes diffraction limited 
optics 

3 Vidicon Beam Resolution 

This MTF applies only to the RBV cameras It is 
supplied by the manufacturer 4 

4 Video Amplifier 


This MTF is applied only m the direction of scan The 
two systems have considerably different bandwidths and 
only the MSS bandpass characteristic is specified, a 
two pole RC filter 5 For convenience a similar filter 
will be assumed for the RBV 


11 System Design Study, Multi -Spectral Scanner (MSS), 11 Hughes Aircraft 
Company, Santa Barbara Research Center, August 26, 1969 

2 

0 T Neill, E L , Introduction to Statistical Optics, Addison-Wesley 
Publishing Company, Inc , Reading, Mass , 1963 

4 

n Vidicon-2 n Diameter Return-Beam Camera Tube, RCA Developmental 
Type Number C23061A, Ir RCA Preliminary Technical Data Sheet, 
January 1968 

5 Ibxd 2 
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1 


t 4 oo 


1 + k 2 /k 0 2 


T^(k) - MTF function 

k - spatial frequency (cycles /mm) 
k Q ■= corner frequency (cycles /mm) 


For the RBV we have 

k Q = 164 cycles /mm 
For the MSS we have 

k Q =15 cycles /mm 


5 Spacecraft Limits Cycle Oscillation 

The MTF defining the effect of spacecraft limit cycle 
motion on the spatial frequency response of the cameras 
is 


Tj.(k) = sinc(Tr a^k) 

Tp.(k) = MTF function 

k = spatial frequency (cycles /mm) 
a^ = image motion during exposure (mm) 

a r - F 6 r 
D 

F = objective focal length 
6 ■= angular rate (radians /sec) 
t = exposure time 
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The anticipated angular rate of the ERTS satellite is 
1 4 x 10-4 radians /sec per axis The nominal exposure 
time of the RBV cameras is 12 msec For the RBV, 


-3 

aj- = 0 212x10 mm 

The nominal exposure time of the MSS is microseconds 
Hence this term is significant only for the RBV cameras 

6 Uncompensated Image Motion 

This MTF is applied m the direction of flight The MTF 
function of UMC is 


T^(k) = sinc(ir a^k) 

T^(k) - MTF function 

k = spatial frequency (cycles /mm) 

- image motion during exposure (mm) 

F 

a 6 “ A VT 

F - focal length 
A - satellite altitude 
v = satellite subpoint velocity 
t = exposure time 

For the RBV 

a, - 1 09 x 10"^ mm 
6 


As before, the MSS exposure time is so short that its 
MTF is unaffected 
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7 Sampling Aperture 


Raster line spacing and aperture size represent the 
same phenomenon, although raster MTF is applied only 
m the cross scan direction The MTF for both is 

- sm (it a^,k) 

T^(k) - MTF function 

k ~ spatial frequency (cycles /mm) 
a.^ ~ aperture size (mm) 

For the RBV 

_3 

a^ - 5 443- x 1 0 mm 
Far the MSS 

a~, - 5 84 x 10 ^ mm 

This list represents all the MTF's believed to affect sensor 
resolution up to the data down link Aperture correction is believed to 
be undesirable for this application 6 

The MTF of each sensor is 

MTF(k) = nT.(k) 

where the product is taken over the appropriate i Table 5-16 lists the 
appropriate i. The individual MTF's for each sensor are shown m 
Figures 5-20 and 5-21 The appropriate products are shown m 
hhgures 5-22 through 5-27 Table 5-22 lists the S/N ratio for viewing 
200' bars and the limiting resolution These values apply to the infor- 
mation output from the sensor Further calculations are necessary to 
determine the S/N ratio of the data used to reconstruct the images 
6 

Laverty, N P , "Resolution of Film Records of Video Data from ERTS 
Return- Beam Vidicon Cameras, " IOC 7244. 3-23, TRW Systems 
28 Aug 69 
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Table 5-22. Sensor System Performance 




Band 1 

Band 2 

Band 3 



RBV 

MSS 

RBV 

MSS 

RBV 

MSS 

— 







S/N 

200' bar 

Along 

Track 

23 dB 

31 dB 

25 dB 

28 dB 

16 dB 

20 dB 

target 

Cross 

Track 

31 dB 

28 dB 

33 dB 

25 dB 

24 dB 

17 dB 

Bar size 

Along 
Track i 

150' 

120' 

150' 

125 1 

160' 

135' 

S/N=0 dB 

Cross 
Track \ 

130' 

125' 

130' 

j 

130' 

140' 

145' 


5 6 1.4 Data Link 

The S/N ratio and the MTF from the two camera systems will be 
degraded by the spacecraft video tape recorder, die RF link and the 
ground based recorder The degradation can be computed by the 
following equation 


S/N = 




2 


Table 5-23 lists the S/N ratios of the various items of the data 
link and the total degradation The S/N ratio for viewing 200' bars and 
the limiting resolution may be calculated and are presented in Table 5-24 
The MTF' s of the various items m the data link are believed flat in the 
frequency region of interest The photo reproducer similarly will not 
degrade the MTF (as it is down 3 dB @ 9000 TVL) The S/N ratio of the 
finished prints is not well defined and is not generally a useful param- 
eter Radiometric quantities must be measured from the tapes 
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Table 5-23. Data Link S/N Degradation 


Band 1 


Band 2 


Band 3 


Sensor S/N 

SC recorder 

RF link 

Ground 

recorder 


RBV 

MSS 

RBV 

MSS 

RBV 

MSS 

257 

126 

318 

87 

117 

38 

126 

225 

126 

225 

126 

225 

153 


153 


153 



10 

error 

rate 


10 

error 

rate 


10 

error 

rate 


Output S/N 


38 

93 

38 

73 

36 

37 

32 dB 

39 dB 

32 dB 

37 dB 

31 dB 

31 dB 

16 dB 

3 dB 

18 dB 

2 dB 

10 dB 

1 dB 


Table 5-24 System Performance 


Band 1 


Band 2 


Band 3 


RBV MSS RBV MSS RBY MSS 


S/N ratio Alcm? 25 5 dB dB 

200 ' bar Track 


Cross 

Track 


Track 

Cross 

Track 


15 dB 


13 dB 16 dB 



5 6 2 1 Basic Considerations 

The transfer of information through an imaging system, whether a 
photographic, television, or direct viewing optical system, is accom- 
panied by degradation of the information content by the inherent 
characteristics of the system Generally, such systems are designed 
to display an image for direct viewing by a human observer Therefore, 
the observer must be considered as a component of the system 
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The system may be designed to compensate m part for the 
degradation of information caused by the characteristics of the observer's 
eye and brain Image enhancement techniques are used to increase or 
change those characteristics of the input information, as required, to 
ease the task of the photointerpreter m detecting, recognizing, and 
identifying objects from their images m the output display 

The return beam vidicon (RBV) three- camera system and the 
multispectral spot scanner (MSS), used m the Earth Resources Tech- 
nology Satellite (ERTS), depend on energy derived from the Sun and 
reflected by objects on the surface of the Earth or within the Earth's 
atmosphere The variations m reflectivity, as recorded at the output of 
the ground station, provide clues to the nature of the object The division 
of the intensity of reflectivity into spectral bands affords further infor- 
mation about the scene to be analyzed 

In most imaging systems, the quality criterion is that the imagery 
be pleasing to the observer Those qualities that present a pleasing 
appearance are determined by the characteristics of tonal reproduction, 
sharpness (as contrasted with resolution), and, m color reproduction 
systems, the color balance or fidelity both m hue and purity In general, 
it may be said that a pleasing reproduction of a scene is one that the 
observer feels is a reasonable likeness of the original 

For scientific work, the criterion is more stringent The data 
received from the imaging system must be capable of interpretation with 
precision and repeatability with regard to position within the image field 
(geometry), intensity of the reflected radiance with regard to object 
characteristics, both spectral and intensity, and the characteristics of 
the radiating source, modified by the characteristics of the transmission 
path (atmosphere, radiometry) It is also necessary that the edges of 
objects seen in the image be sharply defined (MTF) 

Error models of each of these aspects of the RBY and the MSS 
system to be flown aboard ERTS have been constructed Fixed systematic 
errors may be compensated for either in the satellite system or during 
processing of the imagery at the data processing stations 
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A second form of systematic errors is one m which the errors are 
variable but predictable, based on Sun angle, satellite orientation, and 
similar criteria If these variables are known, corrections can be 
applied to the imagery or in the interpretation of the imagery to reduce 
the errors to within close tolerances 

A third form of errors is one m which the errors are random, and 
generally cannot be predicted or corrected Knowledge of the nature of 
these errors permits the calculation of possible deviation from prescribed 
tolerances 

The reflective characteristics of objects on the Earth's surface may 
range from completely diffuse to completely specular Generally, the 
characteristics will lie somewhere between these extremes A com- 
pletely diffuse (Lambertian) surface reflects energy as a function of the 
cosine of the angle between the normal to the surface and line of interest 
Since the projected area of the object as viewed along the line of sight 
varies inversely as a function of the cos me, the radiance of the Lamber- 
tian surface is independent of the angle at which it is viewed Diffuse 
surfaces appear equally bright in all directions 

Specular surfaces obey the law of reflection for mirrors All 
energy is reflected at an angle equal to the angle of incidence. Unless 
the line of sight lies along this line, little energy will be received by the 
sensor Energy reflected specularly and viewed by the sensor along the 
angle of reflection will be very intense and may cause overload of the 
sensor or the video transmission channel 

The sensor measures the intensity of the reflected light from the 
scene by exploring an image formed optically on the photosensitive 
surface, point by point The accuracy of the measurement relative to 
the actual intensity of the target depends on the quality and spectral 
characteristics of the illumination, the attenuation by the atmosphere, 
and the flare light caused by the radiance of the atmosphere and m the 
optical system Atmospheric haze varies with locality as well as 
weather, and its attenuation of information bearing light is a variable 
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Figure 5-28 

60° SOLAR Zenith Angle 
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Rayleigh scattering of radiance caused by molecular characteristics 
of air m the atmosphere varies inversely as the fourth power of the wave- 
length Blue light is scattered more than red Larger atmospheric 
particles on the order of 10 times the wavelength of light, such as smoke, 
dust, water, fog, and other particles, also cause a scattering, which is 
not spectrally variable, known as Mie scattering 

On a clear day, the scattering is predominantly Rayleigh and is 
evidenced by the clear blue sky In hazy weather, Mie scattering is 
predominant with a neutral gray overcast or white clouds that reflect 
radiant energy uniformly throughout the spectrum 

Scattered light reduces the modulation of the target, causing a 
reduction m contrast in the r-eproduction Compensation for this reduc- 
tion of contrast may be applied by increasing the gamma of the reproduc- 
tion In a linear system, a value of luminance may be subtracted from 
the video signal (clipping), increasing contrast of objects that are of 
greater radiance than the flare light It is difficult to bring out details 
of objects that are of lesser radiance than the radiance of flare light 
since the noise (random fluctuations) which is always additive, cannot be 
removed 

An analysis utilizing a Monte Carlo computer program, developed 
by TRW illustrates atmospheric effects A large number of (numerical) 
photons are generated, their paths are followed down through the atmos- 
phere, reflected or absorbed at the surface, and the returns counted 
The output which is the fractional photon return, is normalized by the 
solar constant and plotted in these terms 

The first analysis considers a large uniformly reflecting target 
and both hazy and clear atmospheres No distinction is made between 
target/signal photons and atmosphere/noise photons If one chooses a 
given albedo and wavelength, the radiance at orbital altitude is readily 
obtained from Figure 5-28 

The second analysis contains several refinements The first 
analysis considers a large uniform target Actually, any target has 
small structure and here we consider a small target of variable albedo 
surrounded by a larger area of constant (20 percent) albedo Further, 
we may now define a target return which includes the radiance reflected 
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off the target and an atmospheric return which includes everything else 
The target return contains the functional dependence on albedo which 
interests us, the atmospheric return is an additive noise-type term 
which degrades resolution The radiance at orbital altitude is the sum 
of the two It is important to note that the atmospheric return can vary 
greatly and the results obtained apply only to a particular case Fig- 
ure 5-29 shows these results Note that the two analyses are consistent 
only to about 30 percent as certain assumptions differ 


The results of these analyses, particularly the second, indicate 
that atmospheric radiance is a serious problem when one wishes to 
determine spectral signatures Further work is necessary to model 
these effects and determine the accuracy with which they may be removed 


5 6. 2. 2 In-Flight Sensor Calibration 

Absolute m- flight radiometric calibration of the signal transfer 
function of the RBV and the MSS system will depend on the availability of 
an absolute radiometric source of energy, modulated by known step 
values If the step values are measured at the output of the system when 
they are plotted against the absolute input values the transfer function can 
be determined 


It is proposed that the input energy source be obtained from the 
erase lights, normally used m the preparation of the photo conductor 
surface before exposure Several levels of intensity may be used, each 
exposing the RBV surface uniformly and being read before the next value 
is used If the method of prelaunch calibration is used, the input values 
are known, and the signal transfer function can, m principle, be obtained, 
as previously described 

Aging effects of incandescent lamps can be compensated for by 
using a second set of lamps as a primary calibration standard The 
limited use of these lamps will prevent serious error that results from 
again effects over the life of the onboard systems 
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Figure 5-29 

NADIR SPECTRAL Return Radiance 

E a = Atmospheric Return (For Average Terrain Reflectance of 20%) 
E_ = Target Return (For Target Reflectances 5% 20%) 

(50° Solar Zenith Angle) 
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This concept is similar to that used m conventional photographic 
sensitometry The lamp, used as the exposure source m daily work 
with the sensitometer, is compared on a more or less regular schedule 
with a second lamp which is reserved as a primary standard and used 
only for calibration purposes Adjustments in the voltage and current of 
the everyday lamp are made as the result of each calibration to reduce 
the errors caused by aging effects When the lamp fails, the primary 
standard lamp is used to calibrate its replacement so that no discon- 
tmutity m the calibration need exist 

Sources of error m calibration will be caused by the aging of the 
primary standard lamp, variations m lamp voltage (lamp output varia- 
tions in lamp voltage (lamp output varies as V^), and changes m color 
temperature Changes m the transmission of the optics will not be 
corrected or detected, unless the exposure is made through the total 
system 

A more basic approach to the problem of calibration m flight is to 
use the Sun, which, since the system is operating above the atmosphere 
with all its variables, represents an ultimate m stability both m irradiance 
and spectral content An integrating sphere can be mounted at a convenient 
position on the spacecraft Fiber optics conduct the energy to a photocell, 
placed so that the erase lights can be compared to the sun by chopping or 
alternating the light from one source to the other, and adjusting the intensity 
of the Sun- derived energy automatically to produce a null Readout of the 
amount of attenuation required for balance will calibrate the system 

A better system of m- flight radiometric calibration is using the 
sunlight to provide the exposure source However, since this may 
present severe engineering problems it may be impractical 

Basically, the best system of m-flight calibration is one which 
requires the least dependence on preflxght calibration The use of 
several exposure cycles, each illuminating the total surface of the 
vidicon tc a different intensity (as proposed by RCA) is dependent on 
preflight calibration, and will not necessarily detect variations that 
result from aging or other causes A single exposure of several gray- 
scale step tablets, placed at strategic points over the surface of the 
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photo conductor, would produce more reliable information for the dynamic 
signal transfer function than the proposed single intensity step method 

It should be noted that at least a two-level signal is always required 
because, unless the system is capable of transmitting dc levels accurately, 
a knowledge of the output intensity relative to input is indeterminate 
Normally, the two-level signal contains what is termed a black level 
reference This portion of the signal waveform represents zero radiance 
input This knowledge of intensity value is used to restore dc levels at 
the output end of the system 

5. 6. 2. 3 Sensor Differential Shading 

Shading variations that are acceptable m a black and white 
reproduction may not be acceptable m a three-color system Differences 
m the responses of the three RBV tubes as a function of position on the 
photo conductor when exposed to a uniform spectrally neutral field may 
cause noticeable departures from neutrality m the reproduction when 
printed onto color material Variations noticeable as color shading of 
reproduction of a gray field will indicate the need for further shading 
correction beyond the correction that is inserted m each RBV camera 
systems 

5. 6. 2. 4 Gray-Scale Tracking 

In all three-color systems, a certain ratio of energy inputs to the 
three sensors is specified as neutral or gray As the intensity of the 
input energy is increased while the specified ratio is maintained, the 
luminance of the reproduction increases (lighter grays) toward white 
A change m the ratio will produce color 

Adjustments of the three RBV camera systems to maintain 
gray- scale tracking over the range from black to white is critical if 
spectral signatures are to be reproducible m the false color hard copy 
output of the system 

Differences m nonlinearity of each camera and associated video 
amplifier and differences in black level setting affect the tracking 
capability of a three-color vidicon camera system If the system is 
adjusted for correct tracking of large area gray steps, errors in small 
area gray steps at the black end may still exist because of flare light 
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addition Since the photo conductor of the vidicon may exhibit increased 
transparency toward the red end of the spectrum, this flare increases 
for the red as compared with the blue-green Before flight, the RBV 
system should be calibrated for black level shifts m small area blacks as 
a function of the average irradiation reaching each sensor Some correc- 
tion can be made m the image recorder by resetting the black level of 
the recorder for each spectral band 

5. 6 2, 5 Gamma Correction 

The term gamma, originating m photography, denotes the 
maximum slope or gradient of a characteristic (D-log E) curve plotted 
on log-log paper Also, for signal changes of a small percentage, the 
value of gradient or gamma is the ratio of relative output to relative 
input For example, when a gradient of 2 0 is specified for a given point 
m the characteristic curve, it means that a 2 percent change m input 
will cause a 4 percent change m output If the characteristic curve that 
describes the dynamic transfer function of log radiance input to log 
radiance output were a straight line, the gamma would be a constant at 
all points This does not occur m practical systems, so that gamma 
refers only to the maximum slope of the curve The term gradient 
should be applied to describe the derivative of the curve at all points 
other than the maximum 

Thus, the measure of gamma is a measure of the enhancement or 
reduction of contrast of a reproduction relative to the original If the 
gamma is greater than unity, contrast differences m the object are 
increased m the reproduction If gamma is lower than 1, contrast is 
reduced When gamma is unity, the reproduction system is linear at 
the point of maximum gradient 

The characteristic curve of practically all image reproduction 
systems has an S- shaped characteristic, so that the gradient m the low 
light or toe and the high light or shoulder of the curve is considerably 
less than the maximum gradient or gamma that is measured m the 
middle of the curve This lowering of gradient causes a compression 
of gray step values at these points Generally, this is considered 
beneficial since it permits fitting a relatively large input dynamic range 
into a somewhat more restricted dynamic range of the reproducing system 
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Photographic systems -using negative films that are printed onto 
positive materials for final reproduction have a minus gamma for each 
step, negative and positive The system gamma is the product of the 
two minus gammas, which results m a positive gamma This is impor- 
tant m electronic video systems, since simple polarity inversion of the 
video signal is sometimes used to produce a direct positive recording 
with the result that tonal values are greatly distorted in the reproduction 
A nonlinear amplifier whose output is inversely proportional to the input 
is required to correct the video signal for direct positive film recording 

This discussion is based on the plot of the characteristic curve m 
log -log coordinates This is significant when the human observer is part 
of the system because, as described by the Weber-Fechner law, the 
observer’s response to a change m stimulus is proportional to the 
existing stimulus In other words, the eye’s response is logarithmic, 
to a best approximation, so that equal increments along the abscissa of 
the log-log plot represent equal increases m the gray- scale step value 
For example, the range from 0 1 to 1 foot-lambert has as many percep- 
tible steps of gray as the range from 10 to 100 foot-lamberts Also, 
m a dynamic range of 30 1, there are as many gray- scale values between 
3 3 and 18 percent as there are between 18 and 100 percent 

The RBV light signal transfer function has a maximum gradient 
(gamma) of 0 65 This is a function of the photo conductor surface This 
value is complementary to the maximum gradient of cathode -ray display 
tubes, which is about 2 0 The product gamma of a vidicon and the 
cathode-ray tube is somewhat higher than 1 0, which produces good 
tone reproduction 

The laser beam recoi der (L.BR) is a highly nonlinear device, 
because of the modulator characteristics The modulator nonlinearity 
is corrected by means of a circuit whose sighted function is the inverse 
of the LBR modulator This circuit can be adjusted to provide a linear 
transfer function over a dynamic range greater than 100 1, or to provide a 
further change of gamma to compensate for nonlmearities elsewhere m the 
system 
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If no gamma correction is used in the satellite RBV camera 
subsystem and the vidicon has a gamma of 0 65, gamma correction can 
be provided as part of the L.BR system The amount of correction to be 
used will be based on the gamma of the photographic process (film, 
chemical processing, and reproduction) 

The RBV noise is caused by shot noise inherent m the scanning 
beam The noise varies as the square root of the beam current Since 
the beam modulation is negative, i e , maximum return beam occurs 
for darker areas and minimum beam for bright areas, the RBV noise 
becomes greater m the darker areas of the image as compared to the 
brighter areas 

An increase m gamma increases the signal increments m the 
white end of the transfer function and decreases it m the dark end The 
perceptibility of noise for the RBV will be lessened by this approach 

When the video signal is applied to the image recorder, two 
adjustments are required The gam or amplitude of the video signal 
must be set, thereby setting the exposure range for the film record 
Also, a dc voltage or bias must be added to the video to position the 
exposure at the optimum point On the modulator -film characteristic 

In a linear system, these two _ operations are straightforward and 
each operates independently without affecting the other, if the video 
signal reference level is a black level reference and is clamped to the 
inserted bias The dc bias then sets the value of exposure for the 
darkest part of the video, while the gam adjustment sets the maximum 
exposure 

5 6. 2 6 Tone Reproduction of ERTS Imagery 

A multispectral acquisition system provides a means of enhancing 
differences m recorded imagery by breaking the spectrum of an object 
into several component bands This allows increased information 
discrimination Tonal differences that may exist m narrow spectral 
bands are often masked when integrated over a broad spectral range 
For a multispectral band acquisition system, the system should be 
balanced to acquire the maximum amount of information discrimination 
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in each, individual band Tests have indicated that this condition is 
satisfied when the dynamic range of each band is on the order of 10 1 
and when the individual bands are of equal dynamic range 

The ERTS system will utilize tonal change detection to enhance our 
ability to detect and recognize Earth resources information Since 
multispectral change detection is the mam goal, the system should be 
precalibrated for maximum information acquisition for nominal image 
acquisition conditions Once these conditions have been established they 
should be fixed, and at least m bulk mode reproduction should reflect the 
same system parameters In this way, any tonal or color changes m the 
final output imagery will reflect only changes m the appearance of the 
imagery as seen by the satellite and will not reflect changes m either 
the satellite telemetry or reproduction systems In special cases, such 
as m the precision mode, it would be advantageous to deviate from this 
fixed tone reproduction standard 

Because of atmospheric attenuation and scattering, one would 
expect the dynamic range of the multiband records to be different 
Atmospheric scattering, often referred to as haze or atmospheric 
radiance, is basically a wavelength dependent function and varies from 
being completely selective to being almost totally nonselective depending 
on the types and concentrations of particles in the atmosphere For pure 
Rayleigh atmosphere, the scattering goes inversely as the fourth power 
of wavelength For a hazy atmosphere with large concentrations of 
water vapor and smoke, the scattering is almost wavelength dependent 
Some authorities have stated that a typical atmosphere combining both 
Rayleigh and Mie scattering goes approximately as a function inverse to 
the 1 6 power of wavelength In any case, the nominal conditions for 
scattering are predominantly m the ultraviolet and blue region of the 
spectrum and fall off significantly as one goes to the red and infrared 
regions of the spectrum For the ERTS case, one would expect the 
apparent contrast m the red band to be higher than that m the green, and 
that m the infrared to be higher than that in the red 
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5 6 2 7 Tonal Balancing of Records 

Test results at Itek have indicated that optimum results are 
obtained when the density range of the individual spectrum records are 
approximately equal If the density of the records are equivalent, the 
printing and color balancing of the resultant additive color records are 
simplified A nominal dynamic range on the order of about 10 1 
(equivalent to a density range of 1 0) has been found to be ideal for 
making additive color records Most color reproduction materials have 
a fairly high gamma with a relatively narrow exposure latitude Too 
high a contrast on the separation records can result m a loss of informa- 
tion because the color reproduction material is not able to record the 
entire tonal range of the input imagery The density range should also 
be kept down because it makes exposure and color balancing easier in 
the additive color process 

The individual separation records can be dynamically balanced to 
some extent either m the satellite or m the ground handling equipment 
Some correction may be possible m the satellite by adjusting the gam 
control on the individual RBV cameras Thus, it may be possible to 
telemeter equivalent signal ranges for each RBV This will have the 
effect of altering the apparent scene contrast range in one or more of 
the records to produce a near balance This would be done in a cali- 
brated fashion so that the actual apparent object radiance could be 
reconstructed if desired Once this calibration is set, the RBV 
responses would not be changed during the life of the system Uncon- 
trollable system changes that result from aging and temperature 
variations would be determined by calibrated reference standards within 
the satellite 

The input dynamic range can also be altered m the initial printout 
state Here the taped input data for each channel can be computer 
analyzed and a corrected signal can be fed to the LBR In this manner, 
the dynamic range for each record can be either raised or lowered to a 
prescribed level to produce a nominal standard equivalency Again, 
once this condition has been established for the bulk mode, it should 
remain constant 
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For the precision mode, one might want to deliberately* deviate 
from this standard For example, there might be three possible 
atmospheric conditions, heavy haze, moderate haze, and clear Specific 
pre calibrated correction functions could be stored m the computer for 
these three cases The computer would input these specific input speci- 
fications to the LBR to balance the dynamic ranges of the individual 
records for each specific condition 

A third possible way to correct the contrasts of the multispectral 
black and white records is by processing each individually to different 
gammas However, this is not recommended, since the control and 
repeatability of such a process would be a complex process with a high 
potential error The major tone reproduction corrections should be 
made m the computer and LBR stages of the reproduction process 

Itek's experience has shown that, m subsequent duplication stages, 
a near ideal tone reproduction produces best results This occurs at 
a process gamma of 1 0 at which there is 1 1 correspondence between 
log exposure range and density range Once the film type and processing 
conditions have been established to achieve the optimum nominal condi- 
tion, they are held constant to ensure repeatable results m all future 
reproductions 

If a calibrated reference step wedge is generated in the LBR, the 
best process control and color balancing of the imagery is obtained 
This step wedge can be used as a process control standard for subsequent 
production stages By monitoring the densities reproduced on the wedge, 
proper exposure and photographic processing control can be assessed 
The three integrated calibrated step wedges on the RBV records should 
produce a neutral black and white tone Deviation from a neutral tone 
would indicate an improper balancing of one or more of the multi- 
spectral print channels used to produce the additive color record 

Neutral tone balancing has been used as a standard m color 
printing for many years In spectral cases, balancing to a neutral tone 
for certain types of images may not produce optimum results In the 
precision mode, the balance maybe altered to produce ac records that 
meet the optimum needs of a user However, this would be a special 
process 
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5 6 2 8 Radiometric Accuracy of a Photographic Process 

Since the ERTS imagery will extract radiometric data from the 
ground scene, it is necessary to carefully consider the accuracy of 
photographic materials used as quantitative radiation sensors A well 
controlled photographic process can produce remarkably accurate and 
repeatable radiometric and photometric data 

Several types of errors can occur when film is used as a 
radiometer The first error is caused by the nonuniformity of sensi- 
tivity m the typical photographic material Any nonuniformity m 
sensitivity across the film web will produce errors m the developed 
imagery However, the uniformity of the quality films made today is 
exceptionally good, and sensitivity variations are less than 5 percent as 
a result of manufacturing tolerances As long as the film is stored m 
a suitable environment, uniform sensitivity is maintained until the film 
is used The absolute sensitivity of photographic materials can vary 
with each batch However, batch variations can be determined m a 
quality controlled sensitometric operation Batch variations in 
sensitivity might vary as much as 10 to 15 percent m a normal photo- 
graphic material 

The largest radiometric error generally occurs m film processing 
To achieve good radiometric results, the uniformity of the processing 
must be well controlled Each part of the film must receive the same 
chemical treatment This requires uniform agitation and uniform 
temperature as well as constancy of developing solutions Processing 
uniformity on the order of 5 percent can be achieved in precise sensi- 
tometric processing However, m a normal well controlled production 
processing operation, variations might run as high as 10 percent This 
variation would be the total variation across the format and not the point 
to point variation between adjacent detail The relative radiometric 
accuracy will be significantly higher 

The size of the object or objects being evaluated influences 
radiometric accuracy Photographic film is composed of small grains 
of silver These grains represent the noise m the photographic process 
When a large area is being evaluated, we are well above the noise level 
of the system and the radiometric accuracy will be very good However, 
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as we begin to evaluate smaller objects, we begin to get into the gram 
noise of the system and density fluctuations as a function of this gram 
noise begin to appear We begin to lose accuracy as the signal to noise 
level increases The film that has been selected for ERTS duplication, 
EK-2430, has the lowest granularity (18) of currently available produc- 
tion duplicating materials and it will have a minimal effect on the 
radiometric accuracy of the system 

Exposure effects could also affect photographic materials used 
as quantitative radiation sensors The uniformity of the exposure on 
the film will affect results Most quality printing devices can be made 
to accuracies of within a few percent 

The possible reciprocity effect that results from large exposure 
variations could cause error However, we do not anticipate this prob- 
lem m the ERTS system because it will be basically fixed and extra- 
ordinarily large ranges are not anticipated 

Latent image decay is another potential source of radiometric 
error There will be an essentially constant time m all photographic 
operations so that significant time variations between exposure and 
processing are not anticipated 

After considering these factors, it would appear that a well~ 
controlled photographic process involving good sensitometnc calibration 
and control procedures could maintain a radometnc accuracy on the order 
of 95 percent for objects whose images are m the size range of a milli- 
meter and larger The objects whose images approach the grain noise 
limit of the film will have larger errors 
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APPENDIX A 

COORDINATE SYSTEMS AND TRANSFORMATIONS 

In order to facilitate the discussion of Sections 4, 2 and 4 3 the 
right handed orthogonal coordinate systems used m those two sections 
are defined in this appendix Also the required transformations among 
the various coordinate systems are given. 

1 COORDINATE SYSTEMS 

1 1 Earth Centered Inertia (ECI) 

T 

The earth centered inertia system, x.j = (x^y^Zj) , has it origin 
at the center of the earth, x and y axes lying m the equatorial plan with 
the positive x axis directed toward vernal equinox 

1 2 Earth Centered Fixed (ECF) 

T 

The earth centered fixed system, x:^ = (x^,, y^, z^) , has its 

origin at the center of the earth, x and y axes lying m the equatorial 
plane with the positive x axis directed toward the prime meridian 

1. 3 Orbital Reference (OR) 

T 

The orbital reference system, x^ = Yr* z r) 1S centered at 
the spacecraft position m orbit and oriented such that the positive z axis 
points down to the center of the earth. The positive x axis lies m the 
orbit plane and has a negative dot product with the spacecraft velocity 
vector, and the y axis completes the right hand set. 

1 4 True Body (TB) 

T 

The true body coordinates, = (x^y-g, , are fixed m the 

spacecraft body and are defined with respect to the horizon scanner 
planes where it is assumed that the horizon scanner planes are orthogonal 
to each other The positive z axis points downward along the intersection 
of the A, C, and B, D scanner planes The position x and y directions axe 
as described m Figure A-l, 

1 5 Nominal Orientation (NO) 

T 

The nominal orientation coordinate system, 2 £q- ( x q) Yq> 2 q) * 1S 

centered at the spacecraft position in orbit and oriented such that with 
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ORBITAL VELOCITY HORIZON SCANNER PLANES 



Figure A-l Body Coordinates Defined Based 
on the Horizon Scanner 

the horizon sensors tracking the horizon profile p^ = p^. and Pg = Pp 
where p is the angle the ith sensor makes with the +z axis With the 
positive direction of the z axis so defined, the positive x axis lies m the 
orbital plane and has a negative dot product with the spacecraft velocity 
vector The y axis completes the right hand set 

2 MAGNETIC COORDINATE SYSTEM (MC) 

T 

The magnetic coordinate system x^ = (x^., y^, z^) has its 
positive x axis at longitude -150° and latitude -11° with respect to the 
EFC system The y axis lies m the equatorial plane and z completes the 
right hand coordinate system 

2 1 T ransformation 

The transformation from the ECF system to the ECI system is 

—I = T l— F (A-l) 
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where 


COS OOjt 


-sm u 0 


T i = 


sm cj^t 


COS LO^t 


(A-2) 


where 

= earth's rotation rate 

t = time of passage of the prime meridian past vernal equinox 
The transformation from OR system to the JECI system is 


-R = T 2 — I + — 


(A-3) 


where 


cosQ smv 
■fsinfi cos t cos v 


-cosft COS i COS V 


+ smfi smu 


-sm t cos v 


T - 
X 2 


smS7 sin <. 


-cos Q sm t 


cos l 


(A -4) 


-cos£2 cosv 
+ sm£2 cos i smv 


-smfl cost; 


-cosfi cost smv 


-sin t smv 


and where for the spacecraft orbit (see Figure A-2) 

£2 - longitude of ascending node 
t = inclination of orbit 

v = spacecraft angle orbit from the equatorial plane 
T T 

R = (0, 0,R) , R - distance of spacecraft from the center of the 
earth 
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Figure A-2 Relationship Between ECI and OR Coordinate System 
The transformation from the NO system to the TB system is 

*B = S^ 0 (A-5) 



cos 0 sin 41 


sin sin 8 sin 

+ COS i> COS 4* 


cos sin 0 sin^ 
-sin <j> cos 4* 


-sm 0 


sin <j> cos 0 (A-6) 


cos <f> cos 0 


Here the Euler angles are yaw ( 41 ) about the z axis, then pitch (0) about 
the y axis and finally roll <b about the x axis as shown m Figure A-3 
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TRUE BODY COORDINATES Xg, Y fil Zg 

NOMINAL ORIENTATION 
COORDINATES Y_ Z 

O O O 

EULER ANGLES YAW $ (ABOUT Z AXIS) 

PITCH 0(ABOUT Y AXIS) 
ROLL 0(ABOUTXAXIS) 


EULER ANGLES MEASURED FROM NOMINAL 
ORIENTATION TO TRUE BODY COORDINATES 



Figure A-3 Euler Angles Measured from Nominal Orienta- 
tion to True Body Coordinates 


The transformation 


from the OR system to the NO system 


— O " T 3-R 


(A -7) 


will now be derived From (4-11) we have 


tan (3 = 



AC±B 


A 


(A-8) 


which is used to compute for the two roots of (3 , |3 ^ and |3 ^ For a = 0° we 
compute from (A-8), a ^ 


6 


R 



(A-9) 
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For a = 90° we compute from (A -8) a 5 niri , ru = 6-q 

sr 1 JL Orl ir 


6 P " 2 (,3 2 " 


Given & R and 6p and the definition of the NO and OR system, the 

direction cosine matrix relating these coordinates system can now be 

derived Let x , v , z be the unit vectors m the NO system The com- 
o o “o p 1 

ponents of z - (x _,y „ , z _) m the OR system is solved by the 

O Zj JA. 2 y K Z y 

following relationships (see Figure A-4) 


*z,R = sm5 P 


^Z.R = 31n 6 R 


(A-ll) 


*v p = 1 ~ p - y 


2 _ 2 
z, R 7*. R 


We now have 


-o ~ (x z, R, y z, R, Z z, R* 


(A-12) 


The components of x q = (x^ R , R , z y , R ) and of ^ = (x y ^ R , y y } R , z y > R ) ‘ 
are obtained by solving the following dot product equations respectively 


z * x = y.R 
— o — o J z’ 


z R 
z’ 


1 - X 


(A -13) 


Figure A-4 Z Q Axis of the NO System 
and its Relationship to the 
OR System 
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and 



_X z,R 
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z V = 

o •*’o 

^z,R 
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__ Z z,R_ 


, 2 
L 1 - y y.Rj 
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The equations are quadratic m the unknown and y^ r ^ The positive 

roots are the required solutions The transformation is therefore 



— 

X x,R 

0 

1 - tj 

x, R 

T 3 = 

0 

y y,R 

, 2 
1 - y y,R 


_ X z,R 

y z ,R 

Z, R 


(A -15) 


In the development of the system equations we need to relate the 
body rates w ^ z a bout the body x,y, z axes respectively to the 

Euler rates ($,0 ,i|j) with respect to the rotating NO system rotating at^ Q 
To do this we successively apply the rotation matrix Sj(4>,0,4») given by 
(A-6) to each oi the Eules rates as follows 
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where 


s 2 (<M) = 


1 tane sinij) tan e cos 4 >" 


cos 4 > 


-sin 


sin<j>/cos 6 cos<j>/cosG 


(A-17) 


If we denote the transformation S^, Equation (A- 6 ) by 


s x (<|»,e,»|i) = s 1 


(A-18) 


the transformation from the MC system to the EFC system is given by T 


x_ = T .x. . 
— F 4“ M 


(A-19) 


where 


T 4 = S 1 ( 0 °, 11°, -150°) 


(A -20) 


Note that. 


S 3 S 2 S 1 
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and the derivatives of and are given m terms of Euler angle rates 
as follows 


S 0 =1 
2 


2 2 

0 4 >{tan 0 cos <j> ) + 0 { sec 6 sm 'f’) 0 (sec 0 cos<}>) -<|)(tan 0 sm<f>) 


- 4 > sin<(> 
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and 
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APPENDIX B 

ERTS ERROR ANALYSIS FIGURES OF MERIT 


The absolute and relative error covariance matrices C c and C At 

6 A S 

completely describe the error probability density functions for normally 
distributed error sources 


p e x .i y ) = 


2tr| C|| 1/2 


exp 




{B- 1) 


y > ■ 2,r|C 1/2 e *P 


(B-2) 


Several "figures of merit" can be used to describe ERTS performance in 
terms of the covariance matrices and C AZ The most common ones 
are 


• One and three-sigma position error ellipses 

• CEP 

• Variance m the distance between two points on an ERTS 
picture 

The first two figures of merit can be applied to both absolute and relative 
error (internal consistency), while the distance error variance is only 
applied to the internal consistency 

ERROR ELLIPSE 


The 1 -sigma error ellipse is obtained from the equations 



By substituting equations (B-3) and (B-4) m equations (B-l) and (B-2) it is 
obvious that the 1 -sigma error ellipse is the locus of points with equal 
probability densities 


P <e x . iy) 
P (A 5x- 



(B-5) 
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The integrated probability density inside the 1 -sigma error ellipse 
A is (m principal axis coordinates ^ and 


'lo- JJ d ^i d ^2 z^T z exp { 
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_n__ b 2 
2 + 2 
a l ^2 
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Equation (B-7) can be transformed to the integral over the unit circle 


P l<r = eXp {'I [ w 1 2+w 2 Z ]} 


(B-8) 


where 


^2 

w = — w_ = 

1 1 2 *2 


or, m polar coordinates the 1-sigma error ellipse probability is 


P 
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1 2tt 

= // £ exp [4 p2 l dpde 


0 393468 


(B-9) 
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The 3 -sigma error ellipse probability is, similarly 


P 


3(r 


3 -Zv 



0 988890 (B-10) 


The numerical values (B~9) and (B-10) give the probability of finding an 
ERTS picture point within the 1-sigma and 3-sigma error ellipses 
respectively 

DISTANCE BETWEEN TWO POINTS 

Define two image points as follows 

Point 1 (x^ + AXj, y^ + Ay^) (B-ll) 

Point 2 (x^ + ^ x 2 > y 2 + Ay 2 ) 

the square of the distance between the points is given by 

(S+AS) 2 = ^(x^+AXj) - (x-,+Ax 2 )j * [^1 +A ^1^ _ ^ y 2 +Ay 2^] 2 (B-12) 

or, to first order m AS 

S 2 + 2SAS = (x^x^ 2 + (y i-Yz) 2 + 2 ^ x l“ x 2^ A ^x + 2 ^ y l“ y 2* A ^y 

where 

A§ - Ax. - Ax. (B-13) 

-X 1 Ar 

A£ y = A Yi - Ay 2 
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thus 


(x -x ) (y,-y ? ) 

AS = — A§ + — — A$ 

s x S b y 
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(B-14) 


where 




m 


A^ = col (A£ ,AL ) 
y 


the variance m the distance 



is then 


°AS 2 " <as2 > ~ is C a£ -S 


(B- 15) 


Equation (B-15) is a quadratic form (positive definite) m terms of the 
distance unit vector d^g It is obvious from (B-15) that the distance error 


variance lies between the maximum and minimum eigenvalues of C 




°min - 17 AS “ °~max 


(B-16) 


CEP 

The CEP or circular probable error is defined as the radius of the 
circle around the correct image point for which there is a 50% probability 
of the erroneous point being within Thus, it is the value of R for which 
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An approximate solution to (B-17) which is accurate to better than 3% is 


CEP = 0 589 o* + <r. ] 

l £. i 


(B- 18) 


This equation is accurate for the cases of interest m the ERTS error 
analysis 
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APPENDIX C 

DERIVATION OF THE NOMINAL ROLL, PITCH, 
YAW DIRECTION COSINES 


Let the roll, pitch and yaw axes be denoted by e X) ey, e z unit 
vectors respectively 


e = 

X 1 + X 1 4- X k 

— X 

x— y-* z— 

e = 

Yi+Yi+Yk 

-y 

x— y* z— 

e " 

Zi-fZ j + Z k 

— 2 

x— z— 


where 


x, y, z = earth centered inertial (ECI) coordinates with 
x pointed m the direction of the vernal equinox 

_i , 2 ., k = unit vectors along the ECI axes 


The e^. (yaw) axis is directed downward along the local geodetic vertical 
of the ellipsoid representing the Earth. The &x (roll) axis is oriented so 
that the orbital rate signal is sensed by the roll-yaw gyro Thus, ex lies 
along the e^. vector with positive sense m the direction of flight Hence 


e 

— x 
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with 


B 


the e 

-y 


(pitch) axis is then 



e - e x e 
— y — z — x 


or, m terms of direction cosines 
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The evaluation of the nine direction cosines of the attitude coordinate axes 
proceeds from the determination of the direction cosines of the local 
geodetic vertical The Earth’s figure is defined by 


S = 





(C-6) 


where x e , y e , and z e are points on the Earth’s surface m ECI coordinates, 
a is the equatorial radius, and b is the polar radius The attitude vector 
is 


h = r 


- r 

— s 


= h e 
— s 


(C- 7) 
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where 


r_ = xi. + y _^ + z k, satellite position 

r = xji+yj+zk, sub- satellite position 
s s s — s 

h = scalar altitude 
_ c 

— s = |vSj J um * : vec t° r normal to Earth's surface 
the gradient of S at the sub- satellite point is 



If the scalar parameter, K, is denoted by 
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K = 


2h 

|VS| 


(C-9) 


then the satellite position and sub-satellite position are related as follows 
{from equations (C-7) - (C-9) 

x = (1+K) 

y = (1+K) y (C- 10) 

s 

z = (l+a 2 /b 2 K) z 

s 
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And since r s satisfies equation (C-6). 


S = 


x + y 

(1+K) 2 


2 „Z. 2._.2 
a (b /a +K) 


- a = 0 


(C-ll) 


the altitude vector h is then 


h - ^ x i -j- k y -i -). k q v 

- " l+K - 1+K — 1+K - 
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where 


q = 


1+K 


1 b 2 /a 2 +K/ 


The scalar altitude is 


, _ K . .2 2 . 2 . 2 

h — ^ — x + y + q 


(C- 13) 


The scalar K is computed by iterating as follows The initial value of K 
is that for a circular earth b = a 
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Then 
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n+1 


= K + AK 


n 
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and, from equation (C-ll) 
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so that 


AK = 
n 



(0-16) 


where 
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x + y 

(1+K n )2 
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a (b / a +K ) 
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dK 
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2(x -f v ) 

(1+K )3 
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The direction cosines of the roll, pitch, and yaw unit vectors can 
now be evaluated The yaw axis e_ z is along the normal to the Earth’s 
surface Hence, 
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APPENDIX D 


CALCULATION OF SIGNAL TO NOISE RATIOS 

Using Hughes data* on the radiance of agricultural targets (with 55° 
solar zenith angle) we can calculate similar signal to noise ratios Taking 
healthy oats and barley as an example, the radiance levels are 


Crop 

Band 1 

Band 2 

Band 3 

Oats 

0 402 

0 275 

0 350 

Barley 

0 396 

0 407 

0 447 


-2 - 1 

Crop Radiance (mW cm ster ) 


Crop 

Band 1 

Band 2 

Band 3 

Oats / Barley 

0 6 

-4 dB 

10 

20dB 

4 

12dB 


MSS S/N Ratio -Extended Crop Target 



Band 1 

Band 2 

Band 3 

S/N Ratio 

-15 dB 
- 1-8 dB 

9 dB 
6 dB 

0 dB 
-3 dB 


Cross Track 
Along Track 


MSS S/N Ratio-200 1 Rows Crop Target 


^"System Design Study Multi-Spectral Scanner (MSS), 11 

Hughes Aircraft Company, Santa Barbara Research Center, 26 August 1969 
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1 SCOPE 

This specification establishes the performance, design, development and 
test requirements for the Earth Resources Technology Satellite Program 

2 APPLICABLE DOCUMENTS 

The following documents of the exact issue shown, form a part of this 
specification to the extent specified herein In the case of TRW 
Systems documents, the latest issue shall apply In case of conflict 

between documents referenced herein and the requirements of Sections 
3, 4, and 5, the requirements of Sections 3, 4, and 5 shall apply 


SPECIFICATIONS 


Military 

MIL-D-1000 
01 March 1965 

TRW Systems Group 

D-13353 

D-13354 

D- 13500 

D-13503 

D-13505 


PR12-1 
D- 13700 
D-13701 
D-13750 

D-13800 


Drawi ngs-Engineermg and Associated 
Lists 


Design Qualification Test, ERTS 
Environmental Acceptance Test, ERTS 

System Specification for Earth 
Resources Technology Satellites A & B 

Thor/Delta Launch Vehicle Interface 

ERTS Specification Government Furnished 
Parts List and Ground Support 
Equipment List 

Identification and Marking 

ERTS Ground Data Handling System 

ERTS Personnel Subsystem 

ERTS-Operations Control Center 
Subsystem 

ERTS NASA Data Processing Facility 
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STANDARDS 


Mil itary 

MIL-STD-100 
01 March 1965 

MIL-STD-129D 
11 April 1969 

MIL- STD-1 43A 
14 May 1963 

MIL-STD-803A-1 
27 Jan 1964 


Engineering Drawing Practices 
Marking for Shipment and Storage 


Specifications and Standards, Order 
of Precedence for the Selection of 

Human Engineering Criteria for 
Aircraft, Missile, and Space Systems 
Ground Support Equipment 


MIL-STD-1247A Marking, Functions, and Hazard 

20 Dec 1968 Designations of Hose, Pipe, and Tube 

Lines, for Aircraft, Missile and 
Space Systems 


MS-33586A Metal, Definition of Dissimilar 

16 Dec 1958 


OTHER PUBLICATIONS 

TRW Systems Group 


PAR 700-53 
11 November 1969 


PAR 700-54 
27 November 1969 

PAR 700-55 
13 February 1970 


ERTS Configuration Management Plan 

ERTS Logistics Plan 

ERTS Maintainability Plan 

ERTS Quality Program Plan 

ERTS Reliability Program Plan 

ERTS Test Monitoring and Control Plan 

Procurement Product Assurance 
Requirements, Quality Systems 
Provisions, Project ERTS (Subcontracts) 

Subcontractor Reliability Requirements 
for ERTS 

Procurement, Performance Assurance 
Requirements, Quality, Reliability and 
Maintainability Provisions, ERTS 
Subcontractors (GDHS Equipment) 
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PAR- 700-56 
27 March 1970 

PAR-700-61 
27 March 1970 

NASA 

GMI 8040 1 
NHB 5300 4 (3A) 

NHB 5300 4 (IB) 
April 1969 


NPC 250-1 
July 1963 

GSFC PPL-10 
July 1968 

S-320-G1 
October 1969 


Configuration Management Provisions 
for ERTS Major Subcontractors 

Procurement Performance Assurance 
Requirements, ERTS Subcontractors for 
GDHS Software 


Configuration Management Instruction 

Requirements for Soldered Electrical 
Connections 

Quality Program Provisions 

for Aeronautical and Space System 

Contractors 

Reliability Program Provisions for 
Spacecraft Contractors 

Preferred Parts List 


General Acceptance Test Specification 
for Spacecraft and Components 
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3 REQUIREMENTS 
3 1 Performance 

311 Performance Characteristics 

The Earth Resources Technology Satellite (ERTS) program shall consist of 
satellites, a Ground Data Handling System (GDHS) and a Data Collection 
System (DCS) The GDHS shall consist of an Operations Control Center (OCC) 
and the NASA Data Processing Facility (NDPF) The DCS shall consist of 
ground based collection platforms and transmitters for transmission of 
data to the satellite The DCS also includes receiver equipment located 
in the satellite 

The satellites shall be launched into orbit by the Thor-Delta launch 
vehicle The program shall also include the ground support equipment 
necessary to support satellite functions The satellite shall consist of 
spacecraft structure (including launch vehicle adapter and payload 
mounts), communication and data handling (including antennas), electrical 
integration, active thermal control, orbit adjustment, stabilization and 
control, and solar array power supply subsystems and payload 
provisions 

The program shall provide mission data necessary to satisfy the requirements 
of the Department of Agriculture, Department of Interior, Department of 
Navy, Department of Commerce and other users concerned with agriculture, 
forestry, geography, oceanography, hydrology, geology, and weather 

3111 Missions 

The ERTS Program overall mission is to gather information about the 
natural and cultural resources of the earth and phenomena destroying or 
threatening these resources Specific missions shall include the 
following 

31111 Agricultural Applications 

The ERTS program shall observe and acquire data for major agricultural 
applications which include 

a) Species identification and measurement 

b) Plant growth rate 
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31111 (Continued) 

c) Factors relating to stress on crops and forests 

d) Assessment of crop vigor and health leading to yield predictions 

3 1112 Geological, Geographical, and Hydrological Applications 

The ERTS Program shall observe and acquire data for major geological, 
geographic and hydrological applications which include 

a) Improved classification of areas by geological or geomorphological 
characteristics such as surface composition, water runoff 
patterns, etc 

b) More accurate monitoring of time- van ant phenomena, such as 
population movements, transportai on flow, and environmental 
hazards to man (air pollution) 

c) Measurements of specific hydrological parameters such as soil 
moisture, snow extent and depth, etc 

31113 Oceanographic and Hydrographic Applications 

The ERTS program shall observe and acquire data for major oceanographic 
and hydrographic applications which include 

a) Measurement of sea state 

b) Location and tracking of major ocean currents 

c) Mapping of sea-ice 

d) Detection of specific phenomena of limited anal extent and 
varying locations, such as fish schools, oil slicks. Red Tide, 
etc 

e) Shoreline analyses 
3112 Mission Support 

To support the mission requirements the ERTS Program shall provide the 
following mission support functions 

31121 Data Handling 

The ERTS Program shall be controlled by the Ground Data Handling System 
(GDHS) The GDHS shall convert the spaceborne sensor and relayed Data 
Collection System signals, in accordance with user requests, into 
interpretable data in a form that provides maximum utility to the system 
user 
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31122 Payload 

The payload of ERTS-A, and possibly ERTS-B shall consist of 

a) Return Beam Vidicon camera 

b) Multi -Spectral Scanner 

c) Data Collection Receiver 

d) Wide Band Video Tape Recorder 

31123 Payload Weight and Volume 

The satellite shall contain provisions for a payload of at least 450 lbs 
in a minimum volume of 13 cubit feet 

31124 Launch 

The satellite shall be launched from the Western Test Range by a NASA 
Thor-Delta launch vehicle Launch and range support, facilities and 
services shall be provided during prelaunch and launch activities 

3 112 5 Orbit 

The satellite shall be placed into a 490 n mi circular, 99 degree 
inclination, sun synchronous orbit 

31126 Global Coverage 

The measurements shall be global m coverage in less than three weeks, 
repetitive in observations at the same local time, with the resultant 
images having minimum distortion 

31127 Image Overlap 

The image overlap and sidelap shall be such that continuous photographic 
coverage be produced in the presence of Satellite Attitude Excursions 

31128 On-Orbit Tracking 

On-orbit tracking shall be provided by STADAN and Selected Ground 
Tracking Stations. 

3 112 9 Command Control and Data Acquisition 

The command and control, and data acquisition shall be performed by 
NASA STADAN/MSFN stations. The satellite shall be capable of transmitting 
PCM data and receiving commands from the ground stations for any 
satellite attitude 
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3 1 1 2 10 User Requirements 

The data products required by the users for the various missions shall 
be as listed in Table I The data shall be transformed into user oriented 
media and stored in retrievable and reproducable form 

3 1 1 2 11 Launch Vehicle Weight Capability 

The launch vehicle shall be capable of launching a satellite weighing 
1520 lbs into the specified orbit 

3 1 1 2 12 Satellite Housekeeping Telemetry 

The program shall be capable of acquiring satellite housekeeping telemetry 
data 

3 1 1 2 13 Satellite Height 

The total satellite weight shall not exceed 1520 lbs (1480 lbs spacecraft 
and 40 lbs interface hardware) 

3 1 1 2 14 Satellite Orientation 

The satellite shall be 3 axis stabilized with one geometric axis 
coincident with the local vertical within 0 7 degree and the rate error 
not to exceed 0 04 degree per second 

3 1 1 2 15 Command and Telemetry Software 

The software shall contain a central file of active user sensor coverage 
requests, merge orbit data with this file to determine the predicted 
opportunities for observation by the sensors and schedule the payload 
events The software shall generate valid observatory commands based on 
the selected events and manual events These commands shall be checked 
for validity, safety and compatibility with the observatory 

3 1 1 2 16 Housekeeping Data Storage 

The GDHS software shall record all PCM health telemetry for storage in 
the NDPF, and shall be capable of generating trend data for off-line 
analysis 

3 1 1 2 17 Command History 

The GDHS software shall maintain a current history of all executed 
commands 
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TABLE I 


ERTS PRODUCTS 


Product Type 

Product Form 


Prints 

Negatives 

IA 

o> 

> 

r— 

-M 

r— 

tfi 

O 

D_ 

Digital Tape 

Listings 

Other 


RBV Bulk Imagery 

L 

r 


r 

r 

r 


Black & White 

IB 

B 

B 


B 

IB 

Color Composite 

v/ 

■ 

■ 


B 

B 

MSS Bulk Imagery 


■ 

■ 


■ 

■ 

Black & White 

B 

e 

B 




Color Composite 



B 


■ 


RBV Precision Imagery 

■ 

■ 

B 


B 


Black & White 

IB 

fl 

R 


B 


Color Composite 


■ 



B 


MSS Precision Imagery 


B 

H 

BB 


Color Composite 

B 

■ 


mm 


Computer Readable Imagery 




i 


RBV 




m 

B 

B 

MSS 




B 


J 

Index/Abstract Data 

fl 



■ 

B 


Catalog 

■ 

j 


>/ 



DCS Data 

B 




B 


Catalog 

□ 



s 

a 


Montage 

■1 



■1 

B 


Sheets 

a 



■1 

fl 


Catalog 

n/ 

n 

j 

□ 

□ 

□ 

Spacecraft Performance 


B 

■ 


B 

B 

Master Digital Data 



■ 


El 

B 

SPECIAL USER REQUESTS 


n/ 




s/ 


Tapes, Listings, 
Plots as Required 

IMAGERY, DATA, TAPES 
LISTINGS, PLOTS 
AS REQUESTED 
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312 System Definition 

3121 System Engineering Documentation 

The ERTS Program top flow diagram is shown in Figure 1 and the functional 
schematic diagram is shown in Figure 2 

3122 System List 

The ERTS Program shall consist of the following systems 

a) Ground Data Handling System 

b) Earth Resources Technology Satellite System 

c) Data Collection System 

d) Launch Vehicle System 

e) Remote Tracking Station System 

f) Users 

3123 System Specification Tree 

The system specification tree is shown in Figure 3 

3124 Government Furnished Equipment (GFE) 

The system shall contain the following Government Furnished Equipment 
Nomenclature Part Number 

1 ) Intercom Subsystem 

2) Communications Panel Assembly 

3) Dial Exchange Phone System 

4) Teletype Terminal 

5) Teletype Punch 

6) MSS Status Monitor 

7) MSS Demultiplexer 

8) Console Power Supply 

9) Equipment Rack 

10) Matrix Monitor 

11) Payload Status Panel 

12) Spacecraft Status Panel 

13) RBV Tape Reproducer 

14) MSS Tape Reproducer 

15) RBV Tape Recorder 

16) Multi -Spectral Scanner 

17) Return Beam Vidicon Camera 7BD 


TBD 
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18) 

Yaw Gyro Assembly 

201371-1 

19) 

Gyro Electronic Assembly 

202163-1 

20) 

Horizon Scanner AC 

A 0G0 

21) 

Horizon Scanner BD 

A 0G0 

22) 

Sun Sensor 

292384-1 

23) 

Sun Sensor 

202744-1 

24) 

Sensor Electronic and Logic 
Assembly 

200932-2 

25) 

Attitude Control Assembly 

200833-1 

26) 

Drive Electronics Assembly 

200934-1 

27) 

Rate Gyro Assembly 

200935-2 

28) 

Inverter assembly 

C 207450-2 

29) 

Control Switch Assembly 

228576-3 

30) 

Array Drive Mechanism 

202540-2 

31) 

2 Reaction Wheel 

C218461 -2 

32) 

Array Shaft Transducer 

212075-1 

33) 

Pneumatic Tank 

C- 111 928-1 

34) 

Valve 

243553-1 

35) 

Pressure Switch 

PT2-3003 

36) 

Regulator 

PT2-3004 

37) 

High-Low Pressure Transducer 

PT2-3000/C1 

38) 

Battery Pack 

323820-2 

39) 

Power Control Unit 

232215-2 

40) 

Charge Control Assembly 

204553-1 

41) 

Charge Control Assemlby 

204553-2 

42) 

Converter Number 2 

206566-1 

43) 

Converter Number 5 

206569-1 

44) 

Converter Number 6 

206569-2 

45) 

Converter Number 7 

206805-1 

46) 

Converter Number 8 

206805-2 

47) 

Converter Number 9 

206570-1 

48) 

Command Receiver 

217907-1 

49) 

137-154 MHz Diplexer Coupler 

218264-2 

50) 

Digital Decoder 

218303-1 

51) 

137-154 MHz Antenna 

217530-1 

52) 

Special Purpose Telemetry Unit 

202650-2 
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53) Ten Watt Power Monitor 

54) Analog Data Handling Assembly 

55) Digital Data Handling Assembly 

56) Low Frequency Timing Assemlby 

57) Signal Conditioner 

58) Digital Recorder Electronics 

59) Digital Recorder Transport 

60) Command Distribution Unit 

61) RBV-BTE Target Collimators 

62) RBV-BTE Spacecraft Command 
Simulator 

63) RBV-BTE Spacecraft Clock Simulator 

64) RBV-BTE Power Supply System 

65) RBV-BTE 

66) RBV-BTE Quick Look Display System 

67) RBV-BTE High Resolution Hard Copy 
Recorder 

68) RBV-BTE Vidicon Thermal Controller 
Generator 

69) RBV-BTE Vidicon Thermal Controller 

70) RBV-BTE Linearity Checker 

71) RBV-BTE Video Signal and Sync 
Processor 

72) VTR-BTE Spacecraft Command 
Simul ator 

73) VTR-BTE Spacecraft Clock Simulator 

74) VTR-BTE Power Supply System 

75) VTR-BTE Telemetry Monitoring and 
Recording System 

76) VTR-BTE Video Signal Simulator 

77) VTR-BTE Downlink Simulator 

78) VTR-BTE Signal Analyzer/ 

Morn tor 

79) VTR-BTE Standard Test Equipment 

80) VTR-BTE Test Cables and Breakout 
Boxes 


214252-2 

201140-1 

218486-1 

218487-1 

201738-1 

207463-2 

207464-5 

201529-3 

TBD 


n!/' 

TBD 
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81) MSS-BTE Precision Collimator TBD 

82) MSS-BTE Spacecraft Command /s 

Simul ator 

83) MSS-BTE Power Supply System 

84) MSS-BTE Telemetry and Test 
Point Monitor 

85) MSS-BTE Precision Oscilloscope 

86) MSS-BTE Video Processor and 
Image Recorder (Single Channel 
Quick Look Capability Type) 

87) Standard Test Equipment 

88) Test Cable and Breakout Boxes 

89) Spacecraft Structure (Kit) 

90) Solar Array Drive Shaft 

91) Louver Banks (+X, -X) (Kit) 

92) Insulation (set) 

93) Horizon Scanner Bracket 

94) 137-154 MHz Antenna Support 

95) Release System (Set) 

96) Interstage and Separation V- * / 

(Kit) TBD 

31241 Ground Support Equipment 

Government Furnished Equipment and Government Furnished Parts are 
as specified in TRW Systems Specification No D-13505 
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Figure 1 ERTS Program Top Flow Diagram 
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Figure 2 System Functional Schematic 
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313 Operability 
3 13 1 Rel i abi 1 1 ty 

As a design goal, the satellite shall have a minimum life time of one 
year in stabilized orbital flight The Operations Control Center shall 
be capable of fully supporting the satellites, following launch, on a 
24 hour per day basis The NASA Data Processing Facility shall provide 
availability, in accordance with TRW Systems Document PAR 700-55, for 
processing observatory data for coverage Cases A or B 

3132 Maintainability 

The ERTS Program elements shall be designed for maximum consideration of 
interchangeability, accessibility, and repl aceability concepts The 
GDHS maintainability shall be in accordance with the contractor developed 
GDHS Maintainability Plan 

3 13 3 Useful Life 

The ERTS-A and B shall be designed for a minimum lifetime in orbit 
of one year The GDHS shall have a useful life of 10 years with 
routine maintenance 

3134 Environment 

The satellite and ground support equipment shall be designed to with- 
stand, without degradation of performance, the following ground, 
transporation and handling, storage, launch and ascent, and orbital 
environments The GDHS shall be capable of operating within a NASA 
furnished building 
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3 13 4 1 Ground Environments 
a) Test Facility 

Temperature 25°C + 3°C 

Relative Humidity 55% or less 

Cleanness Normal laboratory air conditioning 


b) Launch Facility In accordance with the requirements of 

Vandenberg Western Test 
Range 


3 1 3 4 2 Transportation and Handling Environments 

When packaged or otherwise prepared for shipment, the satellite and ground 
equipment shall withstand air, land, or sea transportation methods 

3 1 3.4 3 Storage Environment 

When packaged or otherwise prepared for storage, the satellite and ground 
equipment shall withstand the environments of paragraph 31 3 4.1 

31345 Launch and Ascent Environmen t 

a) Temperature - On stand 

The spacecraft shall be exposed to controlled 
temperature humidity and dust environment in an 
air-conditioned room The available environmental 
control for the spacecraft in the payload fairing 
shall be 


1) Relative Humidity 

2) Conditioned Air 

Flow Rate 
Temperature 


10% to 50% 

50 lbs/m n 
10 to 15°C 
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b) Vibration 

Frequency 
10 - 17 Hz 
17 - 23 Hz 
23 - 200 Hz 


5 - 14 Hz 
14 - 200 Hz 


Thrust Direction 

Amplitude 
(Zero to Peak) 

+ 2 3 g 

+ 4 0 g 

+ 1 5 g 


Sweep Octaves /Min 
4 0 

3 0 

4 0 


Lateral (Two axes orthogonal 
with thrust axis) 


+ 1 3 g 4 0 

+ 1 0 g 4 0 


c) Acoustic Levels 


The estimated acoustic spectra, inside and outside the 
aerodynamic shroud during flight are given in Table II 

d) -Loads.. 

The following limit load factors are the maximum expected 
flight load factors expressed in gravity units (g) 

1) Longitudinal (parallel to the thrust axis) aft 
minus 10 2g combined vectorial ly with 2 Og in either 
direction along any lateral axis (perpendicular to the 
thrust) The minus 10 2g aft load factor is composed 
of 6 2g from the vehicle "steady" thrust acceleration, 
plus 4 08g (vector) vibration occurring in the vehicle 
first longitudinal mode at 17 to 23 cps The 4 Og 
(vector) vibration is at the front end of the Delta 
assuming a rigid spacecraft at 17 to 23 cps If the 
spacecraft and adapter cannot be assumed rigid in this 
frequency band, the load factor must be increased to 
account for the additional load generated by dynamic 
amplification over the length of the spacecraft and 
spacecraft interstage adapter 

2) Longitudinally forward 3 Og combined vectonally with 
2 Og in either direction along any lateral axis 
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SPACECRAFT FLIGHT ACCEPTANCE 
DELTA ACOUSTIC NOISE 


Octave Band Center 
Frequency (Hz) 

Sound Pressure Level 
(db, ref 0 0002) 

Exterior of Shroud Interior of Shroud 

Duration 

15 8 

135 

123 

1 minute 

31 5 

136 

123 


63 

140 

125 

for 

125 

142 

130 


250 

145 

133 

comp! ete 

500 

145 

134 


1000 

143 

134 

exposure 

2000 

140 

131 


4000 

136 

127 


8000 

131 

121 


Overall 

151 

140 
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e) Shock 

The shock environment is specified in Figure 4 
31345 Orbital Conditions 

The satellite shall be designed to withstand the following orbital 
conditions 

a) Solar Radiation The intensity of solar radiation outside 
the atmosphere of the Earths mean distance from the sun 

shall be considered to have a mean intensity of 429 BTU/hr ft 2 

b) Earth Albedo The ratio of the total luminous flux incident 
on the earth in all directions to the total flux incident 

on the earth in a beam of parallel light shall be considered 
to be 0 34 

c) Earths Thermal Radiation The intensity of earth radiation 
outside the Earths atmosphere shall be considered to have 

a mean intensity of 68 BTU/hr- ft 2 

d) Earths Shadow The satellites emergence from the Earth shadow 
shall be considered between 29 to 35 minutes from the time the 
satellite entered the shadow 

3135 Transportability 

The ERTS satellite shall be designed to be transported by the ERTS 
transporter Each component of the GDHS shall be designed to be transported 
by common earner with a minimum of protection Special handling or 
packaging techniques shall be used to assure that transportation methods 
do not impose design penalties 

3136 Human Performance 

The design of all ERTS systems requiring man/machine interfaces shall use 
MIL-STD-803A, Part I, as a general guide 
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Figure 4 Delta Shock Spectra Input Spacecraft Flight Acceptance 
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3137 Safety 

31371 Dangerous Materials and Components 

Personnel functions within the hazardous area, following reaction control 
servicing and propulsion separation, or destruct system arming, shall be 
minimized 

Design provisions and procedures for the safe handling, storage, and 
installation of propellants, explosives and exciter devices shall 
conform to the safety requirements as indicated on the applicable 
assembly and installation drawings 

3 1 3 7 2 Noise and Vibration 

The noise and vibration levels associated with the ERTS and its components 
in required combinations as specified in TRW Systems Specifications 
D-13353 and D-l 3354 shall be controlled under its operating conditions 
to levels of tolerance to personnel as specified in MIL-STD-803A, Part I 

3 2 System Design and Construction Standards 

321 General Design and Construction Requirements 

3 2 11 Selection of Specifications and Standards 

Selection of specifications and standards for identification, control, and 
procurement of parts, materials, and processes not specified herein shall 
be in accordance with the provisions of MIL-STD-143A 

3212 Materials, Parts, and Processes 

Only materials, parts, and processes conforming to the applicable - 
requirements of the GSFC Preferred Parts List PPL-10 and of Group I and V 
documents of MIL-STD-143A shall be considered standard and shall be used 
Engineering drawings and associated lists shall be in accordance with the 
provisions of MIL-D-1000 and MIL-STD-100 
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3213 Standard and Qualified Parts 

32131 Electronic Parts 

Electronic parts shall be as specified in the contractor's ERTS Reliability 
Program Plan 

32132 Screening Tests 

Screening tests shall be as specified in the contractor's ERTS Reliability 
Program Plan 

32133 Derating 

Derating shall be as specified in the contractor's ERTS Reliability Program 
Plan 
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3 2 14 Fungus and Moisture Resistance 

Materials that are not nutrients for fungus shall be used whenever possible 
Where the use of fungus nutrients cannot be avoided, treating, packing, or 
other protective means shall be employed to ensure no degradation in 
system performance 

3215 Corrosion of Metal Parts 

System parts, including spares, shall be protected against corrosion 
To avoid electrolytic corrosion, dissimilar metals, as defined in MS-33586A, 
shall not be used in direct contact Protective methods and materials for 
cleaning, surface treatment, and application of finishes and protective coating 
shall be accomplished in accordance with MIL-F-7179C, where applicable 

3216 Interchangeability and Replaceability 

System designs shall meet the requirements of MIL-STD-100 for interchange- 
ability and replaceability 

3217 Workmanship 

32171 Workmanship Standards 

Workmanship shall conform to the requirements of the applicable process 
specifications relanng to fabrication and assembly as invoked by the 
particular assembly drawing Critical steps of fabrication that are item- 
peculiar shall be detailed in drawing notes which shall include appropriate 
criteria of workmanship Workmanship relating to all other aspects of 
fabrication, general handling, and storage shall be deemed adequately 
covered by the quality control program 

32172 Personnel Certification 

Personnel involved in assembly, soldering, welding, or other activity 
requiring special technical skills shall be certified as to their capability 
to perform such duties effectively as specified in the contractor's 
Quality Program Plan 

3218 Electromagnetic Interference 

The satellites shall be designed to minimize or eliminate inductive or 
electromagnetic RF interference effects 
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3219 Identification and Marking 

32191 Article Identification 

System components and ground support equipment components shall be 
identified in accordance with PR 12-1 

32192 Pipelines 

Fluid and gas tubing pipelines shall be identified and marked in accordance 
with MIL-STD-1247A, except for lines which do not have adequate space for 
such marking 

32193 Explosives 

All propellant grains, ignitors, squibs, and ordnance charges shall be 
classified and marked in accordance with applicable TRW procedures m 
compliance with appropriate government requirements 

3 2 1 10 Storage 

The satellites and associated ground support equipment, with the exception 
of batteries, shall be designed to be stored for a period of one year 
without requiring major refurbishment, maintenance or retesting at the end 
of storage 

322 Design Disciplines 
3 2 2 1 Civil 

Facility design shall be in accordance with joint industry conference 
standards and all other civil, architectural, and structural industry 
standards as applicable 

3 3 Perf orman ce All ocati ons 

3 3 1 Ground Data Handling System 

The Ground Data Handling System shall consist of two major subsystems as 
described in TRW Systems specifications D-13750 and D-13800, the Operations 
Control Center (OCC) and the NASA Data Processing Facility (NDPF) 
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3 3 11 Allocated Performance and Design Requirements 

33111 Mission Data Handling 

The Ground Data Handling System shall accept, record, process, and display 
PCM telemetry data and provide such data to users in compliance with the 
requirements of paragraphs 31111, 3111 2, 3111 3, and 31121 

33112 Command and Control 

The GDHS shall generate, translocate and validate spacecraft conmands in 
compliance with the requirements of paragraph 31129 

33113 Satellite Housekeeping Telemetry 

The GDHS shall accept and process spacecraft housekeeping telemetry data 
m compliance with the requirements of paragraph 3 1 1 2 12 

33114 Maintainability 

The GDHS maintainability requirements shall be as specified in the contrac- 
tor's GDHS Maintainability Plan 

3312 Peculiar Performance and Design Requirements 
33121 GDHS Software 

The software requirements of paragraphs 311215,311216, 31 12 17 
and 3 1 1 2 18 shall apply in their entirety to the Ground Data Handling 
System 

331121 Command Generation 

Command generation shall be provided to the appropriate STADAN & MSFN 
ground stations for translocation and uplink, transmission to the satellite 
and payload 

3 3 112 11 Command Modes 

a Real-time command modes shall be provided as defined in 
TRW Systems Specification No D-13750 

b Stored command modes shall be provided as defined in TRW 
Systems Specification No D-13750 
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331122 Command Translocation 

Commands shall be transmitted to the appropriate STADAN/MSFN ground station 
in a 600 bit message block via the NASCOM HSD lines by ADPE 

331123 Command Validation 

Command validation shall be provided by 

a Visual display prior to translocation 

b ADPE comparison with the command library 

c Automated MSFN error coding 

33122 Data Processing 

The GDHS shall have the capability to automatically process data as 
requested by the user Data shall be transformed into user oriented media 
and retained for use by the various users as specified in paragraph 
3 1 1 2 10 

33123 Unified Display Consoles 

The GDHS shall have unified display consoles in both the OCC and the NDPF 
The displays shall be capable of interacting with either the OCC or the 
NDPF data bases 

3313 Functional Interfaces 

33131 Ground Data Handling System/Remote Tracking Stations System 

The GDHS shall originate commands for spacecraft control and relay these 
commands to the remote tracking stations for transmission to the satellite 
The GDHS shall receive mission data and housekeeping telemetry from the 
remote tracking stations 

33132 GDHS/Users 

The GDHS shall accept requests for data from the users and shall process 
and disseminate data in accordance with user requests as defined in TRW 
Systems Specification No D-13800 
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332 Earth Resources Technology Satellite System 

The Earth Resources Technology Satellite (ERTS) System shall consist of 
spacecraft structure (including launch vehicle adapter and payload mounts), 
communication and data handling (including antennas), active thermal 
control, electrical integration, orbit adjustment, stabilization and con- 
trol, and solar array power supply subsystems and payload experiment pro- 
visions The satellite c vstem shall also include the necessary ground sup- 
port equipment to support the satellite 

3321 Allocated Performance and Design Requirements 

33211 Primary Mission Requirements 

The satellite system, in compliance with the requirements of paragraph 
31111,31112, and 31113, shall acquire data through the utili- 
zation of the payload The payload data shall be managed by the Ground 
Data Handling Systems Operations Control Center and shall be processed by 
the NASA Data Processing Facility The data shall be subject to bulk and/or 
precision processing modes as defined in TRW Systems Specification No 
D- 13800 

33212 Satellite Weight 

The satellite weights, in compliance with the requirements of paragraph 
3 1 1 2 11 and 3 1 1 2 13, for ERTS A shall be 1491 lbs , ERTS B shall 
be 1506 lbs 

33213 Payload 

The satellite shall contain provisions for a payload as follows in compli- 
ance with the requirements of paragraph 31122 and 31123 

Sensors 

a) Return Beam Vidicon Camera 

b) Multi -Spectral Scanner 

c) Data Collection Equipment 

d) Wide Band Video Tape Recorder 
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Weight 

The satellite shall be capable of supporting a payload weight including 
harness and adapter structure of 475 pounds 

Vol ume 

13 cubic feet mini mum 

33214 Orbit Injection 

The satellite shall be capable of removing the launch vehicle injection 
errors in compliance with the requirements of paragraph 31125 through 
an orbit adjust and Stabilization and Control function incorporating 
thrusters utilizing Krypton gas and reaction wheels as torque generators 
In plane and cross plane adjustment shall be provided yielding a capability 
of greater than 50 feet per second 

33215 Acquisiti on 

The acquisition of the satellite reference axis shall be initiated at a 
predetermined time subsequent to separation from the launch vehicle and 
at any time when one or more reference axes are lost 

3 3 2 1 5 1 Acquisition Sequence 

After separation from the launch vehicle the acquisition of the sun, 
earth, and orbit plane shall occur in the following sequence 

a) Sun Acquisition Sun acquisition in the sun-acquisition mode 
shall be achieved as early as possible, enabling the positive 
pitch axis of the satellite to be aligned to the sun 

b) Earth Acquisition Earth acquisition in the earth-acquisition 
mode shall occur by rotation of the satellite about the pitch 
axis, for the necessary length of time in orbit, until the yaw 
axis of the observatory sufficiently intersects the earth at 
which time the earth trackers will lock-on the earth 
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c) Orbit Plane Acquisition Orbit plane acquisition shall occur 
on manual command In this mode, the pitch axis shall be 
rotated from sun pointing to a position where the X axis is in 
the orbit plane with the -X axis nominally aligned with the 
satellite orbital velocity This shall be accomplished by 
means of a gyro compass control mode 

33216 Satellite Housekeeping Telemetry 

The satellite shall be capable of processing, storing and transmitting 
spacecraft and payload status information in recognition of the require- 
ments of paragraph 3 1 1 2 12 The downlink bit error rate shall be no 
greater than 1 x 10" 6 bps 

The housekeeping telemetry bit rate shall be 1 Kbs and 32 Kbs 
3 3 2 1 7 Communication and Command 

The satellite shall provide a coherent S-band transponder for ranging 

and tracking by the MSFN stations, in comoliance with the requirements 

of paragraphs 31128 and 31129 VHF equipment shall provide PCM 

data and a carner'for use by the STADAN network The satellite shall be 

capable of receiving, decoding and executing real time ground and stored 

-fi 

commands The uplink error rate shall be no greater than 1 x 10 bps 
Mission data shall be transmitted over two S-band downlinks The MSS 
wideband data link shall have a bit rate of 15 mbs 

3 3 2 1 8 Rel i abi 1 1 ty 

The spacecraft and payload shall have a design goal of one year, minimum, 
of stabilized orbital flight 

33219 Maintainability 

The satellite and GSE shall be designed in accordance with the requirements 
of paragraph 3132 


E-41 



Specification No D-13600 


3 3 2 1 10 Environments 

The satellite shall withstand the environments and operate in orbit 
as specified in paragraph 3134 

3 3 2 2 Peculiar Performance and Design Requirements 

33221 Orbit Plane Orientation 

The satellite shall be 3 axes stabilized and shall be aligned with the 
plane of the orbit within 0 7 degree during normal operation 

33222 Satellite Orientation Rate Errors 

The satellite orientation rate errors shall not exceed those specified 
in paragraph 3 1 1 2 14 

33223 Global Coverage and Image Overlap 

The requirements of paragraphs 3 112 6 and 3 112 7 shall apply in 
their entirety 

3 3 2 2 4 Separation Velocity 

The satellite separation equipment shall impart a velocity to the 
satellite relative to the final stage of the Thor-Delta of 4 2 ft/sec 
+ 10 % 
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33225 Electrical Power 

Primary power for the satellite shall be supplied by an 
onentable array of solar cells mounted on two paddles The array shall 
be supplemented during periods of heavy load and eclipse by two recharge- 
able nickel -cadmium batteries In addition to the solar array and 
batteries the electrical power function shall include charge control 
devices, dc-dc converters, a payload converter and pulse power and 
standby batteries The power subsystem shall be capable of delivering 
a total of 30,000 watt-mi n of energy per orbit to the satellite and 
payload This capability shall be based on a 103 minute orbit with a 
34 minute eclipse time at the end of a one year flight The satellite 
shall require 16,000 watt-mm of energy per orbit, the remainder shall 
be available to the payload and permit 20 minute sensor operation per 
orbit 

3 3 2 2 5 1 Power Profile 

The power required from the bus shall vary in time and an arbitrary 
schedule of sensor operation to produce the power profile is shown in 
Table III 

33226 Satellite Ground Support Equipment 

The satellite ground support equipment shall include equipment 
and instrumentation necessary for the verification and acceptance of the 
satellite for flight Equipment and instrumentation shall include 
capabilities for installing, aligning, and calibrating subsystems, and 
for verifying the electrical and functional integrity of components, 
subassemblies, and interrelated or cooperative end systems Equipment 
shall be provided to handle, service, and transport the spacecraft and 
satellite as well as to prevent malfunctions caused by dust, moisture, 
frost, airborne or applied chemicals, or other potentially degrading 
environmental factors 


E-43 



-44 


Table III Time History of Bus Current 


Time 

0-10 

10-20 

20-30 

30-80 

80-90 

90-103 

SUN/ECLIPSE 

SUN 




ECLIPSE 

34 mm 

Base Load, Amperes 

5 29 

5 29 

5 29 

5 29 

5 29 

5 29 

Data Collection System 

03 

03 

03 

03 

03 

03 

RBV, 155 W 

-» 

5 65 

5 65 

-- 


— 

MSS, 60W 

— 

2 15 

2 15 


— 

— 

W B Recorders, Z Record 


— 

6 0 




W B Recorders, Z P B 


— 

__ 

— 

5 4 

— 

W B Transmitter, 72W 


2 8 

— 


2 8 

““ 

-24+1/2 V Conv Losses 15°/ 

- 

1 5 

1 9 


1 0 


TOTAL, Amperes 

5 32 

17 42 

21 02 

5 32 

14 52 

5 32 
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3 3 2 3 Functional Interface 
3 3 2 3 1 ERTS System/Data Collection System 

The satellite shall receive data from the data collection platforms 
when the satellite is within the line of sight of the platforms The 
satellite shall furnish power to operate the airborne DCS equipment and 
shall transmit DCS data to the remote tracking stations for transmission 
to the GDHS 

3 3 2 3 2 ERTS System/Launch Vehicle System 

Detailed interface requirements between the satellite and the launch 
vehicle shall be as specified in TRW Specification D-13503 

332321 Satellite/Launch Vehicle Adapter 

The satellite to launch vehicle adapter shall be designed to align and 
attach to the fitting as described in paragraph 3 3 4 2 2 

3 3 2 3 2 2 Volume Constraints 

The satellite shall be designed to fit within the payload fairing 
envelope described in paragraph 3 3 4 2 1 
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332323 Separation 

The satellite separation signal and electrical power shall be furnished 
by the launch vehicle for separation from the launch vehicle 

33233 ERTS System/Remote Tracking Stations 

The satellite shall be capable of receiving ground commands from the 
remote tracking stations The satellite shall transmit data from the 
airborne sensors and data collection platforms to the ground stations for 
retransmission to the GDHS 

333 Data Collection System 

The data collection system (Figure 5) shall consist of data collection 
platforms which include sensors on the ground to collect data and a transmitter 
capable of transmitting the data to the satellite and equipment on the 
satellite to receive and transmit DCS signals The system shall operate 
in a random mode to eliminate the need for timing signals and interrogation 
commands and shall allow for a finite number of message collisions 

3331 Allocated Performance and Design Requirements 
33311 Mission Data 

The data collection system shall collect data in support of the primary 
ERTS missions m compliance with the requirements of paragraph 31111, 

3 1 1 1 2, and 3 1 1 1 3 

3 3 3 1 2 Ai rborne Equi pment 

The DCS shall contain airborne receiving equipment in compliance with 
the requirements of paragraph 31122 The DCS flight receiver 
shall be a double frequency conversion superhetrodyne design to operate 
on a frequency of 400 to 406 MHz with a noise figure limited to two db 
maximum 

3332 Peculiar Performance and Design Requirements 

3333 Functional Interfaces 

33331 Data Collection Platforms/satellite 

The data collection platforms shall be capable of inputting data to the 
satellite DCS receiver on a 117 KHz bandwidth when the satellite is with- 
in the line of sight of the platforms 
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The data shall be converted to an IF frequency and routed to a subcarrier 
oscillator in the spacecraft The data collection airborne receiver 
equipment shall receive the necessary power to operate from the spacecraft 


E-47 



E-48 


D 11590 



Figure 5 


DCS Specification Tree 













Specification No D-13600 


334 Launch Vehicle System 

The launch vehicle system shall consist of a Thor-Delta launch vehicle, 
the launch pad, and other launch support facilities at the Western Test 
Range (WTR) 

3 3 4 1 Allocated Performance and Design Requirements 
3 3 4 1 1 Payload Capability 

The launch vehicle payload capability shall be no less than 1520 lbs 
in compliance with the requirements of paragraph 31 1 23, 31 124, 
and 3 1 1 2 13 

3 3 4 1 2 Orbit Injection Errors 

The launch vehicle shall place the satellite in orbit within the following 
nominal parameters in compliance with the requirements of paragraph 
3 112 5 


a) Inclination of orbit 

b) Equatorial Altitude 

c) Right ascension of ascending 
node at vernal equinox 

d) Eccentricity 

e) Eclipse duration 

f) Period 

33413 Separation Rate Errors 


99 08° + 0 05° 

492 +10 nautical miles 
142 5 minus zero plus 7 5° 

less than 0 006 
29 to 34 minutes 
103 3 + 0 3 minutes 


The launch vehicle rate errors at satellite separation shall be limited 
to the following 3a values 

Pitch and Yaw 0 1 degree/sec 

Roll 05 degree/sec 


3 3 4 1 4 Ascent Loads 

j t 

The launch vehicle shall support the satellite during lift-off, ascent, 
and insertion into orbit 
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3 3 4 2 Peculiar Performance and Design Requirements 

The launch vehicle system shall provide the launch facilities and support 
functions required during all prelaunch and launch activities 

3 3 4 2 1 Payload Fairing Envelope 

The payload fairing envelope shall be as described in Figure 5 
33422 Alignment and Attachment 

The launch vehicle shall provide alignment and attachment provisions for 
the satellite as shown in Figure 6 

3 3 4 2 3 Satellite Separation 

The launch vehicle system shall provide the signal and power required to 
initiate satellite separation when the orbit is achieved 

3343 Functional Interfaces 

33431 Launch Vehicle System/ERTS System 

The detailed interfaces between the launch vehicle and the satellite shall 
be as specified in TRW Specification D- 1 3503 The launch vehicle system 
shall provide facilities and equipment necessary for the operation of the 
satellite ground support equipment 

335 Remote Tracking Stations System 

The Remote Tracking Stations System shall consist of the Gilmore (Alaska), 
Corpus Christ! (Texas), Rosman (North Carolina), and NTTF (GSFC/Maryland) 
stations dedicated to the ERTS program and other stations accessible 
through the NASCOM network for backup and launch and ascent support 

3351 Allocated Performance and Design Requirements 

33511 Orbit Determination 

The ground stations shall track and acquire range and range rate data to 
permit precision orbit determination in compliance with the requirements of 
paragraph 31125 and 3 1.1 28 

33512 Satellite Telemetry 

The tracking stations shall acquire satellite housekeeping data in 
compliance with the requirements of paragraph 3 1 1 2 12 
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NOTES: 
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33513 Satellite Command, Control and Data Acquisition 

The ground stations shall transmit commands to the satellite and acquire 
mission data form the satellite in compliance with the requirements of 
paragraphs 31111.31112, 31113 and 3 112 9 

3352 Peculiar Performance and Design Requirements 
3 3 5 2 1 Commanding and Control 

The Remote Tracking Station shall transmit commands required to operate 
the satellite 

33522 Satellite Telemetry Data Receiving 

The Remote Tracking Station shall be capable of receiving status and 
mission data telemetered from the satellite 

3353 Functional Interfaces 

33531 Remote Tracking Station System/GDHS 

The Remote Tracking Stations shall receive commands originated in the GDHS 
for transmission to the satellite The remote tracking stations shall relay 
mission and housekeeping telemetry received from the satellite to the GDHS 
The availability of the long lines for transmission of data from remote 
tracking stations to the GDHS shall be at least 0 995 

33532 Remote Tracking Station System/ERTS System 

The Remote Tracking Station shall transmit commands to the satellite and 
shall receive mission and housekeeping data from the satellite 

3 3 6 Users 

The users shall consist of the Department of Agriculture, Department of 
Interior, Department of Navy, Department of Commerce and other governmental 
and private agencies concerned with the earth's natural and cultural 
resources The users shall request data from the GDHS and furnish inputs 
to the GDHS for the format and types of data desired 
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4 PERFORMANCE ASSURANCE 

The contractor shall develop a performance assurance program in compliance 
with the provisions of NASA document NPC-250-1 , NHB 5300 4 (IB and 3A) 
as further defined by TRW documents Nos PAR 700-52, PAR 700-53, PAR 700-54, 
PAR 700-55, PAR 700-56, PAR 700-61 and the applicable TRW plans as shown in 
Section 2 The objectives of this program shall include 

a Configuration Management in compliance with GMI 8040 1 
b Verification of the total program performance 
c Verification of individual system's performance 
d Verification and accumulation of natural and cultural resource 
data for all missions 

e Demonstration of ERTS System performance repeatability and overall 
program operability and reliability 
f Isolation of failure modes and evaluation of resulting 
corrective action 

g Verification of ground data handling and data collection system 
performance 

Verification of performance, design, and construction of all program 
elements shall be accomplished by physical inspection, review of analytical 
data, demonstrations, tests, and reviews of test data Compliance with 

# 

program requirements which cannot be validated by inspection, 
demonstration or test shall be validated by engineering analysis 

The test program shall be made up of Phase I and Phase II tests as defined 
in Section 6 

4 1 Phase I Program Test Requirements 
4 1 1 Engineering Test and Evaluation 

Integrated program tests shall be conducted on the developmental model 
systems in direct support of design and development activity Results 
of component, equipment, and subsystem tests performed on this and other 
programs shall be used to support the development test program 


E-54 



Specification No D-13600 


4 1 2 Qualification Testing 

Formal qualification testing of the program elements shall be accomplished 
at the environment levels required to assure that the systems will comply 
with the requirements of paragraph 3134 herein 

Formal qualification testing of the ERTS satellite and the ground 
data handling system shall be accomplished at the worst probable combination 
of environmental levels specified in paragraph 3134 Formal qualifi- 
cation testing of ground support equipment shall be accomplished at 
applicable environmental levels 

413 Electrical and EMI Tests 

Electrical and EMI tests shall be conducted to validate the performance 
of the program elements to satisfy the requirements of paragraph 3218 
when operating in a simulated system configuration Testing shall be 
accomplished by operating all electrical equipments in normal prelaunch 
and launch sequence except when more severe environments are required 
to establish operation with adequate safety margins 

4 1 4 Reliability 

Tests performed for the purpose of obtaining reliability and maintain- 
ability confidence shall be performed on components as approved by the 
procuring agency 

415 Acceptance Testing 

Formal acceptance testing of the program elements shall be accomplished at 
the environmental levels specified in the contractor's environmental test 
specifications, procedures and plans 

4 2 Phase II Program Test Requirements 

The Phase II test program shall validate the performance and operability 
of the system The validation shall be accomplished through the analysis 
of test data acquired from integrated system ground tests and flight 
tests 
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4 2 1 integrated System Ground Tests 

These tests shall validate the compatibility and performance of the 
integrated program airborne and ground equipment prior to, and for the 
purpose of, flight test initiation The tests to be performed shall 
include but not be limited to the following 

a Compatibility verification of physical and functional interface 
between systems 

b Verification of system performance in the presence of the EMI 
environment 

422 Flight Test 

The flight test program shall provide flight test data to verify that 
the performance requirements of the system have been satisfied The 
flight test data and analysis of the data shall validate the airborne 
equipment performance requirements of Section 3 of this specification 

4 2 3 Failure Criteria 

The inability of the flight test program to demonstrate (in accordance 
with the test capabilities available) the performance requirements of the 
program shall be defined as a system failure Individual system 
equipment failures which may cause system failure include but are not 
limited to the following 

a Failure of the launch vehicle to insert the satellite into 
orbit 

b Failure of the satellite to acquire data 
c Failure to demonstrate system operability 
d Failure of the remote tracking stations to support the 
satellite or failure to acquire satellite data 
e Failure of the ground data handling system to acquire 
and process satellite data or to display information 
required by users 
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5 PREPARATION FOR DELIVERY 
5 1 Packaging 

Packaging of each item for delivery to the procuring agency shall be 
as specified below 

511 Containers 

Individual containers shall be so constructed as to allow removal of the 
article for inspection without destruction of the container or of the 
wrappers affixed to the article If paper wrapping is used on the article, 
acid free paper shall be used As an objective, the container shall 
provide equal protection, without use of special tools, to articles 
repackaged following inspection 

512 Special Instructions 

If the article requires special attention during receiving, inspection, 
installation and operation or if non-obvious characteristics require 
that special handling be used, the procuring agency shall be notified 
under separate cover and a removable instruction tag shall be attached 
Attachment shall be to the shipping container or to the article as 
appropriate 

5 2 Marking 

5 2 1 Marking for Shipment 

Exterior shipping containers and non-carrier packages and separately 
shipped items shall be marked in accordance with MIL-STD-129D 

522 Re- Inspect! on 

Articles requiring periodic re-inspection shall be marked with the next 
inspection date 
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6 NOTES 
6 1 Defimti ons 
Phase I Program Test 

Development test and evaluation of individual components, assemblies, 
and subsystems and in certain cases, the complete system, which is 
predominantly conducted by the contractor 

Phase II Program Test 

Testing and evaluation spanning the integration of configured 
items into a completed system in as near an operational configuration 
as is practicable, which is a joint contractor - NASA effort under 
NASA control 

6 2 Abbreviations 

BTE - Bench Test Equipment 

ERTS - Earth Resources Technology Satellite 

DCS - Data Collection System 

GDHS - Ground Data Hnadling System 

G$E - Ground Support Equipment 

GSFC - Goddard Space Flight Center 

MSFN - Manned Space Flight Network 

MSS - Multi -Spectral -Spectrometer 

MTBF - Mean Time Before Failure 

MTTR - Mean Time To Repair 

NDPF - NASA Data Processing Facility 

OCC - Operations Control Center 

RBV - Return Beam Vidicon 

STADAN - Satellite Tracking and Data Acquisition Network 
VHF - Very high frequency 
WTR - Western Test Range 
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